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PRÉFACE 


Le rôle pédagogique des techniques du calcul est reconsidéré 
dans les écoles techniques supérieures d'U.R.S.S. D'objet d'étude, 
le calculateur se transforme en outil permettant une assimilation 
efficace des processus et phénomènes de la simulation mathéma- 
tique. La matérialisation de cette nouvelle fonction « pédagogique » 
des calculateurs implique l'enseignement précoce de la programma- 
tion dès le second ou le troisième semestre et ne se justifie que si l’on 
dispose d'une littérature didactique adaptée au niveau des étudiants 
du second cycle. Mais la situation est loin d’être satisfaisante, 
notamment en ce qui concerne la simulation mathématique sur 
calculateurs analogiques. 

Les manuels destinés à l'initiation des élèves de première et 
deuxième année aux calculateurs analogiques faisant défaut, le 
Cours pratique en question tente de remédier à cette carence. Par 
son contenu, sa méthode et le caractère de l'exposé, ce cours poursuit 
le même objectif que les « Eléments de simulation sur calculateurs 
analogiques » du même auteur, savoir familiariser les élèves de 
première et seconde année avec la programmation et la manipulation 
des calculateurs analogiques dans une perspective d'application 
à l'étude de disciplines scientifiques générales. 

Ce cours généralise les acquis de dix ans de travaux de l'Institut 
de l’Acier et des alliages sur la mise au point d'un système d'utili- 
sation permanente des calculateurs à tous les degrés de l’enseigne- 
ment. Nous exprimons l'espoir que l'édition de cet ouvrage fera 
profiter les instituts similaires du fruit de notre travail. 


S. Emélianov, 


mombre correspondant de 
l'Académie des sciences de l’U.R.S.S. 


AVANT-PROPOS 


Le cours pratique est articulé sur quatre chapitres. 

Le premier est introductif. I] renferme un grand nombre d'exer- 
cices sur l'établissement d'une correspondance entre les schémas 
fonctionnels, les circuits électriques et les expressions mathématiques 
réalisées. Ces exercices ont pour but de faire assimiler à l'étudiant 
les symboles et le langage de la programmation des calculateurs 
analogiques (langage des organigrammes). 

Le second chapitre est consacré aux problèmes de la programma- 
tion des calculateurs analogiques, relatifs au câblage d'équations 
différentielles ordinaires données sous forme de problèmes de Cauchy. 
Les problèmes abordés se rapprochent au maximum du niveau des 
premiers semestres et portent sur l'étude, avec la machine, de fonc- 
tions et de courbes connues des cours d’analyse mathématique, de 
géométrie analytique et de géométrie descriptive. 

Le troisième chapitre se fixe deux objectifs: d'abord rompre 
l'étudiant au maniement d'un calculateur analogique, en l'occur- 
rence le MH-7M, et ensuite développer les principales méthodes de 
simulation mathématique de processus physiques décrits par des 
problèmes de Cauchy, des problèmes aux limites et des équations 
différentielles aux dérivées partielles. 

Le quatrième chapitre est appliqué. Il contient une vingtaine 
d'exercices de simulation mathématique empruntés à diverses dis- 
ciplines. 

Le diagramme suivant nous semble le plus favorable à l'assimi- 
lation du « Cours pratique »: 


7 Ch. 2, K$3à SNS 
Ch. 1 + Ch. 2, $1 + Ch. 3,1, 2—Ch.2,82 3Ch. 4 
NCh. 3,#3à5/ 
où Ch. désigne le chapitre. 


L'auteur tient à exprimer sa profonde gratitude à S. E mé- 
lianov, membre correspondant de l'Académie des sciences de 
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l'U.R.S.S., dont l'attention, l’intérêt et l'aide ont contribué à la 
formation et l'orientation des principales conceptions pédagogiques. 

L'auteur est également profondément reconnaissant aux profes- 
seurs B. Anissimov et L. Presnoukhine dont les 
critiques bienveillantes et les conseils ont considérablement influencé 
l'essence des « Eléments de simulation sur calculateurs analogiques » 
et du « Cours pratique de simulation ». 


CHAPITRE PREMIER 


OPÉRATEURS ET SCHÉMAS 


Dans les exercices des planches P1 à P29, on demande d'établir 
une correspondance entre les circuits électriques, les schémas fonc- 
tionnels et les expressions mathématiques. 

Pour faciliter la résolution des exercices, on a donné sur la 
figure 1.1 un tableau des principaux opérateurs avec les circuits 
électriques et les schémas fonctionnels correspondants, ainsi que les 
expressions mathématiques réalisées. 


ExEMPLE 1.1. Etablir une correspondance entre les circuits 
électriques des opérateurs de la planche P1 et les schémas fonction- 
nels des opérateurs de la planche P2. | 

En utilisant le tableau de la figure 1.1 et en calculant les coef- 
ficients de transfert des opérateurs, on trouve les équivalents fonc- 
tionnels des circuits électriques donnés. La solution du problème 
est représentée sous forme du diagramme de correspondance de la 
figure 1.2. Les six numéros supérieurs figurent les circuits de P1, 
les huit numéros inférieurs, les schémas de P2. L'équivalence est 
traduite par un segment reliant les numéros supérieurs aux infé- 
rieurs. À l'examen du diagramme on remarque qu'à deux reprises 
deux circuits électriques distincts (2; 4 et 1 ; 3) ont le même équi- 
valent fonctionnel (3 et 8). Cela veut dire que le passage d’un schéma 
fonctionnel à son équivalent électrique n'est pas univoque. 

1.1. Trouver sur P4 les équivalents électriques des schémas 
fonctionnels représentés sur P3. 

1.2. Trouver sur P6 les schémas fonctionnels et sur P7 les expres- 
sions analytiques correspondant aux circuits électriques représentés 
sur PS5. 

Indications. 1) On suppose que les potentiomètres des 
circuits de P5 portent tous le même coefficient de transfert «. 

2) Les tensions électriques u, et u, (les variables machine) des 
circuits électriques de P5 sont assimilées aux variables mathéma- 
tiques zx, et x, des schémas fonctionnels de P6. 

3) La solution de l'exercice 1.2 doit contenir un diagramme de 
correspondance composé de trois groupes de numéros reliés par des 


segments orientés de P5 vers P6, puis vers P7. D 
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Amplificateur opérationnel z ÿ 
y—=kKz 


Amplificateur d'échelle zx ce y 
y= —QAT 

Sommateur " Q, 

n 2 1 

2:00 

ne \ . ÿ 

y=— Sax SE 

ji j 


Intégrateur (sommateur intégrateur) 
t 


n 
y = — D ae j(t) dt+ y (0) 


0 j=1 
Potentiometro 5 y 
y =$z —@-— 
Multiplieur Tr Z 
= +} y T— 


Z 
Diviseur T D 
2: z/y 
ÿ 
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Ue Ua Non K>4.10t—-gain de l'amplificateur opé- 
rationnel 


Uy = —(Re.r/R)uzx, = Re.r./ R—coeffi- 
cient de transfert de l'amplificateur d'e- 
chelle 


n 
Re.r. 
My D Ux 
_ "4 j 
4j = Re.r./Rj—coefficient de transfert du 
sommateur pour la /-ème entrée 


{ n 
= —: Ù Uxj (£) ( 
um O=—) De d+ue (0. 
0 2-1 
&j—1/RjC—coefficient de transfert de 
l'intégrateur pour la j-ème entrée 


uy=(r/R)ux, 
B=—r/R—coefficient de transfert du po- 
tentiomètre, OSB<1 


u::=Uxu,,/100, 1/100—échelle constructi- 
ve du multiplieur 


uz—=10ux/u,. 10—eéchelle constructive du 
diviseur 


Fig. 1.1 
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Planche 1 
! ! 
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025 05 
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a! 05 0,05 
05 025 0025 


Planche 2 
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Planche 


Planche 4 


Planche 5 


Planche 6 
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&) Dans certains schémas fonctionnels de P6, l’amplificateur 
opérationnel (AO) est considéré comme un opérateur isolé. En écri- 
vant les expressions mathématiques délivrées par ces circuits, on ne 


Fig. 1.2 


perdra pas de vue l’importante propriété de l’amplificateur opéra- 
tionnel à contre-réaction, savoir la somme algébrique des variables 
d'entrée est nulle. 


7 


1 


7 
- (sr +0x2) 


- (x, +0,14x;) 


1 


Gal (x, +0,12) 


—(01x,+ax:) 


Planche 7 


Montrons sur un exemple simple comment on peut décrire le 
fonctionnement d'un circuit en se servant de cette propriété. 


ExemPLe 1.2. Trouver l'expression de la variable z,, délivrée 
par le schéma de la figure 1.3, a. 


2—0112 
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_ Ce schéma met en jeu un AO à contre-réaction par l’intermé- 
diaire d’un potentiomètre de coefficient de transfert f$. La propriété 
formulée entraîne 


Bit + Aloe + A3B32a =. 
En explicitant z,, on obtient 


__ GB1T1 + Goo 
abs 


5) Certains schémas fonctionnels de P6 renferment un sommateur 


Za = 


L, æ, 
Z4 

T, œ, 

y 


Fig. 1.3 
à réaction. Montrons sur un exemple concret comment on trouve 
l'expression mathématique délivrée par un tel sommateur. 


EXEMPLE 1.3. Trouver l’expression de la variable z, (cf. fig. 1.3, b). 
Le schéma de la figure 1.3, b étant celui d’un sommateur, on a 


2p = —(@Birs + Gore + G3B328). En explicitant z,, on obtient 
2 AB1T1 + Lee 
co 1+ ss 


Le schéma envisagé câble, comme celui de l'exemple 1.2, le 
produit d’une somme par un nombre. Dans ce dernier exemple, ce 
nombre est une fraction propre. Dans l'exemple 1.2, le schéma 
effectue un produit par le nombre positif 1/œ&:B3. 

1.3. Trouver sur P9 les circuits électriques et sur P10 les expres- 
sions analytiques correspondant aux schémas fonctionnels de P8 
(cf. indications pour l'exercice 1.2). 

1.4. Etablir une correspondance entre les expressions analytiques 
de P11, les schémas fonctionnels de P12 et les circuits électriques 
de P13 (cf. indications pour l'exercice 1.2). 

1.5. Les expressions analytiques de P14 sont réalisées par trois 
procédés différents (par des sommateurs sur P15, par des sommateurs 
à réaction sur P16 et par des AO à réaction sur P17). On demande 
de construire le diagramme de correspondance entre P14, P15, P16 
et P17. 


Planche 9 
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4 


1 1 1 
Gros 090223) —(0r,+222) oi ET +) 
Planche 10 
1 4 7 
—(2,5x,+ax,) | -(x,+10ax,) |- £ (2x, +x;) 


] 


Planche 11 


Indications. 1) Monter trois schémas fonctionnels câblant 
les expressions analytiques de P14 en utilisant un sommateur, un 
sommateur à réaction et un AO à réaction. 

2) Chercher les schémas analogues dans P15, P16 et P17 et cons- 
truire le diagramme de correspondance. 

3) Choisir entre les schémas fonctionnels réalisant une même 
expression analytique celui qui semble le plus favorable. 

1.6. Trouver les équations différentielles des schémas fonction- 
nels de P18. Construire le diagramme de correspondance entre P18, 


Planche 12 


Planche 13 


22 


OPÉRATEURS ET SCHÊÉMAS 


2 N] 
-Éz+5æ) 


Planche 14 


[CH. 1 


Planche 15 


Planche 17 
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J 
y"+G7g"-08y-4c=0 y"+01y -3y -c=0 y”+0,2y-6y-2 c=0 


Planche 19 


P19 et P20. Les équations différentielles sont données dans P19, 
la condition initiale y” (0) et les constantes a, b et c, dans P20. 

Indications. Déduisons l'équation différentielle câblée 
par le schéma fonctionnel donné. 
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02a 


Planche 20 


EXEMPLE 1.4. Trouver l'équation différentielle câblée par le 
schéma fonctionnel de la figure 1.4, a. Ce schéma comporte le som- 
mateur intégrateur [1], sérié à l’intégrateur [2]. Cette chaîne d'inté- 
grateurs est embrassée par une réaction par l'intermédiaire de l’in- 
verseur [3] et d’un potentiomètre de coefficient de transfert B, 


Fig. 1.4 


De plus, dans le sommateur intégrateur [1] est réalisée une réaction 
par l'intermédiaire d’un potentiomètre de coefficient de transfert 
Ps 
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Pour trouver l'équation différentielle cherchée on remontera la 
chaîne d'’intégrateurs, i.e. on ira de la sortie de l'intégrateur [2], 
qui débite la variable y ({), aux entrées du sommateur intégrateur 
[1] (fig. 1.4, b). L'opérateur [2] est un intégrateur qui ne délivrera 
la variable y (t) que s’il reçoit la variable —y” (t)/&,, &,; étant son 
coefficient de transfert. Or la variable —y” ({)/x&, est élaborée à son 
tour à la sortie du potentiomètre de coefficient de transfert $;, 
y'(4) 
aps * 
laquelle provient du sommateur intégrateur [1]. On peut donc 
affirmer (en décomposant mentalement les opérations de sommation 
et d'intégration réalisées par le sommateur intégrateur) que dans la 
section verticale des entrées du sommateur intégrateur [1], on aura 


AUDE droite, et la somme 
ab 


donc à l'entrée de ce dernier est appliquée la variable — 


la variable 


Labs 


CŒ— y" (t) — moy (1) 


qui lui est égale à gauche, soit 


g" (0 . ŒnBa _ 
ah Ua, Ÿ (£) — Bay (4). 


Toutes réductions faites, on obtient l'équation différentielle 


y" (4) + ab” (1) + œomiBeBay (1) — maaBsc = 0. 


La détermination de la valeur initiale de la variable y (t) est aisée, 
car y (0) — b. Lorsqu'on cherchera la valeur initiale de y” ((t), 
on ne perdra pas de vue que le sommateur intégrateur [1] délivre 


la variable — CRC (et non —y’ (t)) qui pour { = 0 est égale à a, 


d'où y’ (0) = —a;B;a. Pour les schémas de la figure 1.4, on a en 
définitive 


y" (t) + œafay (1) + œonsBoBay (1) — maaBic = 0, 
y" (0) = —aifsa, y (0) = 6. 


1.7. On demande les équations différentielles câblées par les 
schémas fonctionnels de P21. Trouver la valeur initiale y’ (0) de 
la variable y’ (t). Construire le diagramme de correspondance à 
l’aide des schémas de P21, P22, P23 (cf. indications pour l’exer- 
cice 1.6). 

1.8. On demande les systèmes d’équations différentielles dont 
les solutions x (t) et y (t) sont réalisées par les schémas fonctionnels 
de P24. Construire le diagramme de correspondance entre P24 et 
P25. 
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Planche 21 


22 


ô 
y“ 02y*9y-0,2y=20|y"+ Oly'-15(y")$ 4y0 


J 6 9 
y0ly*03yy-0,3y-0|y#01y+02y-4y=0|y 0. ly+2y y+64 =0 


Planche 22 


1.9. Trouver les équations différentielles câblées par les schémas 
fonctionnels de P26. Construire le diagramme de correspondance 
entre P26 et P27. 

Indication. Se servir de la propriété de l’amplificateur 
opérationnel à contre-réaction pour composer l'équation différen- 
ticlle. 
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Planche 24 


1.10. On demande les équations différentielles cäblées par les 
schémas fonctionnels de P28. Construire le diagramme de corres- 
pondance entre P28, P19 et P18. 
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7 


ô 
z'+0,17 +2y=0 


5 
z'-9xy-0,29x=0 |r-0,2zy =0 
y -9zy+01z=0 |y-zy+0,3y=0 
9 


y“+2zy"-0,3y=0 
6 


—y"+y=0 
4 Gzy"-2y-02z=0|y"-yx"+3x =0 


3 
z-2xy"+0,1y=0 
y"+0,6zy"-0,2y=0 


Planche 25 


Planche 26 


29 
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27 


1 
J'Y 44) +270 | y y +07g + 3y-0 


Planche 28 
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1.11. Composer les équations différentielles câblées par les 
schémas fonctionnels de P29. Construire le diagramme de corres- 
pondance entre P29, P22, P21. 

1.12. La figure 1.5 représente trois couples de schémas fonction- 
nels donnant respectivement les solutions d'équations du second, 
du troisième et du quatrième ordre. Montrer que chaque couple de 


Planche 29 


schémas fonctionnels résout la même équation différentielle malgré 
la difference de structure. Trouver les conditions remplies par 
a, b, c, d'et a*, b*, c*, d* pour que chaque couple de schémas 
réalise les mêmes solutions particulières des équations différen- 
tielles. 

1.13. La figure 1.6 représente quatre schémas fonctionnels. 
Trouver les équations différentielles résolues ; résoudre ces équations 
analytiquement et étudier l'influence sur lim y; (t) (j = 1, ..., 4) 

t 


des erreurs £&,, £+, €3, €, Commises pendant l'initialisation des 
intégrateurs; montrer que lorsque e; = O (j = 1, ..., 4), on a 
yy (t) = exp (—t) (j = 1, ..., 4); interpréter les résultats obtenus 


(CH. ! 


OPÉRATEURS ET SCHEMAS 


Fig. 1.5 


CH. !] OPÉRATEURS ET SCHEMAS 33 


Fig. 1.6 


en faisant intervenir la notion de « stabilité des solutions des équa- 
tions différentielles ». 


1.14. Trouver l'équation différentielle réalisée par le schéma de 
la figure 1.7. 


Montrer en résolvant PAR EE l'équation obtenue que: 


1) si z(4) > 0, lim y (t) = TU indépendamment de la con- 


dition initiale y, et la vites sse de convergence croît avec p; 


Fig. 1.7 
9) si z(t) < 0, la quantité lin y (t) est infinie. 


Les multiplieurs de sombreus calculateurs analogiques sont 
construits sur la base du schéma de la figure 1.7. 
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CHAPITRE 2 


EXERCICES DE PROGRAMMATION 
DES CALCULATEURS ANALOGIQUES 


$ 1. Résolution d’équations différentielles ordinaires 


EXEMPLE 2.1. Monter un schéma fonctionnel donnant la solution 
de l'équation différentielle ordinaire linéaire à coefficients cons- 


tants 
y" (t) + ay" (€) + aoy (4) = f (+) (2.1) 
qui vérifie les conditions initiales 
y" (0) = y, y (0) — yo. 


Suivant la méthode générale *}), la mise au point d’un schéma fonc- 
tionnel résolvant cette équation différentielle comporte six étapes. 

Première étape. On résout l'équation (2.1) par rapport à la dérivée 
supérieure 


y" (€) = —ay" (t) — aoy (t) + f (t). (2.2) 


Deuxième étape. On représente une chaîne d’intégrateurs montés 
en série dont le nombre est égal à l’ordre de l'équation. Ici cette 
chaîne comporte deux intégrateurs (fig. 2.1, a). 

Troisième étape. On applique mentalement la dérivée supérieure 
y” (t) à l'entrée de la chaîne et on note de gauche à droite les varia- 
bles de sortie des intégrateurs (fig. 2.1, b). D'où que la variable 
y” (t) provienne (nous laisserons ce problème de côté pour l’instant), 
on aura la variable —y’ (t) à la sortie de l'intégrateur [1] et la va- 
riable y (t) à la sortie de l’intégrateur [2]. 

Quatrième étape. On construit le schéma fonctionnel qui réalise 
le second membre de l'équation (2.2) sous réserve qu'il existe un 
opérateur générant la fonction f (f), second membre de l’équation 
(2.1). Ce schéma est représenté en traits pleins sur la figure 2.1. c. 
Il comporte l'inverseur [1], le sommateur [2] et un générateur de la 
fonction donnée —f (t). Le sommateur [2] délivre le second membre 
de l'équation différentielle (2.2) 


—@ÿ" (4) — ay (t) + f (6). 


*) Voir A. Ourmaev, Eléments de simulation sur calculateurs analogi- 
ques, Mir, 1978, à partir du chapitre 5. 
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Cinquième étape. La sortie du sommateur [2] est appliquée à 
l'entrée de l’intégrateur [1] (fig. 2.1, d). Mathématiquement cela 
veut dire que le second membre de l’équation (2.2) est égal au premier. 

Sirième étape. On donne les conditions initiales aux intégrateurs 
en tenant compte des variables de sortie de ces derniers (fig. 2.1, e). 

Remarques. 1) Le schéma de la figure 2.1,e génère les 
variables y” (t), y’ (t), y (t) que l’on peut étudier avec les techniques 
du calcul analogique. Si la variable y” ({) ne présente pas d'intérêt 


Fig. 2.1 


pour le problème posé, le schéma 2.1, e peut être grandement sim- 
plifié par groupement des fonctions du sommateur [2] et de l'inté- 
grateur [1] en un sommateur intégrateur comme le montre la figu- 
re 2.1, f. Ceci permet de réaliser une économie d'opérateurs mais 
en revanche réduit l'information sur la solution, car la variable 
y” (t) est éliminée du schéma. 

2) Une fois que le schéma est construit, nous devons nous assurer 
que la machine reproduira correctement y (ti). Pour contrôler le 
fonctionnement de la machine on fait appel à des relations de con- 
trôle. Dans l'exemple envisagé, pour relation de contrôle on peut 


3* 
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prendre l'équation différentielle initiale sous la forme suivante: 
y" + ay + aoy — f (t) = 0. Le premier membre est une somme de 
quatre variables. Sa réalisation implique un seul sommateur à quatre 
entrées. Ces variables sont élaborées par le schéma de la figure 2.1, e 
et leur sommation avec le schéma de la figure 2.2 ne présente aucune 
difficulté. A la sortie du sommateur de contrôle est branché un 


D: (g"+08"+074 


-f(t)) 


Fig. 2.2 


voltmètre qui devrait indiquer le zéro si la machine fonctionne 
normalement, le schéma est correct et les coefficients de transfert 
exacts. 

2.1. Construire les schémas fonctionnels des équations diffé- 
rentielles linéaires représentées sur P19. Construire le diagramme de 
correspondance entre P19 et P28. 

2.2. Construire les schémas fonctionnels des équations diffé- 
rentielles non linéaires représentées sur P22. Construire le diagramme 
de correspondance entre P22 et P29. 

2.3. Construire les schémas fonctionnels des équations diffé- 
rentielles linéaires et des systèmes d'équations différentielles linéaires 
suivants: 


1) y" + 4y° + 3y = 0, 14) y" + 4y = 0. 

2) 2y" + Sy" + 2y = 0. 15) y" + 4y° + 5y = 0. 
3) y” + 2y° + 10y = 0. 16) yIV — y = 0. 

4) y” — 8y = 0. 17) yYI + y = 0. 

5) ylv — 8y = (. 18) 4y” + 4y° + y = 0. 

6) y” + 2y° + y = 0. 19) yY — 10y°" + 3y° = 0. 
7) yY + 6ylV + Gy°” = 0. 20) y” + 3y" + 3y° — y = 0. 
8) yIV + 8y" + y = 0. 21) y — 57" + 4y = 0. 
9) y" + y — y + y = 0. 22) y — 3y + 2y = 0. 
10) y + 4y° + 3y = 0. 23) x" — +3y" + 2y, 

11) yIY + 4y°" + Sy = 0. y" = —3x" — 2x. 

12) y" + y — 2y = 0. 24) x° + x — 8y = 0, 


43) y” + 2y" = 0. y —zxz+y=0. 
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25) z° — —zx + 3y, 


y = 3x — y. 
26) x’ = 2x + y, 
y" = 4y — x. 
27) x° = 2y — 3x, 
y" = y — 2x. 


28) z' = rz+z—7y, 
y = ZT+y—2z, 


z = 2x + y. 
29) z' = zx — y, 

y" = y — 4x. 
30) z'=zr + y, 

y" —= 3y — 2x. 


31) z° + x + 5y = 0, 
y —z—y—=t. 

32) zx = 3x — y, 

y" = 4x — y. 

33) z' — 5x — 3y = 0, 
y" + 3x + y = 0. 

34) x = x — 2y — 2, 
y =y—z+z, 


ZzZ = ZI — 2. 
35) z' = 3x — y + 2, 
y =z+y+z, 


2" = 4x — y + 4z. 
36) z' = 2x — y + 2, 

y = x + 2y — 2, 
Z = Zz— y + 2z. 
37) zx° = 4y — 2z — 3x, 


y =2z+2, 
z' = 6x — 6y + 5z. 
38) z' = 2z + y, 


y = x + 3y — 2, 


z' = 2y + 3z — x. 
39) x’ = 4x — y — 2, 
y = z + 2y — 2, 
2 = z— y + 2z. 
40) x° = y — 2x — 22, 
y = zx + 2y — 2, 
z' = 2y + 3z — x. 
41) z'=zxz— y +2, 


y =z+y—2, 


z = 22 — y. 
42) x = 27 + y, 

y" = 2y + 4z, 

2 T2; 


43) x' = 4x — y, 
y — 8x + y — 2, 


Z'=xz+z. 

44) x = x — y — 2, 
y =r+ y, 
2" = 3x + z. 

45) x° = 2x + 2x — y, 
y" = x + 22, 


Z = y — 2x — 2. 
46) x° = 2x — y — 2, 
y" = 3x — 2y — 32, 
Z2 =22—x+7y. 
47) x° = 3x — 2y — 2, 
y" = 3x — 4y — 32, 


z" = 2x — 4y,. 
48) zx = y — 22 — x, 
y = 4z+y, 


2 = 2x + y — 23. 
49) x° = 2x — y — 2, 
y = 2r — y — 2z, 
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2 = 22—xr+ y. 52) z° = —y — 2, 
50) x” — 2r — 3y, post 
= à 2 = —Zz — y. 
en. 53) z' = y + 2, 
51) x — 2y, y" = 3x + 2, 
y” = —)9r. z = SZ + y. 


2.4. Construire les schémas fonctionnels des systèmes d’équa- 
tions différentielles représentés sur P25. Construire le diagramme de 
correspondance entre P25 et P24. 


EXEMPLE 2.2. Mettre au point un schéma fonctionnel de l’équa- 
tion non résolue par rapport à la dérivée supérieure : 


y" (y'}2+ ay" + by = 0. (2.3) 


Moyennant une augmentation de l’ordre de l'équation (2.3), 
on obtient l'expression explicite 


y= —p{y" (y) + ay" +by]sign [(y'}?1], 
où p est un grand nombre positif. 

En définitive on obtient l'équation différentielle du troisième 
ordre résolue par rapport à la dérivée supérieure 

y"= —p{y"(y'}+ ay" + by]. (2.4) 
L'application de la méthode générale ne présente aucune difficulte. 

Remarque. D'une façon générale, les équations (2.3) et 
(2.4) ne sont pas équivalentes. Au mieux on ne peut parler que de 
l'équivalence asymptotique de leurs solutions. Se référer aux « Ele- 
ments de simulation sur calculateurs analogiques » pour plus de 
détails. 

2.5. Construire les schémas fonctionnels des équations diffe- 
rentielles représentées sur P27. Construire le diagramme de corres- 
pondance entre P27 et P26. Indiquer la fonction de l’amplificateur 
opérationnel dans les schémas de P26. 


$ 2. Génération de fonctions données 


On a à générer des fonctions données d’une variable indépendante 
lorsqu'on résout des équations différentielles à seconds membres, 
à coefficients variables et à non-linéarités spéciales. Par ailleurs la 
génération d’une fonction peut présenter de l'intérêt en soi, lors- 
qu’il est nécessaire d'étudier cette fonction avec les techniques du 
calcul analogique: tracé du graphe, recherche des points d'’extre- 
mum, d'inflexion, des Zéros, etc. 

La méthode de génération de cette fonction s'appelle méthode 
des équations différentielles déterminantes. Une équation différentielle 
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est déterminante pour une fonction donnée y ({) si y (t) en est solu- 
tion. La recherche de l'équation déterminante pour une fonction 
donnée se ramène à une dérivation successive de cette fonction et 
a l'établissement d’une relation entre les diverses dérivées. 


EXEMPLE 2.3. Composer le programme complet *) de génération 
de la fonction y = { exp (—t). 

La programmation comporte plusieurs étapes distinctes. 

1) Déduction de l'équation différentielle 
déterminante. Une dérivation successive de la fonction 
y = lexp (—t) donne 


y" = —texp (—{) + exp (—t) 
ou 
y" = —y + exp (—{); (2.5) 
y" = t'exp (—t) — 2 exp (—t) 
ou 
y" = —y" — exp (—t). (2.6) 


En ajoutant les équations (2.5) et (2.6) on obtient l'équation diffé- 
rentielle déterminante cherchée : 


y" = —2y" — y, y(0) = 0, y'(0) = 1. (2.7) 


Les valeurs numériques des conditions initialés s’obtiennent par 
report de { — 0 dans l'expression de la fonction et de ses dérivées. 

2) Schéma fonctionnel. Le schéma fonctionnel est 
composé pour les variables mathématiques. Îl n'est orienté sur 
aucune machine concrète. On dit que le schéma fonctionnel est 
banalisé. Le schéma fonctionnel exprime le principe de commutation 
des opérateurs en un système de calcul simulant le problème posé. 
Ainsi J'éequation (2.7) est câblée par le schéma fonctionnel de la 
figure 2.3, a qui est composé d'après la méthode générale. 

3) Adaptation des échelles. Les variables mathématiques sont 
assimilées en calcul analogique à des tensions électriques par l’inter- 
médiaire de coefficients dimensionnés appelés échelles. 11 faut choisir 
l'échelle des seules variables mathématiques qui sont données expli- 
citement sur le schéma fonctionnel, i.e. sont délivrées à la sortie 
d'un amplificateur opérationnel quelconque, dans notre cas, les 
variables y et y’. On ne touchera pas à la variable y” puisque celle-ci 
n'est restituée par aucun amplificateur opérationnel du schéma. 
Les variables mathématiques y” et y sont représentées par les ten- 
sions électriques u,, = my" et Uyo — Myoÿ, OÙ My et Myo Sont 
les échelles des variables y” et y. Les échelles sont assujetties à la 
condition 

—100 VLu<100 V 


_*) Un programme complet est un programme composé jusqu'au niveau du 
schéma de commutation. 
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qui traduit le fait qu'aucune variable machine ne doit quitter la 
plage +100 V en dehors de laquelle les opérateurs des calculateurs 
a tubes fournissent des résultats entachés d'erreurs. Aussi l'échelle 


est-elle choisie d'ordinaire d’après la formule m, = 100 V/| Yax |, 
car la variable machine u,, ne peut être supérieure à 100 V. En 
effet, 


_ = y 
Uyo = Myoy = 100 TEEN V <1007. 
L'adaptation des échelles implique la connaissance des valeurs 
maximales des variables. Calculons-les pour y et y. 
La variable y(t) prend sa valeur maximale dans l'intervalle 
t E[0, co], puisque y (0) = lim y (t) — 0. Cherchons donc la valeur 
{—00 


t—=t* qui annule y’. On a 
—t exp (—i) + exp (—t) = 0, 


d’où t* — 1 et Yinax —= 71. 
Pour t — O on a y’ (0) — 1 et lim y’ (t) = 0. Vérifions mainte- 
{—00 
nant la valeur de la fonction y’ (t) à l’intérieur de l'intervalle [0, œl|. 
La condition y” (t) = 0 donne l'équation t exp (—t) — 2 exp (—t) = 
— 0, d’où t** — 2 et y’ (1**) — e-°. Donc, y’ prend son maximum 1 
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pour { — 0. Les échelles sont adaptées dans des tableaux dits fa- 
bleaux d'échelle. 

Expliquons le contenu de la troisième colonne du tableau 2.1. 
Lorsqu'on adapte les échelles on arrondit toujours par excès (pour- 
quoi ?) jusqu’au plus proche nombre de la forme 10+7, 2 X 10*, 


Tableau 2.1 


Valeur maximale 
Symbole de 


Maable Echelle l'échelle 
exacte approchée 


100 

—— — 9 
0.4 250 
100 


1 1 = 100 


ei —0,368 0,4 


4 X 10+* où nest un entier, ceci dans le but de passer sans peine des 
variables mathématiques aux tensions électriques et inversement. 

4) Mise au point du schéma de commu- 
tation. Le schéma de commutation est composé pour les varia- 
bles machine. 11 est orienté vers une machine concrète, plus précisé- 
ment vers celle retenue pour résoudre le problème. Îci c’est un cal- 
culateur MH-7M dont le panneau de commutation est représenté 
sur la figure 3.3. Le schéma de commutation traduit le principe de 
connexion des opérateurs. Îl porte les numéros des opérateurs et de 
leurs entrées. Il sert à composer un problème sur le panneau de com- 
mutation de la machine (cf. fig. 2.3, b). Il se distingue légèrement 
des schémas fonctionnels par la composition et la représentation des 
opérateurs. Îl comporte notamment des potentiomètres supplé- 
mentaires destinés à adapter les échelles. Aïnsi le potentiomètre 
disposé à l'entrée (22) *) de l’intégrateur [6] possède un coefficient 
de transfert égal à 0,25; le produit de ce coefficient par le coefficient 
de transfert fixe de l’intégrateur [6] donne un coefficient résultant 
égal à 2,5. Ce nombre n'est autre que le rapport des échelles m,,'m,1=— 
— 250/100 = 2,5. I] faut adapter les échelles, car la variable —y" 
est délivrée à l'échelle m,, par l’intégrateur [5] et la variable y 
à l'échelle m,, par l’intégrateur [6]. Donc, en choisissant le coeffi- 
cient de transfert résultant de l’intégrateur [6] sous forme du quo- 
tient Myo/Myx, nous appliquons la variable —y" à l'entrée de l’inté- 
grateur [6] avec l'échelle m,, qu’elle a à la sortie de l’intégrateur [6]. 

2.6. Composer et câbler les schémas de commutation des fonc- 
tions suivantes de la variable indépendante t. Utiliser le circuit de 


*) Les parenthèses désigneront le numéro de l'entrée, les crochets, celui de 
l'opérateur. 
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commande préprogrammée pour déterminer les zéros, les extrema et 
les points d'inflexion de ces fonctions. 


Dy=1+t+é +. 22) y = (1 + cos wt) sin ot. 
2) y = a exp (—«t). 23) y = (exp (—1). 
3) y = a sin wt. 24) y = exp (—1/2) (1 + t). 
4) y = b cos wt. 25) y = cos { + sin t — 2tsin ét. 
5) y = a!. 26) y = € — 2 + 1 + exp (—(). 
6) y = t exp (—t). 27) y = cos° wt + a exp (—t). 
7) y = exp (—t). 28) y = cos* wt + 21 — 1. 
8) y = t* exp (—t). 29) y = 1 + f* + exp (—(t). 
9) y = (at + b) exp (—«t). 30) y = exp (—t) + sin t. 
10) y = a sin wt + b cos wt. 31) y = exp (1 — t) + cos t. 
11) y = a sin wt cos wt. 32) y = exp (1 — ft) + 
12) y = exp (—at) sin wt. + sin (1 — t). 
13) y — exp (—at) sin (ot +4). 33) y = * + sin £. 
14) y = exp (—at) cos wt. 34) y = t* + (* exp (—t). 
15) y = sin* wt. 35) y = 1 — exp (—t) — 
16) y = cos* wt. — t'exp (—t). 
17) y = cos t + t sin t. 36) y = sint—tocost. 
18) y = a sh wt. 37) y = sh wt + exp (—wt). 
19) y = b ch wt. 38) y = sin° w{ + sin 2wt. 
20) y = a sh wt + b ch wt. 39) y = cos* wt + sin 2wt. 
21) y = (1 + cos wt) cos wt. 40) y — 1 — 3sin? t. 


Indications et explications. 

1) On mettra en jeu des opérateurs linéaires: des sommateurs 
et des intégrateurs. 

2) Les fonctions sont simulées pour £{ € [0, ol]. Si la fonction 
ou ses dérivées ne sont pas bornées sur l'intervalle (0, œ] il faut 
se limiter à l'intervalle [0, 3] pour l'adaptation des échelles et la 
simulation. Dans ce cas il est nécessaire de prévoir l'arrêt de la 
machine à la condition t — 3 à l'aide d’un circuit de commande 
préprogrammée. 

3) Le schéma doit élaborer les variables y, y” et y” pour la déter- 
mination des zéros, des extrema et des points d’inflexion de la fonc- 
tion. Lorsque ces variables s'annulent le circuit de commande pré- 
programmée doit automatiquement stopper la machine, fixant ainsi 
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les zéros, les extrema ou les points d’inflexion de la fonction. La 
valeur numérique de la valeur indépendante correspondant à tel 
ou tel point « singulier » de la courbe peut être figée à l’aide d’un 
intégrateur élaborant la variable é. 

4) Les schémas de commutation sont génériques. Îls sont réalisés 
pour des valeurs numériques de tous les paramètres égales à l'unité. 
La durée de résolution du problème peut être modifiée pour faciliter 
l'observation sur l'indicateur électronique. 

2.7. Composer les schémas fonctionnels pour la simulation des 
équations différentielles linéaires et non homogènes et des systèmes 
d'équations différentielles linéaires et non homogènes suivants. 

Indication. Se servir de la méthode des équations diffé- 
rentielles déterminantes pour ramener les équations non homogènes 
données à des équations homogènes linéaires d'ordre plus élevé. 


1) y” — 2y" — 3y = exp (—t). 
2) y — y = 2 exp (—t) — À. 
3) y” — 3y" + 2y = sin t. 
4) y” — 5y" + 4y = At* exp (—21). 
5) y” + 3y" — 4y = exp (—4t) + t exp (—t). 
6) y” + 2y° — 3y = (* exp (—t). 
7) y” — 4y° + 8y = exp (—2t) + sin 2t. 
8) y” — Jy = exp (—3t) cos t. 
9) y” + y = exp (—t)sint. 
10) y” — 5y° = 3t° + sin t. 
11) y” + y = 4t exp (—+). 
12) y” + y" — 2t = 3t exp (—t). 
143) y” + y = 4sint. 
14) y” — 3y" + 2y = t'cost. 
15) y” + 4y = cos t + cos 3t. 
16) y” + 2y° + y = 6! exp (—{). 
17) y” + 4y° + 4y = t'exp (—2#). 
18) y” + 2y° + 2y = exp (—t) + t cos t. 
49) y” + 6y° + 10y = 3t exp (—3t) — 2 exp (—3t) cos £. 
20) y” + 8y' + 20y = 51 exp (—4t) sin 2t. 
21) y” + y = 4 exp (—t)sht. 
22) y” + 4y = sh t sin 2t. 
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23) y” + 7y° + 10y = t exp (—2t) cos 5é. 
24) y” + 2y° + 5y = 2t exp (—t) + exp (—t) sin 24. 
25) y” + 2y° + y = 2t exp (—2t) + exp (—t) sin t. 
26) y” + 8y° + 17y = exp (—At) LE — 3t sin t]. 
27) y" + 2y" + 4y' — 8y = exp (—2t) sin 2t + 26. 
28) y” + 2y = sint +t cost. 
29) y” + 6y° + 8y = 21 exp (—2t) + exp (—4t) sin t. 
30) y” + 2y° + y = t [exp (—t) — cos t]. 
31) y" +y"— y + y = 3 exp (—t) + 5tsint. 
32) y” + 6y' + 3y = exp (—3t) — 3 cos* 21. 
33) y” — Jy = exp (—3t) [1° + sin 31]. 
34) yIWN + y" = 7t — 3 cos° t. 
35) y'" + 4y° + 3y° = LE + t'exp (—21). 
36) y” — 4y° + 5y = exp (—2t) sin° t. 
37) y” + 3y° + 2y = exp (—t) cos* t. 
38) y" + 2y'" + 2y = [t + exp (—t)] sin t. 
39) yIV + 5y° + 4y = sin t cos 2t. | 
40) y” + 3y” + 2y = 2t. 
41) y" + 4y° + 3y = cht. 
42) y" + 2y + y = chtsint. 
43) x° = 3x — 4y + exp (—2t), 
y" = x — 2y — 3 exp (—2{). 
44) x = 2y — x + 1, 


y" = 3y — 2x. 

45) x° = y + 2 exp (—31), 
y = zx + À. 

46) z' = 3x + 2y + 2 exp (—t), 
y = x + 2y. 

47) z' = 4x + y — exp (—2{), 
y" = y — 2x. 


48) z° = 5x — 3y + exp (—31), 
y =z+ y +5 exp (—(). 
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49) z° = x + 2y + exp (—t), 
y = x — 5sint. 
50) x° = 2x — y + t'exp (—t), 


y = —xz + y + 3. 
91) z° = 2x + 4y — * exp (—t), 
y" = 3x + Gy. 


52) z' = 2x + 3y + 51, 

y" = 3x + 2y + 8 exp (—{). 
53) z' = x + y + 1 — exp (--t), 
y = 3x — y + 1 — sin t. 

54) z' = —y + sint, 
y —= ZI — COS f. 
55) x' = 2r — 4y, 
y = x — 3y + 3 exp (—t). 
56) x° = 2x — 3y + tsint, 
y" = x — 2y + sin° t. 
97) zx = 2x — y + t'exp (—21), 
y" = x + 2 exp (—t). 
58) x° = 2x + y + 2 exp (—{), 
y = x — 2y — cos* t. 
09) z° = 2x — y + t'exp (—t), 
y" = 2y — x — 5 exp (—t) sin t. 
60) z° = 4x — 3y + sin° t, 
y" = 2x — y — 2 sin 21. 
61) z° = x — y + 4t, 
y" = 5x — y + t exp (—2t). 
62) z° = —x + y + z + exp (—t), 
y =rz—y+2+ exp (—31), 
2 =xz+y+z+esint. 
63) x° = 3y — x — exp (—at), 
y =y+zx + exp (—Bt). 
64) z° + y" — y = exp (—{), 
2x" + y" + 2y = cos* t. 


45 
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65) x” + 2x” + x° = 1 — sin° t. 
66) x” + 3x’ + 2x = t exp (—t). 
67) x” + x = t* exp (—t). 
68) x” + x = sin“ t. 
69) zx” + x = t cos 21. 
EXEMPLE 2.4. Composer le schéma de commutation de la fonc- 
tion y — exp (—1°) sur l'intervalle [0, 41. 


1) Déduction du système d'équations dif- 
férentielles déterminantes. La fonction donnée 


Arrêt programmé 
à la cond rai 


Fig. 2.4 


est composée (est la superposition de deux fonctions). Soit z (£) = {*. 
Une dérivation de la fonction y (t) donne y” = —exp (—2) =’. 
D'autre part z° — 2t, z” 
D'où l’on déduit le système d'équations déterminantes cherché 
y” — —yz", y (0) = 1, (2.8) 
2" — 2, z° (0) = z (0) = 0. 
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2) Schéma fonctionnel. Le schéma est représenté 
sur la figure 2.4, a. Il contient un circuit de commande préprogram- 
mée pour stopper la machine à la condition { = 4. Il en résulte 
que pour reproduire la variable y (t) il n’est pas nécessaire de repro- 
duire z ({) — t*, mais seulement 2” ({). Ce cas est très fréquent lors- 
qu'on a affaire à une superposition de fonctions. 

3) Choix des échelles. Les variables —z2" et y sont 
à adapter. Leurs valeurs maximales se déterminent sans peine. 
L'adaptation est conduite dans le tableau 2.2. 


Tableau 2.2 


Valeur maximale 


Variable Echelle gi A 


exacte arrondie 


4) Schéma de commutation. Il est représenté sur 
la figure 2.4, b. Les coefficients de transfert &, B et y sont encore 
inconnus. Leurs valeurs numériques se déterminent à partir des 
conditions d'adaptation des échelles des variables. Le schéma de com- 
mutation comprend un multiplieur B-ÎT monté sur la base des ampli- 
ficateurs opérationnels [1] et [2]. Le coefficient & est donné par la 
formule suivante 


æ = 2m,,/100 = 20/100 = 0,2. 


Le produit By doit être tel qu’à la sortie de l’intégrateur [5] la va- 
riable y (t) possède l'échelle m, compte tenu de l'échelle 0,01 du 
multiplieur : 


Prenons les valeurs y = 10 et B = 1 qui correspondent à l'entrée 
(17) de l’intégrateur [5] (cf. fig. 3.3. qui représente le panneau de 
commutation de la machine MH-7M). 

2.8. Composer et câbler les schémas de commutation des fonc- 
tions suivantes : 


1) y = sin exp (—t). 6) y = cos exp (—{). 
2) y = sin cos t. 7) y = cos cos t. 
3) y = sinsin t. 8) y = cos sin t. 
4) y = sinsht 9) y = cos sh t. 


5) y =sinché. 10) y = cos ch t. 
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11) y = exp (—exp (—t)). 21) y = cos sin° t. 

12) y = exp (—sh t). 22) y = cos cos* t. 

13) y = exp (—ch t). 23) y = cos exp (—{*). 
14) y = shsht. 24) y = cos (t exp (—t)). 
15) y = sh cht. 25) y = cos ch* t. 

16) y = sin sin* t. 26) y = sin ch* t. 

17) y = sin cos* t. 27) y = cos sh* t. 

18) y = sin exp (—t°). 28) y = sin t. 

19) y = sin (t exp (—t)). 29) y = cos (*. 

20) y = sin sh*° t. 30) y = sh {*. 


Indications et explications. 

1) Lorsqu'on cherche les équations différentielles déterminantes 
d'une fonction du type sin f (t) ou cos f (ft) on a intérêt à considérer 
en même temps l’autre fonction. Si, par exemple, l’on a à simuler 
la fonction sin f ({) on peut envisager parallèlement la simulation 
de la fonction cos f (t). La simulation simultanée de fonctions sim- 
plifie l'adaptation des échelles et facilite la recherche de l'équation 
différentielle déterminante. En effet, si y, = sinf(t) et y. — 
— cos f (t), le système d'équations différentielles déterminantes 
s'écrit 

yi= val (), ys= vil (0). 


Ceci vaut également pour les fonctions sh f (t) et ch f (t). à 

2) Si l'une des fonctions ou ses dérivées ne sont pas bornées sur 
l'intervalle [0, co], pour la simuler il faut se limiter à l’intervalle 
[0, 3], ce qui implique l'interruption de la résolution à l'aide du 
circuit de commande préprogrammée. 

3) Certaines équations différentielles déterminantes sont suscep- 
tibles de posséder des solutions instables, c’est pourquoi on a inté- 
rêt à les résoudre sur machine et à évaluer l'influence de petites 
variations des conditions initiales sur le comportement de la solu- 
tion. 

2.9. Composer les schémas fonctionnels des équations différen- 
tielles suivantes. 


1) y" = —ow* sin y — ay. 

2) y” = —a* sin* y — by. 

3) y” = —0* sin y + sin t. 

4) y” = —exp (—t) sin y. 

5) [(4 + sin £) y”)]" = exp (—t). 
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6) [1 + exp (—t)] y" + 2y° + y = 0. 
7) [4 + exp (—t)) y] = a sin ot. 
8) [y*/(4 + t)]” = exp (—+). 

9) [y°/(4 + t1))}" = 1/1 + 1). 
10) y” = —ay* + sin wt. 
11) y = 1 + y° + exp (— at) sin wt. 
12) y = 1 — y* + exp (—at) cos wt. 


$ 3. Simulation de courbes planes 
données sous forme paramétrique 


2.10. Etant données les courbes suivantes sous forme paramé- 
trique, on demande les schémas de commutation 


1) délivrant le graphe des courbes F, (x, y) = 0, F, (z, y) = 0 
sur l'écran de l'indicateur électronique; 


2) déterminant les points d’intersection de la courbe avec les 
axes de coordonnées : 


3) déterminant les points en lesquels la courbe admet des tan- 
gentes parallèles aux axes de coordonnées. Construire les courbes si 
la machine comporte un traceur de courbes planes. 


1) x(t) = acost, 
y (t) = bsint 
2) z(t) = acht, 
y () = bsht 
3) z (t) = at, 
y (t) = bÉ 
4) x(t) = cost +tsint, 
y(t) =sint—toeost 
3) æ(t) = 2 cos t — cos 21, _ 
_. () = 2 sin t — sin 2 (cardioïde). 
6) z(t) —=1t—sint, 
y(t)=1—cost 
7) z(t) = cos° t 
y (t) = sin° t 
8) z (t) = cos t + 3 cost, 
y (t) = cos t sin { + 3sint 
4—0112 


(ellipse). 
0 < t < 10 (hyperbole). 
0 < t << 10 (parabole). 


< t L 4x (développante). 


0 < t < 8x (cycloïde). 
| (astroïde). 


(limaçon de Pascal). 
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9) z (t) = sin 3t + sint, | 
y (t) = cos 3t — cost; 
10) x (t) = sin 4t + sin 2t, 
y (t) = cos 4t — cos 21; 
11) x (t) = sin 5t + sin àf, 
y (t) = cos 5t — cos 3f; 
12) x (t) = sin 6t + sin 4t, 
y (t) = cos 6f — cos 4t; 
13) x (t) = sin 7t + sin 5f, 
y (t) = cos 7t — cos 5t 
44) x (t) = (cos t + sin t) exp (—t), 
y (t) = (sin t — cos t) exp (—t) 


l = 
” "7 — re : (spirale d’Archimède). 


16) z(t) = (1 +t)"! t + 1). | 
| Un == Û : =: : se 1) (spirale hyperbolique). 


17) z (t) = t — 1,25 sin t, 
y (1) = 1 — 1,95 cost, 0 < t< 67; 


SQ 
18) z(t) —t—1,5sint, 
(t) = 1 — 1,5 cos t, <t< 67; 

< 


(roses). 


(spirale logarithmique). 


y = 
19) x (t) = t—2sin t. : : “=. trochoïdes. 
y(t)=1—20cost, ° > < 87; 
20) x (t) = t — 0,75 sin t, 
y (D —1—0,75 cost, 0 SE 7n ; 
21) zx (t) =3 cos t — cos 34, 
y ) = 3 sin t — sin J!; 
22) x (t) = 4 cos t — cos 4t, 
y (t) = 4 sin t — sin 4t; (épicycloïdes). 


23) x (t) = 5 cos { — cos 5t, 
y (t) = 5 sin { — sin 51; 
24) x (t) = 6 cos t — cos Gt, 
y (t) = 6 sin { — sin 6t 
25) x (t) = 6 cos t — 1,25 cos6t, 
y (t) = 6 sin { — 1,25 sin 6; 
26) x (t) = 6 cos t — 1,5 cos Gt, (épitrochoïdes). 
y (t) = 6 sin t — 1,5 sin 6t; 
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27) x (t) = 6 cos t — 2 cos 64, 
y (t) = 6 sin { — 2 sin 6f; 
28) x (t) = 8 cos t — 2 cos 8t, 
y (t) = 8 sin t — 2 sin 8t 


29) x (t) = 2 cos t + cos 21, 
y (t) = 2 sin t — sin 2t; 
30) x (t) = 3 cos t + cos 3t, 
y (t) = 3 sin { — sin 3t; 
31) x (t) = 4 cos t + cos 4t, (hypocycloïdes). 
y (t) = 4 sin t — sin 4t; 
32) x (t) = 5 cos t + cos 5, 
y (t) = 5 sin t — sin 5t; 
33) z (t) = 6 cos t + cos 64, 
y (t) = 6 sin { — sin 6t 
34) zx (t) = 2 cos t + 1,25 cos 21, 
y (t) = 2 sin t — 1,25 sin 2t; 
35) x (t) = 2 cos t + 1,5 cos 24, 
y (t) = 2 sin t — 1,5 sin 2f; 
306) x (t) = 3 cos t + 1,25 cos 34, 
y (t) = 3 sin { — 1,25 sin 3t; 
(hypotrochoïdes). 


37) z(t) = 3 cos t + 1,5 cos 34, 
y (t) = 3 sin t — 1,5 sin 3t; 
38) zx (t) = 4 cos t + 1,5 cos 4t, 
y (t) = 4 sin t — 1,5 sin 4t; 
39) x (t) = 5 cos t + 1,5 cos 5t, 
y (t) = 5 sin t — 1,5 sin 5t 


co 


Indication. On recommande de programmer la machine 
pour le cas générique comme dans l’exemple 2.5 suivant. Une fois 
le schéma générique obtenu, sa réalisation pourfdes valeurs numé- 


riques données ne présente aucune difficulté. 


EXEMPLE 2.5. Simuler l'hypotrochoïde 
x (t) = © cos { + À cos wf, 
y (t) = o sin t — A sin wf, À > 1. 


4 


(2.9) 
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Au promier examen il semble plus commode de simuler x ({) et y (t) 
par sommation des fonctions cos {, sin {, cos w{ et sin wf. Mais 


dans ce cas la machine ne délivrera pas les variables x (£) et y (t) 


indispensables au tracé de l'hodographe *) F, (x, y) = 0. Nous 
allons procéder autrement. Nous allons chercher les équations diffé- 
rentielles déterminantes des fonctions x (t) et y (t). 

1) Recherche des équations différentiel- 
les déterminantes. Une dérivation de (2.9) donne 


z = —6 sin { — Ào sin wt, 
| (2.10) 
y = & cos  — À cos wtf. 


En portant dans (2.10) les expressions de w cost et wsint, on 
obtient 


z = —y — À (1 + ©) sin ot, 


(2.11) 
y = z— À (1 + ©) cos wt. 
Une dérivation de (2.11) donne 
ZT = —} — Ào (1 + o) cos ot, 
É ) (2.12) 


y = 2 + Ào (1 + o) sin of. 


En portant dans (2.12) les expressions Àw (1 + «) cos wt et 
ko (1 + o) sin ot tirées de (2.11), on obtient en définitive le sys- 
tème d'équations différentielles déterminantes 


D (2.13) 
y = —(o — 1) x — oy 


avec les conditions initiales 
z(0O) =w+X, zx(0) = 0, 


y (0) = 0, y (0) = w (1 — À). 


2) Schéma fonctionnel. Il est composé pour (2.13) 
d'après la méthode générale (fig. 2.5). Les variables y (t) et zx (t) 
sont appliquées respectivement sur les plaques de déviation verticale 
et horizontale de l'indicateur électronique sur l'écran duquel est 
représentée la courbe F, (x, y) — O0 qui n’est autre que l’hypotro- 


choïde allongée donnée. Si l’on applique la variable y (t) sur les 


*) En Afécanique on définit l'hodographe de la vitesse comme suit. A partir 
d’un. point arbitraire du plan (appelé pôle) on trace des vecteurs équipollents 
à la vitesse du point en mouvement. Les extrémités de ces vecteurs décrivent une 
courbe appelée hodographe de la vitesse. Ici, le pôleest l'origine des coordonnées. 
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plaques de déviation verticale et De (t) sur les plaques de déviation 


horizontale, on obtiendra l’hodographe de F, (x, y) = 0. Pour 
stopper automatiquement la machine aux points d'intersection de 
l'hypotrochoïde avec l'axe des z, on applique la variable zx (t) à 
l'entrée du circuit de commande préprogrammée. On peut stopper 
automatiquement la machine en tous les points où x (ft) = 0, en 


Z;Y; E;ÿ CCP 


Fig. 2.5 


modifiant la condition d'arrêt automatique (arrêt lorsque zx (t) 
passe par zéro en décroissant ou en croissant). De façon analogue, 
pour déterminer les points d'intersection de l'hypotrochoïde avec 
l'axe des y, on applique la variable y (t) à l'entrée du circuit de 
commande préprogrammée; pour déterminer les points en lesquels 
l'hypotrochoïde admet une tangente horizontale ou verticale, on 


applique respectivement les variables —y (t) et —zx (t) à l'entrée 
du circuit de commande préprogrammée. 

3) Choix des échelles des variables. On 
adapte l'échelle des variables qui sont « explicitement » représentées 


Tableau 2.3 


Valeur maximale Symbole de 
Variable de la variable Echelle l'échelle 


& (1 + À) 


mi 


« (1 + À) 
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sur le schéma fonctionnel. Dans le schéma de la figure 2.5, on adap- 


tera l'échelle de zx (£), y (t), z (£) et y (t). Cette procédure est con- 
duite dans le tableau 2.3, où les échelles des variables z (t), y (t) 


d'une part et z (4), y (t) d'autre part sont égales. Les valeurs maxi- 
males des variables x (t) et y (t) sont tirées de (2.9), celles de z (4) 


et y (t) de (2.10). Le choix des échelles correspond au cas générique, 
i.e. il est réalisé pour une hypotrochoïde quelconque. 

4) Schéma de commutation. Il est représenté sur 
la figure 2.6. Contrairement au schéma fonctionnel, il comporte des 


Gk:{ek: 5} V7] 


Fig. 2.6 


potentiomètres supplémentaires pour l'adaptation des échelles. 
Les. coefficients de transfert de ces potentiomètres dépendent des 
échelles des variables et des coefficients de transfert du schéma fonc- 
tionnel : 

@({+A) 7 = o+A 


nine == 9 WA 27 ES 


G3 = 9% 0 1. 


Le schéma de commutation est générique, i.e. il simule une hypo- 
trochoïde quelconque. Pour concrétiser le problème il suffit de pren- 
dre des coefficients de transfert &,, &«. et &«; appropriés. 

Trouvons les valeurs numériques des coefficients de transfert 
Qi» +, 3 Correspondant à w = 8 et À = 2, le temps de simulation 
étant ralenti de 10 fois pour permettre le tracé de la courbe avec un 
traceur. On obtient 


8(1+2) _ = 
410.10 — 0,24, a = 5 = 0,333, = —— — = 0,1. 


$ 4] SIMULATION DE COURBES DONNÉES IMPLICITEMENT 55 


Fig. 2.7 


Le schéma de commutation correspondant à ces coefficients est 
représenté sur la figure 2.7. 


$S 4. Simulation de courbes planes données implicitement 


Pour simuler une courbe plane d’équation F (x, y) = 0, il faut 
passer de la forme implicite à la forme paramétrique zx (t), y (t) où 
t est un paramètre. Il existe un procédé qui permet de trouver immé- 
diatement les équations différentielles déterminantes de xt) et 
y (t) sans en calculer les expressions analytiques. Etudions ce pro- 
cédé. 

Soit une courbe plane F (x, y) = 0. Si le point de coordonnées 
x (t) et y (t) appartient à cette courbe, la dérivée totale de la fonc- 
tion F (x, y) par rapport à t s’annule en ce point. Donc, quelle que 
soit la représentation parameétrique de la courbe, on a 

0F dr , 0F dy 

x Toy à — 
Cette égalité a visiblement lieu si la représentation paramétrique 
est telle que 


dr oF dy aF 
AE — METE HT: (2.14) 
où o est un nombre arbitraire différent de zéro. 

Les équations différentielles (2.14) sont équations différentielles 
déterminantes pour les fonctions x (t) et y (t). Il faut se donner les 
conditions initiales x (0) et y (0). L'une des conditions initiales est 
choisie arbitrairement, par exemple x (0) — x,. L'autre condition 
est la solution de l'équation F (x,, yo) = 0 par rapport à y,. Les 
conditions initiales x, et y, signifient que la courbe est simulée à 
partir du point (x,, ÿ.). Le point initial (x,, y,) partage la courbe 
en deux branches, l'une simulée pour w > 0, l'autre pour w << 0. 
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Ainsi, en principe, le problème posé se ramène à la résolution 
sur la machine du système d'équations différentielles 


T= —0——, z (0) = %o, 

[= 
F. (2.15) 
= 0 Ts y (0) = Yo, 
où F(xzo, yo) = 0. 

Mais la simulation de la courbe plane par câblage de l’équation 
(2.15) est susceptible de conduire à un résultat erroné en raison de 
la grande influence des petites erreurs des conditions initiales sur 
la solution. Cette influence est particulièrement sensible lorsque le 
système (2.15) possède une solution instable. Aussi est-on tout natu- 
rellement poussé à « corriger » le système (2.15) de telle sorte que la 
courbe donnée soit simulée assez bien quelles que soient les con- 
ditions initiales. Ceci est possible si l’on remplace le système (2.15) 
par le système: 


TZ = —0 PAF, 
(2.16) 
0F ôoF 
JR APE 


où À © 1, et les termes FF et FS décrivent le mouvement du 
point représentatif dans la direction du minimum de la fonctionnelle 
5 F (x, y), un minimum nul réalisé sur la courbe F (x, y) = 0 *). 
Donc, tant que le point se trouve sur la courbe F (x, y) = 0 les 
termes AF- et AFS sont nuls et (2.16) s’identifie à (2.15). 


EXEMPLE. 2.6. Simuler la courbe plane 
F(z, y)= a —1—y =0. 


En vertu de (2.16) on a le système d'équations différentielles 
— oÿ + ÀtF, z (0) = ze, 


= —@x + ÀyF, y (0) — yo. 


Sautons les étapes intermédiaires de la programmation et indiquons 
directement le schéma de commutation simulant le mouvement du 
point sur un cercle. Le schéma de commutation de la figure 2.8 
correspond au cas —1<r<<1, —1<y<1, 0O<a<i; les 
échelles sont adaptées: a = 0,5, w — 0,5. Les quantités À, w et a 
peuvent être modifiées à l’aide des potentiomètres portant des in- 


(2.17) 


œe Re 


*) Le procédé exposé se base sur la méthode des fonctions de pénalisation, 
cf. Eléments de simulation sur calculateurs analogiques, chapitre 10, & 4. 
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Fig. 2.9 
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scriptions appropriées. Les opérateurs [5], [6] et [13] possèdent des 
sorties plus puissantes. Ce fait est signalé sur le schéma par un noir- 
cissement du sommet désignant la sortie de l'opérateur. Toutes les 
échelles ont été omises pour alléger le déchiffrage de la figure 2.8. 

Sur la figure 2.9 sont représentés les résultats de la simulation 
du système (2.17) à l’aide du schéma de la figure 2.8 pour diverses 
valeurs des conditions initiales x, et y, et une valeur relativement 
petite de À. On constate que la trajectoire x ({), y (t) dégénère en un 
cercle lorsque { —> co. Ceci étant, si le point initial (x,, y.) de la 
trajectoire se trouve à l’intérieur du cercle, la trajectoire tendra 
asymptotiquement vers le cerle en l’'enroulant de l’intérieur. Si le 
point initial se trouve à l’extérieur du cercle, la trajectoire l'enrou- 
lera de l'extérieur. Plus À est grand, plus la trajectoire dégénère 
rapidement en cercle. On peut modifier le rayon du cercle à l’aide 
du potentiomètre en agissant sur l'entrée (40) du sommateur [13]. 
À noter que les trajectoires de cette nature s'appellent cycles limites. 

La notion de cycle limite a été introduite par H. Poincaré. Les 
cycles limites sont étudiés dans la théorie qualitative, ou géométri- 
que, des équations différentielles. Les cycles limites et les applica- 
tions techniques qui y sont rattachées sont très importantes en théo- 
rie de la commande automatique, en radiotechnique et dans nombre 
d’autres disciplines. Dans notre exemple, nous avons délibérément 
choisi une trajectoire qui soit cercle limite, i.e. toute trajectoire 
issue d’un point quelconque du plan z0Oy dégénérera en ce cerle. 
Théoriquement, le point x, = y, = 0 est un point d'équilibre pour 
le système (2.17). Mais ce po nt est instable, donc en raison des 
erreurs inévitables des opérateurs, la machine ne peut délivrer la 
solution æ(t) = y(t) = O vérifiant les conditions initiales x, — 


———— 
—— — 


2.11. Simuler les courbes données sous forme implicite 
1)(c—a) +(y—b,) = R° 6) rt + y -- 2ary = 0. 


2) EME + UÈE 2 4 7) (et + y — ax) -= b (sy). 
) a? 2 b2 Sn : . se A 
5 2 8) (x° + y*)* — 2c° (z* — y*) = 
D ve re À = at — 6, 
4) ( — a) = k (y — bi). 9) (x° + y°)* — 2c* (x° + y°) =0. 
= 10) y* = x? az : 
5) a + y = 3azy. ——. 


Indications. 

1) Choisir les paramètres a, a;, b, b,, R, etc. tels que la courbe 
soit simulée dans le carré —1 7 <1, —1 << y <1. 

2) Composer le schéma de commutation à l’aide de l'équation 


(2.16). En cas de besoin utiliser deux machines ou plus du type 
MH-7M. 
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Fig. 2.10 


3) Simuler les courbes avec les équations (2.15) et (2.16). Rele- 
ver et expliquer les causes des divergences; tracer les courbes avec 
un traceur. 


4) Déterminer les points d'équilibre instable de chaque trajec- 
toire d'équation (2.16). 

2.12. Réaliser les polygones réguliers suivants: 
1) le carré |x|+|y|= 1; 
2) le carré [z—y|+|z+y]|= 2; 


3) l'héxagone | 2y — 1 | + y+1+- sl 4; 


4) l'octogone [z|+|y| + (ls 4 + z+ y) =V2#1. 
ExEMPLE 2.7. Simuler la courbe plane 
Fa y=1z|+lyl—-a=0. 
En vertu de (2.16) on obtient le système d'équations différentielles 
z = —@ sign y — ÀF signz, zx(0) = x, (2.18) 
y = © sign z — ÀF sign y, y (0) = ÿo- 
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Comme dans l'exemple 2.6, omettons les étapes intermédiaires de 
la programmation et donnons immédiatement le schéma de commu- 
tation (fig. 2.10). Les opérateurs [14] et [13] sont des relais parfaits 
(cf. « Eléments de simulation sur calculateurs analogiques », chap. 2, 
$ 10). Les opérateurs [14] et [13] délivrent la fonction sign dont la 
valeur correspond à la tension électrique +100 V à l’échelle adoptée. 
Les sorties des opérateurs [9], [13] et [14] sont plus puissantes. 

La figure 2.11 représente la solution machine du système (2.18) 
d'après le schéma de la figure 2.10. On remarque l'existence d’un 
cycle limite. La] figure 2.12 représente le graphe, en gros plan, de 
la solution du système (2.18) privé des termes « pénalisants » : 


Z = —0 Sign y, y = = @ Sign z 


au voisinage de z5 — ÿo — 0. La trajectoire est une spirale dérou- 
lante formée de segments de ligne polygonale. 


$ 5. Représentation graphique de courbes gauches 


ExEMPLE 2.8. Composer le schéma de commutation reproduisant 
la courbe gauche 
Z, = exp (—Bt) cos wf, 2, = at, z4 = exp (—Bt) sin wt. (2.19) 
La programmation s’amorce par la recherche des équations diffé- 
rentielles déterminantes qui se trouvent sans peine pour les expres- 
sions considérées : 


a= is z, (0) =1, 


Z3 = re z, (0) = 0. 

Le schéma de commutation est celui de la figure 2.13. Les variables 
Vi =Z — 0,523, Ya = Z9 — 0,52; qui assurent la représentation 
de la courbe (cf. Annexe) sont élaborées par les sommateurs [3] 
et [4]. La figure 2.14 représente la courbe et sa projection délivrées 
par le schéma de Îa figure 2.13. Pour obtenir le tracé des projections 
on a commuté les entrées des sommateurs [3] et [4] en s'inspirant 
du tableau de commutation de la figure A.3 (page 177). La figure 2.15 
représente diverses courbes réalisées avec le schéma de commutation 
de la figure 2.13. En haut sont représentées quatre spirales ayant 
chacune son origine sur un axe de coordonnées. Au milieu — leurs 
projections axonométriques sur le plan z,0z, (pour un observateur 
situé au sommet du cône qu'elles engendrent) ; en bas — les projec- 
tions sur le plan z,0z, (pour un observateur se trouvant à la base 
du cône). Si l’on rompt les liaisons entre les intégrateurs [5] et [6] 
aux entrées (17) et (21), les équations initiales de la courbe (2.19) 
deviennent 

Z, = exp (—Bt), z, = at, 23 — exp (—Bi). (2.20) 


Fig. 2.13 


2 


Fig. 2.14 
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Fig. 2.16 


Les courbes correspondantes sont les génératrices du cône engendré 
par la révolution d'une exponentielle autour de son asymptote Oz, 
(cf. haut de la figure). Si l’on ouvre les liaisons des intégrateurs [5] 
et [6] aux entrées (20) et (24), on obtient la ligne hélicoïdale 


Z1 = Sin Of, Zo = Al, Z3 = COS Gé. (2.21) 
La figure 2.16 reproduit les courbes (2.19) et (2.21). Le point repré- 
sentatif décrit une trajectoire qui commence en À et se déplace dans 
le sens de la flèche. Au début c’est une ligne hélicoïdale située sur 
un cylindre ;:elle se transforme ensuite en une spirale portée par 
un cône dont la génératrice a la forme d’une exponentielle; la spirale 
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devient une ligne hélicoïdale qui, à son tour, se transforme en spirale 
en même temps que « change de signe. Toutes ces transformations 
ont été effectuées, machine en marche, par rupture et rétablissement 
des liaisons mentionnées entre Jes intégrateurs [5] et [6]. 

Voir le chapitre 4 pour la simulation d'autres courbes gauches. 

2.13. Tracer sur l'indicateur électronique les courbes gauches 
suivantes. 

Indications. Les courbes 1) à 10) sont simulées à l’aide 
du schéma du chapitre 4, $ 10. Les courbes 11) à 15) à l’aide du sché- 
ma du chapitre 4, $ 8. 


1) z, = sin 0,36, Z3 = (1 + 0,4 cos 106) X 
za — Cos 0,34 cos 1014, X sin 0,34. 
z4 —= cos 0,34 sin 104. 9) z, = sin 0,34, 
2) z, = sin 10, Ze = (1 + 0,3 cos 0,34) X 
z, — cos 10f cos 0,3t, X cos 104, 
z, = cos 10t sin 0,34. za = (1 + 0,3 cos 0,3) X 
3) z,,— sin 0,34, x sin 104. 
ze = (1 + cos 0,3t) cos 104, 10) z;, — sin 10, 
Za = (1 + cos 0,34) sin 104. Ze —= (1 + 0,3 cos 101) X 
4) z, = sin 104, X cos 0,3, 
ze —= (1 + cos 104) cos 0,36, Za = (1 + 0,3 cos 10£) X 
z23°= (1 + cos 10) sin 0,34. x sin 0,34. 
5) ‘z, = 0,7 sin 0,3, 11) z, = exp (—t) cos 104, 
z, = cos 0,3{ cos 10, Z, — exp (—{), | 
z, —=_cos 0,3t sin 104. Zs = exp (—i) sin 104. 
6) z, — 0,7 sin 406 12) z, = exp (—t) cos 104 
? , me 94 | 
z, = cos A0t cos 0,34, ns Cp 


z3 —= exp (—t) sin 104. 


43) z, = exp (—t) cos 104, 
za — exp (—0,3t) sin t, 


z, = cos 104 sin 0,3t. 
7) z, = 0,4 sin 0,34, 


24 = (1 + 0,4 cos 0,5) X Z. — exp (—t) sin 10f. 
x cos 10, 44) z, — exp (—t) cos 106, 
za = (1 + 0,4 cos 0,3t) X Zo = Sin t, 
x sin 404. zs —= exp (—t) sin 10. 
8) _z, — 0,4 sin 10t, 45) z, = exp (—t) cos 10f, 
ze = (1 + 0,4 cos 104) X 2, = 1/4 +), 


X cos 0,3f, zs = exp (—t) sin 104, 


CHAPITRE 3 


TRAVAUX PRATIQUES 


$ 1. Généralitéssur le calculateur analogique MH-7M 


Le calculateur analogique MH-7M est une machine modulaire de 
table, composée de 18 amplificateurs opérationnels dont 16 servent 
à former des opérateurs: sommateurs, intégrateurs, dériveurs, géné- 
rateurs de fonctions, etc. Il permet de résoudre des systèmes d'équa- 
tions différentielles non linéaires (jusqu’au sixième ordre) de la 
forme 


25 = y (Lis Los eZ) T(0)=Zzyy, 1=1,2,...,06, 


où x; (t) sont des fonctions inconnues de la variable indépendante t. 

Si le problème posé nécessite plus d'opérateurs que n’en dispose 
la machine, on peut adjoindre à cette dernière d’autres machines 
pour accroître ses performances. 

Les éléments de base de la machine sont le bloc principal (résol- 
vant), le bloc d'alimentation et l'indicateur électronique. La figu- 
re 3.1 représente le bloc principal et l’indicateur électronique; la 
figure 3.2, le panneau avant du bloc principal sur lequel sont dispo- 
sés les éléments destinés à la commande et au contrôle de la machine ; 
la figure 3.3, le panneau de commutation de la machine MH-7M. La 
description détaillée de la machine MH-7 est donnée dans les « Ele- 
ments de simulation ». Les multiplieurs, diviseurs et générateurs 
de fonctions du calculateur MH-7M sont plus perfectionnés que ceux 
du MH-7. Ces opérateurs ne se connectant pas comme dans le cal- 
culateur MH-7, des modifications appropriées ont été apportées 
dans le panneau de commutation. Par ailleurs, le calculateur MH-7M 
utilise le bloc d'alimentation B3C-1. 


Nature et succession des opérations 


1. Mise sous tension. Avant de brancher la machine il est indis- 
pensable de lire la notice d'emploi, les instructions d'exploitation 
et les mesures de sécurité. 

Mettre la source d'alimentation B3C-1 et l'indicateur électro- 
nique sous une tension de 220 volts. 
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Fig. 3.3.— Panneau de commutation 
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Mettre en marche la source d'alimentation. Mettre pour cela en 
position CETB (SECTEUR) le tumbler inférieur gauche, situé sur 
le panneau avant de la source d’alimentation B3C-1: les deux lam- 
pes de signalisation des coins supérieurs gauche et droit s’allumeront 
aussitôt. Au bout d’une ou de deux minutes, mettre le tumbler 
droit dans la position AHOJI (ANODE). S'assurer à l’aide du volt- 
mètre à aiguille incorporé, en modifiant la position du commutateur 
d’échelles, de l'existence des tensions électriques suivantes: 
+100 volts, +350 volts. —100 volts, —350 volts. Si les valeurs 
indiquées par le voltmètre diffèrent des valeurs nominales il est 
possible de les régler à l’aide des potentiomètres dont les fentes 
aboutissent au panneau de la source d'alimentation. 

Toutes les manipulations ultérieures s'effectuent à l’aide des 
tumblers et des commutateurs disposés sur le panneau de commande 
(fig. 3.2) du bloc principal et sur le panneau de commutation 
(fig. 3.3). Fermer les tumblers —220 B (220 volts) et —26 B (26 volts) 
situés en bas et à gauche du panneau du bloc principal: la lampe de 
signalisation disposée au-dessus du tumbler 220 B s'allumera. 
Ceci achève la mise sous tension de la machine. Toutes les opéra- 
tions nécessitées par la préparation de la machine à la résolution 
d'un problème quelconque s'effectuent machine branchée. 

Pour mettre la machine hors tension on refait le chemin inverse. 
On ouvre d’abord les tumblers — 26 B et — 220 B, ensuite le tum- 
bler AHOÏ (ANODE) et enfin le tumbler CETB (SECTEUR). 

2. Réglage des zéros des amplificateurs opérationnels. Les ampli- 
ficateurs opérationnels de la machine MH-7M possèdent une dérive 
du zéro de la tension de sortie. La dérive qui est une source d'erreurs 
est déterminée par une mesure de la tension de sortie, bornes d'entrée 
court-circuitées. A défaut de l’éliminer on peut en atténuer les 
effets pendant un certain temps. Pour ce faire on vérifie périodique- 
ment la valeur de la tension de sortie de l’AO, puis on règle au besoin 
l'amplificateur opérationnel. Cette opération constitue ce qu’on 
appelle le réglage des zéros de l’AO. 

Le réglage des zéros se fait dans les positions YCTAHOBKA 
HYJIA (RÉGLAGE ZERO) et TNOJATOTOBKA (PRÉPARATION) 
respectivement des tumblers T4 et T5 (cf. fig. 3.2). Il est réalisé 
avec les 18 boutons P12 des résistances variables. Les boutons sont 
disposés dans le coin bas gauche du panneau de commande du bloc 
principal. Chaque bouton porte le numéro de l’amplificateur opé- 
rationnel qu’il commande. Pour mesurer la tension de sortie des 
amplificateurs on utilise le voltmètre V1 (avec le zéro central). 
Le tumbler T13 doit être en position H3MEPEHME (MESURE). 
La sortie de l’amplificateur contrôlé doit être reliée à la douille +V1 
(située en bas et à droite du panneau de commutation) du voltmètre 
VI par un cordon de connexion. (Les sorties des amplificateurs 
aboutissent au milieu de la partie inférieure du panneau de commuta- 
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tion.) En tournant le bouton approprié P12 on effectue le réglage du 
zero d'abord sur l'échelle 2,5 V (1114 est en position basse), ensuite 
sur l'échelle 0,1 V (I114 en position haute). Le réglage est terminé 
lorsque l'aiguille du voltmètre V1 indique le zéro de l'échelle 0,1 V. 

3. Etude de la précision de réglage des tensions constantes à l'aide 
du diviseur de tension. Le réglage et la mesure des valeurs constantes 
arbitraires des tensions électriques constituent une importante phase 
dans la manipulation de la machine. En calcul analogique, les ten- 
sions électriques jouent le rôle de variables mathématiques. Aux 
tensions constantes correspondent des variables mathématiques 
constantes. 

1) Cäbler sur le panneau de commutation le schéma représenté 
sur la figure 3.4. Le rond contenant la lettre JL désigne un potentio- 
mètre à trois décades dont la borne d'entrée (la douille BX) est 
raccordée, avec un cordon de connexion, à une douille quelconque 
du panneau de commutation, portant l'inscription + 100 volts, et la 


Commutateur ”Diviseur 
de tension 


Fig. 3.4 


borne de sortie (douille BbIX) est reliée par le répétiteur des sorties 
à la douille + V1 du voltmètre V1 et à la douille + V2 du volt- 
mètre V2 (zéro à gauche). Les douilles BX (ENT), BbIX (SOR), 
Vi et V2 sont toutes situées à côté du rond symbolisant le diviseur 
de tension. Le commutateur AEJIHTEJE HATPAKEHHA (DIVI- 
SEUR DE TENSION) est placé sur le panneau avant du bloc prin- 
cipal, sous le voltmètre V1. Le commutateur en question est consti- 
tué de trois commutateurs à lamelles: celui de gauche possède une 
graduation égale au dixième (la tension appliquée étant de 100 volts, 
la valeur d'une graduation est de 10 volts); celui du milieu, une 
graduation égale au centième, celui de droite, une graduation égale 
au millième. Le coefficient de transfert du diviseur de tension est 
égal à la somme des valeurs indiquées par les trois commutateurs 
à lamelles. Cette somme est inférieure à l’unité. 

2) Afficher les valeurs numériques des tensions u,; indiquées 
dans la deuxième colonne du tableau 3.1 à l’aide du DIVISEUR 
DE TENSION; inscrire au tableau les indications des voltmètres 
V1 et V2, ï.e. les valeurs v,;, v.,. Calculer les erreurs quadratiques 


…] 
tŸ 
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Tableau 3.1 


moyennes 


1U : 1u 


; A : 
Ji 70 D (Un, — 015}; PT V 70 > (UD, — V5). 
: — 


J=! 


3) Enumérer et expliquer les causes éventuelles des écarts entre 
les indications du diviseur de tension et celles des voltmètres V1 
et V2. 

4. Etude de la précision d'affichage des tensions constantes avec la 
source ATAJIOHHOE HATPAKEHHE (TENSION ÉTALONNÉE). 
La source TENSION ÉTALONNÉE de la machine MH-7M est 
construite sur la base de l’amplificateur opérationnel [17] suivant 
le schéma de la figure 3.5, a, où la fonction de réaction est assumée 
par un potentiomètre à trois décades. La tension d'entrée est ordi- 
nairement égale à +100 volts ou —100 volts selon le signe de la 
tension étalonnée. La tension étalonnée de sortie est réglée par modi- 
fication du coefficient de transfert du potentiomètre à l'aide du 
commutateur TENSION ÉTALONNÉE situé sur le panneau avant 
du bloc principal sous le voltmètre V2. 

1) Monter le schéma de la figure 3.5, b. L'entrée BX (ENT) 
de TENSION ÉTALONNÉE est disposée à droite des douilles 
+V1 et V2 à côté du tumbler STAJJOHHOE HAIÏPAREHTME, 
TIPOTPAMMHEIH PEXKHM (TENSION ÉTALONNÉE — MODE 
PROGRAMMÉ) qui doit être mis en position TENSION ÉTALON- 
NÉE. L'entrée BX est raccordée à la plus proche douille +100 volts, 
entre les douilles + V1 et V2. La sortie BbIX (SOR) de TENSION 
ÉTALONNÉE est reliée par le répétiteur de sorties aux entrées 
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des voltmètres V1 et V2. Le tumbler T4 doit être mis en position 
PABOTA (CALCUL). La position du tumbler T5 importe peu. 
2) Régler à l’aide du commutateur TENSION ÉTALONNÉE 


Tableau 3.2 


Sonme 


(uéty —t'15)° 


(ut; —v23)° 


les valeurs des tensions us: figurant dans la deuxième ligne du ta- 
bleau 3.2, inscrire au tableau les valeurs v,; et v,; indiquées par les 


- Tension 


NS d L,7 Sor L. 
Mode u 
programmé Le 


Commutateur "Tension 
eétalonnée” 


Fig. 3.5 


voltmètres V1 et V2. 
Calculer les erreurs quadratiques moyennes 


g 10 /[ 1U 
| 1 
O1 — V = 70 D (ets — 015), O2 = V 70 2 (Uetÿ — V2). 
3=1 j=1 


3) Enumérer et expliquer les causes éventuelles des écarts entre 
les indications du commutateur TENSION ÉTALONNÉE et les 
voltmètres V1 et V2. 
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0. Détermination de la valeur et de la plage de variation des coef- 
ficients de transfert du sommateur (1] et de l’'intégrateur [5]. On déter- 
mine le coefficient de transfert d’un opérateur à une entrée quelcon- 
que en appliquant à ladite entrée une tension constante connue 
(par exemple, la tension wst de la source TENSION ÉTALONNÉE) 
et en mesurant la tension de sortie u,.. de l'opérateur. Le quotient 


est numériquement égal au coefficient de transfert cherché. 


Usor 

U 

EE se des coefficients de transfert des opérateurs s’effectue 
dans la position CALCUL du tumbler T4 et PRÉPARATION du 
tumbler T5. Dans cette position des tumblers T4 et T5, les conden- 
sateurs (1 uF) sont automatiquement déconnectés dans les circuits 
de réaction de tous les intégrateurs et à leur place sont branchées 
des résistances équivalentes (sur le plan du coefficient de transfert) 
de 1 MQ. 

1) Déterminer le coefficient de transfert à l'entrée (1) du somma- 
teur [1]. À cet effet on organise d’abord l'entrée (1) du sommateur 
[1], i.e. l'entrée de l’AO! est reliée par un cavalier ou un cordon 
de commutation à la douille de la résistance d'entrée de 0,1 MQ 
portant le numéro 1. Ensuite à l'entrée (1) qui est située au centre 
du champ de composition (les 64 entrées sont comprises dans le rec- 
tangle portant l'inscription BXOJIbI) on applique la tension éta- 
lonnée us — 10 volts. Avec le voltmètre V2 on mesure la tension 


Tableau 8.3 
Sommateur [1] Intégratcur sommateur [5] 
Contre-réaction Contre-réaction 
ne MR | o.1mo nie iuF | 1Mmo 
(4) | (20) | 


à la sortie du sommateur [1] pour diverses valeurs des résistances 
de contre-réaction (1 MQ, 0,1 MO). Calculer les valeurs des coeffi- 
cients de transfert correspondant aux résultats des mesures et les 
porter dans le tableau 3.3. 

2) Répéter 1) pour les entrées (2), (3), (4), (17), (18), (19), (20) 
indiquées dans le tableau 3.3. 
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Modifier si nécessaire la tension us appliquée aux entrées des 
opérateurs. Indiquer les valeurs maximale et minimale des coeffi- 
cients de transfert aux entrées paires. 

6. Mise au point des coefficients de transfert des opérateurs. On 
demande de monter le schéma de commutation de la figure 3.6, a 
et d'afficher les coefficients de transfert indiqués dans le tableau 3.4. 


Tableau 3.4 


Variante 


0 1 2 3 4 5 6 7 8 9 
Coef- 
ficient 
Lo.1e|0.18 Lo | 0,51 Lo. Lo. or 0.68 0.28 0,24 
| o.6|0,ss | 0,2 Lo.12 0.13 o- -33 oi [0.08 0.39 [0,3 
0.72 o.s|0. lo 612 36 0,04 | 0.607 lo 4 0,897 lo. 


La composition du schéma commence par l'organisation des 
opérateurs. S'agissant des sommateurs et intégrateurs, le terme 
« organisation » signifie : le branchement des entrées indiquées dans 


Fig. 3.6 


le schéma à l'entrée de l’AO (organisation des entrées), la connexion 
des éléments constituant les circuits de réaction (organisation des 
réactions). Toutes les connexions électriques sont effectuées sur le 
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panneau de commutation à l’aide de cavaliers ou de cordons de com- 
mutation. 

Le câblage ultérieur du schéma de commutation comporte le 
branchement des entrées de certains opérateurs aux sorties d’autres 
(organisation des connexions des opérateurs). Cette opération est 
effectuée sur le champ de composition qui occupe la partie gauche 
du panneau de commutation. Le champ de composition est l’abou- 
tissement des entrées et sorties des opérateurs. Il est compris dans le 
rectangle portant l'inscription BHIXOJHEI (SORTIES). A l'inté- 
rieur se trouve un autre rectangle portant l'inscription BXOJIEI 
(ENTRÉES), qui contient toutes les 64 douilles d'entrée des opé- 
rateurs. Entre le rectangle intérieur et le rectangle extérieur se 
trouvent les sorties des 16 opérateurs. Celles-ci sont reproduites en 
quatre exemplaires pour faciliter la commutation. Les entrées et 
sorties sont reliées par des cordons de commutation. 

Les coefficients de transfert des sommateurs et intégrateurs sont 
mis au point dans la position CALCUL du tumbler T4 et PRÉ- 
PARATION du tumbler T5. Les coefficients de transfert des som- 
mateurs 0,1 ; 0,2 ; 0,5; 1; 2; 5; 10; 50 et les coefficients de transfert 
des intégrateurs 1, 2, 5, 10, 50 sont établis par une organisation 
appropriée des entrées et des circuits de réaction des opérateurs. Si 
l’on a besoin d’autres coefficients de transfert, on les établit à l’aide 
des potentiomètres. Pour faciliter la mise au point des coefficients, 
toutes les entrées des potentiomètres ont été disposées dans la partie 
gauche inférieure du panneau de commutation; ceci permet d'établir 
les coefficients sans rompre les liaisons organisées sur le champ de 
composition. Les douilles numérotées 2, 6, 10, . .., 46 permettent 
d'afficher un coefficient de transfert quelconque compris entre 0 
et 10, les douilles 4, 8, 12, ..., 48, un coefficient de transfert 
compris entre O0 et f. 

Supposons que l’on veuille attribuer le coefficient 6 à l’intégra- 
teur [5]. 11 faut pour cela utiliser l'entrée (18) et appliquer à la douille 
correspondante une tension étalonnée de 10 volts. On mesure la ten- 
sion de sortie de l’intégrateur [5] avec un voltmètre et on déplace 
le curseur du potentiomètre à l’aide d'un tournevis placé dans la 
fente jusqu’à ce que la tension régulée soit égale à —60 volts. Si 
l’on veut obtenir un coefficient de transfert égal à 0,6 il est alors 
plus rationnel de se servir de l'entrée (20) et d'appliquer à la douille 
correspondante non plus 10 mais 100 volts pour obtenir les mêmes 
—60 volts à la sortie de l'intégrateur [5]. 

Organiser les entrées, les circuits de réaction et les connexions 
entre les opérateurs d’après le schéma de commutation de la figu- 
re 3.6, a, afficher les coefficients de transfert &, B, y correspondant 
à la variante choisie dans le tableau 3.4. Si le schéma est correctement 
monté et les coefficients de transfert correctement mis au point, le 
voltmètre indique une tension nulle même sur l'échelle 2.5 V. 
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La précision de la mise au point des coefficients de transfert 
dépend de celle de la mesure des tensions électriques. Cette précision 
peut être accrue par l'application de la méthode de compensation. 
Soit mes la tension électrique à mesurer (cf. fig. 3.6, b). Cette tension 
est appliquée à l'entrée +V1 du voltmètre V1. Le voltmètre V1 
est mis au mode KOMITEHCALHA (COMPENSATION) par action 
sur le tumbler T13, puis son entrée —V1 est raccordée à la sortie de 
TENSION ÉTALONNÉE. La source de tension étalonnée doit 
délivrer une tension électrique 4. de même signe que uns. En action- 
nant les commutateurs à décades on amène l’aiguille du voltmètre 
Vi au zéro d’abord sur l'échelle 100 volts, puis sur l'échelle 2,5 volts 
et enfin sur l'échelle 0,1 volt. Ce qui traduit l'égalité u,jes = Ut. 
La valeur numérique deusest déterminée par la position du commu- 
tateur TENSION ÉTALONNÉE. 

7. Comportement de l'opérateur 15] en dehors de la plage +100 volts. 
La tension de sortie des AO équipant le calculateur MH-7M varie 
entre — 100 volts et +100 volts. En dehors de cet intervalle, les 


U 


Fig. 3.7 


calculs sont entachés d’erreurs supplémentaires dues à des pertur- 
bations non linéaires se manifestant dans le circuit électronique. 
La sortie de la variable machine de la plage +100 volts est la con- 
séquence d’une erreur dans le choix des échelles. L'opérateur humain 
contrôle la plage de fonctionnement de chaque AO à l’aide des lampes 
à néon branchées à la sortie de ce dernier. Les lampes s'allument 
dès que la tension de sortie des AO est de +105 volts. Les 18 lampes 
au néon L15 sont disposées dans le coin supérieur gauche du panneau 
avant du bloc principal (cf. fig. 3.2). Elles portent l'inscription 
CHTHAJIH3ALHNA TEPETPY30OK (SIGNALISATION SUR- 
CHARGE). 

On se propose de déterminer la valeur de la tension de sortie de 
l’intégrateur [5] à l’instant précis où la lampe s'allume. On se servira 
du schéma de la figure 3.7, a. Appliquons une tension de +100 V 
à l'entrée (20) de l’intégrateur [5] par un potentiomètre de coefficient 
de transfert 0,1. 
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L'intégrateur délivre la variable machine u (t) = u, — 104, 
où u, est la condition initiale. 

Quelle que soit la condition initiale u,, il arrive un instant 
t — t* tel que | u (t) | >> 100 volts (cf. fig. 3.7, b). Si u, = 0 cette 
inégalite est réalisée au bout de 10 secondes. 

A l'instant où la lampe de signalisation s'allume on ne peut 
mesurer la tension de sortie de l’intégrateur [5] avec V1 ou V2 puis- 
qu’elle est supérieure à 100 volts. La variable uv (t) est appliquée à 
l'entrée (4) de l'amplificateur d'échelle [1] par un potentiomètre de 
coefficient de transfert 0,5, de sorte qu’à la sortie de [1] on mesure la 
variable —u (t)/2 avec Vi. 


La donnée des conditions initiales de tous les intégrateurs et, 
en particulier, de la condition initiale u, de [5] s'effectue dans la 
position PABOTA (CALCUL) des tumblers T4 et T5. Les douilles 
ŸY5, V6, Y7, Y8, Y15, Y16 des conditions initiales des intégrateurs 
[5], (6], [7], [8], [15], [16] disposées dans la partie inférieure droite 
du panneau de commutation sont reliées par des cavaliers ou des 
cordons aux douilles de sortie 1, II, III, IV, V, VI des potentio- 
mètres d'affichage des conditions initiales. Les sorties des inté- 
grateurs sont connectées à tour de rôle par des cordons extérieurs à 
l'un des voltmètres. En déplaçant les curseurs des potentiomètres 
à l’aide des boutons P16 HAUAJIBHEIE YCJIOBHA (CONDI- 
TIONS INITIALES), on obtient la valeur désirée de la tension initiale 
à la sortie des intégrateurs. Les six boutons P16 sont disposés dans 
le coin supérieur droit du bloc principal. Les manettes des poten- 
tiomètres portent les mêmes chiffres romains Ï, 1], 111, IV, V, VI 
que les sorties des potentiomètres. Sous chaque bouton P16 se trouve 
un tumbler T17 pour la donnée du signe de la condition initiale. 

Dans le schéma de la figure 3.7, a les sorties de l’intégrateur (5] 
et de l’amplificateur d’échelle [1] sont reliées aux douilles B1 et B2 
des entrées de déviation verticale de l'indicateur électronique 
qui sont situées dans la partie inférieure du panneau de commutation 
a côté et légèrement à droite des douilles BX (ENT) et BbIX (SOR) 
du diviseur de tension /] 

Les variables u (t) et —u (t)/2 sont observées sur l’écran de l’in- 
dicateur. Celui-ci est branché lorsque le tumbler supérieur droit se 
trouve dans la position SECTEUR. La lampe de signalisation placée 
au centre du panneau, au-dessus de l'écran, s’allume ; au bout d’une 
ou de deux minutes on amène le tumbler supérieur gauche dans la 
position JIYY (FAISCEAU). Pour résoudre un problème en une 
seule fois, il faut mettre le tumbler T11 dans la position OHHO0- 
KPATH. (MONOCOURSE). Il faut veiller à ce que le commutateur 
PEXHM (MODE) occupe la position MONOCOURSE. Pour observer 
simultanément sur l'écran de l’indicateur deux variables machine 
appliquées aux entrées Bi et B2, on utilise le tumbler KOMMYTA- 
TOP (COMMUTATEUR) qui se trouve en bas du panneau de l’indi- 
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cateur. Si le tumbler occupe l'autre position, seule la variable appli- 
quée sur B1 est visualisée sur l'écran. 

Lorsqu'on observe une variable machine, la durée de balayage 
borizontal doit être légèrement supérieure à la durée de résolution 
du problème. Pour modifier la durée de balayage on dispose de deux 
interrupteurs BPEMfA PA3BEPTKH (DURÉE DE BALAYAGE): 
celui de gauche opérant par bonds, celui de droite, de façon continue. 

Composer le schéma de commutation (fig. 3.7, a) et établir les 
coefficients de transfert requis (T4 — CALCUL, T5 — PRÉPARA- 
TION). Brancher l'indicateur, le mettre au mode du bloc principal, 
brancher le commutateur. Faire varier les conditions initiales de 
l'intégrateur [5] (T4 — CALCUL, T5 — CALCUL), observer les 
variations sur l'écran et sur le voltmètre V1. Régler la durée de 
balayage en procédant à quelques essais de mise en marche de Ia 
machine par action sur le bouton K6 — IIYVCK (MARCHE). Faire 
Uuo = 0, us — +50 volts, u, — —50 volts, déterminer la tension 
d'allumage de la lampe au néon. Visualiser les courbes uv ({) et 
—u (t)/2 sur l'écran. 

8. Etude du circuit de commande préprogrammée (CCP) en mode 
d'arrêt. Le circuit de commande préprogrammée est organisé sur 


Circuté de commande 
préprogrammeé 


Fig. 3.8 


la base du 17-ème amplificateur opérationnel. Pour mettre cet ampli- 
ficateur en circuit on amène le tumbler STAJIOHHOE HAIIPAXKE- 
HE — IPOTPAMMHEIM PEKHM (TENSION ÉTALONNÉE — 
MODE PROGRAMMÉ) dans la position MODE PROGRAMMÉ. 
Le tumbler OCTAHOB — IEPEKJIIOUEHHE (ARRÊT — COM- 
MUTATION) qui se trouve à côté et à droite est mis dans la position 
ARRET. Le CCP stoppe alors automatiquement la machine à l’aide 
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des relais à contacts (qui en sont une partie intégrante) à l'instant 
où la variable appliquée aux entrées BX1 et BX2 du CCP change 
de signe. Avec le tumbler « + — » qui se trouve au-dessus du tumbler 
T11 sur le panneau avant du bloc principal on peut changer la pola- 
rité de la diode dans le circuit de réaction de l’AO [17] (cf. fig. 3.8). 
On distingue deux cas de changement de signe de la variable d'entrée 
du CCP : la variable s’annule et devient négative, la variable s'annule 
et devient positive. 

Le schéma de la figure 3.8 stoppe automatiquement la machine 
à un instant {* numériquement égal à la condition initiale À de 
l'intégrateur [5]. En effet (cf. graphe), la variable u, (t) — À —t 
devient négative à partir de {* — A. 

Composer le schéma de commutation de la figure 3.8, relier la 
sortie de [5] à l'entrée BX1 du CCP. Observer sur l’indicateur la 
variable u, (t) pour différentes valeurs numériques de la condition 


Tableau 3.5 


3 


{; Somme 


â E 


1EEIRINI 
CIRE 
| 1 | | 
tt 


initiale us (0) = À; (elles sont indiquées dans la deuxième ligne 
du tableau 3.5). Mesurer avec un chronomètre le temps t{,; écoulé 
entre la mise en marche et l'arrêt automatique de la machine. Com- 
pléter le tableau 3.5. 

Calculer les erreurs quadratiques moyennes avec la formule 


(4j —t})° | 


(1) 
1 
O — 5 D (4; —t}. 


j=1 


Enumérer et expliquer les causes éventuelles de l'écart entre les 
valeurs À; et ft}. 

9. Etude de la précision de la solution machine d'une équation 
différentielle. La figure 3.9, a représente un schéma de commutation 
comprenant un intégrateur [7] soumis à une réaction à l'entrée (28) 
par l’intermédiaire d’un potentiomètre de coefficient de transfert 0,3. 
On s’assure sans peine que la variable machine u, (f) vérifie l’équa- 
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tion différentielle 
du 


dont la solution analytique est 
u, (t) = 100 exp (—0,3f). 


La courbe intégrale est représentée sur la figure 3.9, b. On se propose 
de comparer la solution exacte u, (t) avec la solution machine. Les 


uy(t) 


dv./ 


f00exp (-036) 


à) b) 


Fig. 3.9 


valeurs exactes de la solution u,; correspondant aux valeurs discrètes 
suivantes 0, 1, 2, , 9 de la variable indépendante t, sont repro- 
duites dans la deuxième ligne du tableau 3.6. 


Tableau 3.6 


Monter le schéma de commutation de la figure 3.9, a. Mesurer 
avec le voltmètre V2 les valeurs v,; correspondant aux valeurs de f; 
indiquées dans la première ligne du tableau. Calculer l'erreur qua- 
dratique moyenne avec la formule 


o— _ : (ua, HE on) 


6—0112 
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Enumérer et expliquer les causes éventuelles des écarts entre u,; 
et Uzj. 

Il est recommandé d'utiliser le procédé suivant pour stopper la 
machine MH-7M au bout d’une seconde pour mesurer les variables 
machine. Presser le bouton K6 — IIVCK (MARCHE). Presser le 
bouton K7 — OCTAHOB (ARRÊT) dès que la lampe de signalisa- 
tion L9 s'allume. Le schéma de commande du calculateur MH-7M 
est conçu de telle sorte que la machine s’arrête au bout d’une seconde. 
Les valeurs numériques prises par les variables au bout d’une seconde 
sont mésurées avec un .:voltmètre. On remet la machine en marche, 
on l’arrête en pressant sur-le bouton K7 dès que la lampe LS s'allume 
puis on mesure les valeurs des variables machine pour t—=28s,et 
ainsi de suite: s 

10. Jntégration d'une équation finie avec un circuit de commande 
préprogrammée. Le schéma de commutation de la figure 3.10, a 
permet de trouver la racine .de l’équation 


| 100 exp (—0,1t) = t. 
En effet, l'intégrateur [8] résout l'équation différentielle 


8 — _Oius, us(0)—100 V. 


dt 
La ue analytique est ug (t) — 100 exp (—0,1t).: L'intégrateur 
[5] délivre la variable u, (t) — —t. L'addition des variables u, (t) 


2 
-#(5] , 
2 CCP 
Hbexpt01t\\6] 
a) b) 
Fig. 3.10 


et u, (t) provoque la réponse du CCP et l’arrêt de la machine à l’in- 
stant où u, (t) + us (t) = 0, i.e. lorsque —t + 100 exp (—0,1t) = 0. 
La valeur de la racine de l'équation est mesurée avec le voltmètre V2. 
Les variables u, (t) et u, (f) sont représentées sur l'écran comme l’in- 
dique la figure 3.10, 6. 
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7 Monter le schéma de commutation de la figure 3.10, a et déter- 
miner la valeur numérique de Ja racine. Vérifier si la solution est 
exacte en utilisant la table de la fonction exp (—t). 

11. Etude de la précision de la multiplication. Composer à cet 
effet le schéma de commutation de la figure 3.11, a. L’AO [2] du 
multiplieur B-II délivre la variable -u, (t) = 10-?ui (t); l’AO [4] du 
multiplieur B-IV la variable u,(t) — 10-tu (t). Les variables 
us, (t), us (1), u, (t) sont reproduites sur l'écran (cf. fig. 3.11, b) pour 


J 

y (t) B-N À 

(61-2100 AS, 0° 
(9) D : DT L 


Fig. 3.11. 


faciliter le contrôle visuel du fonctionnement du schéma. Le CCP 
permet de limiter à 20 secondes la durée de reproduction des variables. 
Stopper et remettre la machine en marche toutes les secondes, mesu- 
rer et porter dans le tableau 3.7 les valeurs des variables v,,, L.;, vyy 
correspondant aux valeurs {; indiquées dans la première ligne. Les 
lignes 2, 5, 8 reproduisent les valeurs exactes (théoriques) des 
variables us, ue, ua. Calculer les erreurs quadratiques moyennes 
avec les formules 


TU TU 
1 14 

N V + 2 (ue vas 0= + Di (es — vas), 
j=i j=1 


Enumérer et expliquer les causes éventuelles des erreurs. 

Le multiplieur-diviseur (MD) du calculateur MH-7M est commuté 
sur le panneau de commutation. Il utilise deux amplificateurs opé- 
ç® 


L'E nvalqv EE 
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rationnels supplémentaires. Ainsi pour organiser le multiplieur 
B-IV (fig. 3.11, a et 3.12) il faut effectuer les connexions indiquées 
sur la figure 3.13, a. Le multiplieur B-IV.est construit sur la base 
des AO [3] et [4]. Le panneau de commutation (cf. fig. 3.3) comporte 


6-1 


La : 
DE — 1 


Fig. 3.12 


un rectangle CETKH, BBIXOJIDI, B-III, B-IV :(GRILLES, SOR- 
TIES, B-I1II, B-IV) avec 4 paires de douilles. Les MD sont disposés 
en bas du bloc principal sous le panneau de commande. Ils font 
bloc avec les générateurs de fonctions non linéaires: ]ls portent des 
numéros pairs et les générateurs de fonctions, des numéros impairs. 
Le MD peut multiplier ou diviser. Pour choisir entre ces deux opé- 
rations, on dispose d'un commutateur ou de deux .tumblers placés 
à l’arrière du MD à côté des prises. Pour commuter Jes tumblers il 
faudra donc retirer le MD de la machine. Les entrées 1 et 2 du MD 
et ses sorties dédoublées (fig. 3.13, b) sont placées sous le champ de 


Fig. 3.13 


composition dans un rectangle portant l'inscription B-1V. Par ailleurs 
toute sortie de l’AO [4] est sortie du B-IV. Le MD est tel que l’on 
peut attribuer n'importe quel signe au produit effectué. Le signe 
est affecté par commutation des douilles sur la face du MD. 

Pour diviser, on envoie le dividende à l'entrée (1) du MD et le 
diviseur à l'entrée (2). 

12. Etude de la division. Commuter le bloc B-IV sur la division 
(cf. fig. 3.11, a). L'échelle du diviseur étant égale à 10, l'AO [4] 


86 TRAVAUX PRATIQUES [CH. 3 


doit théoriquement délivrer la variable 


1Our(t) 
Uy () = oo, y — lus (E). 

Visualiser les variations de u4, u+, u, en changeant alternative- 
ment (par commutation des opérateurs) les signes du dividende et 
du diviseur. Interpréter les résultats obtenus et les comparer avec 
les théoriques. Formuler la règle du signe du quotient en fonction 
des signes du dividende et du diviseur. Interpréter le résultat de la 
division dans le cas où le diviseur est une variable machine com- 
prise entre —10 et +10 volts. 

13. Remise à zéro des amplificateurs. Mettre à zéro les amplifica- 
teurs opérationnels et expliquer la « dérive» des zéros. 


& 2. Technique de composition des schémas de commutation 
sur le calculateur MH-7M 


Les schémas de commutation sont câblés sur le panneau de com- 
mutation, machine branchée (T4 — RÉGLAGE ZÉROS, T5 — 
PRÉPARATION). On recommande la succession suivante. 

1) Organiser les entrées et les circuits de réaction des opérateurs 
sur le panneau de commutation comme l'indique le schéma de com- 
mutation. Il est conseillé de signaler sur le schéma l'organisation de 
chaque entrée pour éviter toute confusion. 

2) Organiser les connexions entre les opérateurs sur le champ de 
composition. Prendre les mêmes précautions que précédemment. 

3) Régler les zéros des 18 amplificateurs opérationnels. 

4) Etablir les coefficients de transfert exigés des opérateurs 
(T4 — CALCUL, T5 — PRÉPARATION). 

5) Organiser les connexions omises dans 2) et assurant l’alimen- 
tation en tensions constantes. 

6) Donner les conditions initiales des intégrateurs (T4 — CAL- 
CUL, T5 — CALCUL). 

7) Régler l'indicateur électronique: mode de travail, durée du 
balayage et amplification des entrées verticale et horizontale. 

Nature et succession des opérations. Composer trois schémas 
de commutation : le premier étant un des schémas des figures 3.14 
à 3.17 pour une variante quelconque du tableau des variantes; le 
second, une variante du schéma de la figure 3.18 pour. toutes les 
valeurs des paramètres des courbes; le troisième, une variante du 
schéma de la figure 3.19. 

Chaque schéma reproduit une courbe plane. La visualiser sur 
l'écran, puis la tracer à l'aide d'un traceur à deux coordonnées. 

Les schémas des figures 3.14 à 3.17 ne contiennent que des opé- 
rateurs linéaires, des sommateurs et des intégrateurs, et reproduisent 
des courbes appelées cycloides. 
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g{t){f1) 
(=) 
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cost [5] 
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Fig. 3e15 


On appelle hkypocycloïde la courbe décrite par un point d’un cercle 
roulant intérieurement, sans glisser, sur un autre cercle. 

On appelle épicycloïde la courbe décrite par un point d’un cercle 
roulant extérieurement, sans glisser, sur un autre cercle. 
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Fig. 3.17 


On appelle hkypotrochoïde allongée la courbe décrite par un point 
situé en dehors d’un cercle roulant intérieurement, sans glisser, sur 
un autre cercle. 


Paramètres des courbes 
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Fig. 3.18 
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Fig. 3.19 
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On appelle épitrochoïde allongée la courbe décrite par un point 
situé en dehors d'un cercle roulant extérieurement, sans glisser, 
sur un autre cercle. 

Dans chacune de ces variantes, les cycloïdes possèdent un nombre 
fini de branches. Donc, si le schéma est correct et les coefficients de 
transfert assez exacts, le point décrivant la courbe revient à l’endroit 
d'où il est parti après avoir effectué un ou plusieurs tours du cercle 
fixe. Si les coefficients de transfert ne sont pas exacts, les cycloïdes 


><? EL 
LP 


Fig. 3.20 


ne reviennent pas à leur point de départ mais à un point situé légè- 
rement à côté, de sorte qu'après chaque tour de cercle elles sont 
décalées en avant ou en arrière. De tels cas sont représentés sur la 
figure 3.20. 

Le schéma de la figure 3.18 simule des courbes planes appelées 
courbes de Habenicht, du nom du mathématicien naturaliste alle- 
mand qui a cherché des expressions analytiques pour décrire de 
nombreux modèles géométriques du monde végétal. Le schéma com- 
prend un grand nombre de multiplieurs pour la réalisation du rayon 
vecteur r et de ses projections sur les axes de coordonnées x et y. 
Les « fleurs » de Habenicht citées dans le tableau des variantes 
possèdent un nombre entier de pétales. Donc, le point décrivant la 
fleur doit revenir à son point de départ si les coefficients de transfert 
sont correctement mis au point. Certaines courbes de Habenicht sont 
représentées sur la figure 3.21. Les valeurs numériques des coefficients 
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Fig. 3.21 


de transfert du sommateur [2] doivent être empruntées au tableau 
des paramètres des courbes, variantes 1 et 6. Toutes les autres varian- 
tes s’obtiennent par élimination et rétablissement de certaines liai- 


sons à l'entrée du sommateur {2] sans remise au point des coefficients 
de transfert. 
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Le schéma de la figure 3.19 renferme des opérateurs dont une 
partie sert au contrôle par la méthode de la variable redondante. 
Donc, si le schéma est juste et les coefficients de transfert correcte- 
ment mis au point, le sommateur [10] délivre une variable nulle. 


$ 3. Résolution de problèmes de Cauchy 


Nous allons étudier la résolution d’un problème de Cauchy *} 
en simulant le mouvement d’un point matériel lancé obliquement 
dans un milieu résistant. Dans le chapitre « Divers problèmes de 
simulation » on trouvera la position et les méthodes de simulation 
d’autres phénomènes physiques se ramenant à des problèmes de 
Cauchy. 

1. Position du problème. Un point matériel de masse m est lancé 
d’un point © dans une direction faisant un angle & avec l'horizontale, 
avec une vitesse initiale v, (fig. 3.22, a). Le point se déplace dans un 
milieu plan opposant une résistance Q proportionnelle à la vitesse 
du mouvement. On demande de simuler le mouvement du point 
matériel en supposant que le point O et le point de chute À sont situés 
à une même hauteur. 

2. Mise en équations du problème. Soit un système de coordonnées 
rectangulaires zOy dont l’axe des x passe par O et À. La loi de New- 
ton donne 


X=mx, Y=my, (3.1) 
où À, Ÿ sont les projections sur les axes de coordonnées de toutes 
les forces appliquées sur le point matériel, x et y les accélérations du 
point matériel le long des axes de coordonnées. 

Le point matériel est soumis à l'action de deux forces : son poids 
P — —mg et la résistance du milieu Q — —kmv, où k est un coef- 
ficient de proportionnalité. En portant les expressions des projec- 
tions des forces P et Q sur les axes de coordonnées 
X = Àp + Xo Y=Yp+ Yo 


« 


ou 
Xp=0, Xo=—kmr, Yp = mg, Yo = —kmy 


dans (3.1), on obtient le système d'équations différentielles du mou- 
vement 


re — kz, z (0) =vocosa, x(0)=0, 
é ; (3.2) 
y—= —ky—g, y(0)=vosina, y(0)—-0. 


*) Le chapitre 11 d'Eléments de simulation sur calculateurs analogiques 
contient un grand nombre de problèmes de Cauchy rattachés à la simulation de 
la cinétique de réactions chimiques. 
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Pour contrôler la machine, on appliquera la méthode de la variable 
redondante. Introduisons à cet effet la variable redondante 2 (t) 


Arrét programme 
a da condition 
y=0 


b) 


Fig. 3.22 


et exigeons qu'elle vérifie la relation de contrôle 


so ke + x + ky + y —3=0. (3.3) 
Pour réaliser (3.3) il est nécessaire de chercher d'abord une méthode 
de! reproduction :de la.:variable.redondante z (t). En dérivant. (3.3) 
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par rapport à { on trouve aisément l'équation différentielle . LE 
minante | 


B= ke 2 Hg + ÿ. 
C.mpte tenu de (3.2) on obtient finalement les expressions tathé- 
matiques indispensables à la simulation du problème: 


2= — g, z(0)=tv% (C0s +5 &); 
T — LT, z (0) =vcos a, “ x(0)= 0: ". ; (3.4) 
y —hy—g, “y(0)=vésie, y(0)=0 “| 


pas 


3. Schéma fonctionnel. Les schémas fonctionnels sont élaborés 
pour les variables mathématiques ct indiquent le principe de con- 
nexion des opérateurs pour câbler les expressions ‘mathématiques 
dorinées. Le-schéma fonctionnel réalisant (3.4) est celui de la figu- 
re 3.22, b. É 

Le schéma prévoit un arrêt de la machine à la condition y = 0, 
ce qu: traduit le fait que le point de lficement O et le point de 
chute À sont situés à la même hautéür<-Les Sorties des intégrateurs 
délivrant les variables y (t) et x (t) sont reliées aux entrées de dévia- 
tion verticale et horizontale de l'indicateur électronique. La tra- 
jectoire y —  (r) du mouvement est reproduite sur l'écran de l'in- 
dicateur. À la sortie du sommateur de EC qui élabore l’« erreur » 
de la relation de contrôle 


= — (kr + 2 + ky + y— 2) 


est branché un voltmètre pour |’ observation visuelle de cette 
« erreur ». 

4. Choix des échelles. Les échelles sont des coefficients de propor- 
tionnalité dimensionnels entre les variables mathématiques et les 
variables machine. Il y a lieu de faire une distinction entre le choix 
des échelles des variables dépendantes et de la variable indépendante. 
On ne modifie l'échelle que des variables mathématiques dépendan- 
tes qui figurent « explicitement » dans le schéma fonctionnel. Autre- 
ment dit, les variables dont on modifie l'échelle sont élaborées par 
les opérateurs comprenant un amplificateur opérationnel. Ici il 


s'agira des variables x (t), x (t), y (t), y (), z(t). Il n'est pas néces- 


saire de choisir les échelles des variables z (6) et y "(#) puisqu'elles 
ne sont délivrées par aucun opérateur contenant un amplificateur 
opérationnel. On a le tablean de correspondance suivant entre les 
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variables mathématiques et les tensions électriques . . 


V'ariable Tension 


rt) | Uxy = Moz (1) 
z () l Uxi = MZ (#) 
y (t) .' Uyo = Moy (4) 


y (t) Uri — my (4) 
z (4) | U, = M2 () 


Uxos Uxis Uyor Uyrs U; SOnt des tensions électriques, ne. ma, Mo 
les échelles. Pour faciliter la comparaison, les échelles des variables 
de même dimension ont été prises deux à deux égales. 

Les échelles doivent être telles que les variables machine restent 
dans la plage +100 V (pour les calculateurs à. ‘ampes). Cette con- 
dition sera respectée si l’on calcule l'échelle d'une variable mathé- 
matique, par exemple, z(t), avec la formule 

| m0 gs» 


| | Zmax | ° 
En effet, il vient alors | u, | — m,|2z(t) | << 100 V. 
Dans notre Cas, ROUE calculer les échelles il faut connaître 


Tmax;: his Ymax» TR ns 
Où les prendre puisque la solution analytique de (3.2) est inconnue? 
I1 nous fäut dônc chercher les valeurs maximales des variables sans 
tecourir à à la solution analytique de (3.2). On faït alors appel au 
principe physique du processus. Au lieu du processus donné, décrit 
par .des équations différentielles « compliquées », on choisit un pro- 
cessus similaire, régi par des équations simples dont la solution 
analytique permet d'évaluer la plage de variation des variables 
däns le processus simulé. Ce processus auxiliaire s'appelle problème 
d'échelle. Dans notre cas, comme problème d'échelle on peut considé- 
rer le mouvement d'un point matériel lancé obliquemént dans le 
vide. Le vide n'offrant pas de résistance, le point matériel se dépla- 
cera plus vite et tombera plus loin. Donc les valeurs maximales des 


vitesses x, y et des déplacements x, y du point dans le vide convien- 
dront pour le choix des échelles des variables décrivant le mouve- 
ment d'un point dans un milieu résistant. Au mouvement dans le 
vide correspond la valeur k = 0, donc le système d'équations diffe- 
rentielles (3.2) devient 


z= 0, z (0) =VpCOS&, ZT (0) = 0; 
.. (3.6) 
y= —g, y(0)=wsina, y(0)=0. 
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Les équations (3.6) admettent la solution analytique simple 


3 
z(t)=(vocosaæ)t, y(t)=(vosin a)t—Æ, 
d'où il suit que: 

1) quel que soit l’angle de lancement «, les projections z, y 
de la vitesse v sur les axes de coordonnées x et y ne sont pas supé- 
rieures à la vitesse v,; 

2) le point de chute le plus éloigné zmax Correspond à un angle 
de lancement &« = x/4 et se trouve à une distance de vi/g; 


3) la hauteur maximale ymax atteinte par le point matériel cor- 


respond à un angle de lancement « = < et vaut vy/2g; 


4) la valeur maximale de Zzmax << 2Uo. 
Le calcul des échelles est conduit dans le tableau 3.8. 


Tableau 3.8 


Valeur maximale Symbole de 
Variable de la variable Echelle l'échelle 


95. Schéma de commutation. Le schéma de commutation est 
composé pour les variables machine et orienté vers le calculateur 
retenu : ici le calculateur MH-7M. La figure 3.23 représente des sché- 
mas de commutation. À la différence du schéma fonctionnel de la 
figure 3.22, ils contiennent des potentiomètres supplémentaires aux 
entrées des opérateurs [6], [7], [8], [9], [15], qui servent à l’adapta- 
tion des échelles. 

Le circuit de commande préprogrammée du schéma de commuta- 
tion de la figure 3.23, b interrompt la résolution aussi bien à la con- 


dition y = 0 qu'à la condition y = 0, ce qui correspond à l'instant 
où le point matériel atteint la hauteur maximale y... L'intégrateur 
116] est réservé à la mesure du temps de vol du point matériel. 
Signalons encore que dans les schémas de commutation de la 
figure 3.23, a et b, les coefficients de transfert des opérateurs et les 
conditions initiales des intégrateurs sont représentés génériquement 
sous forme d'expressions mathématiques renfermant les paramètres 
a, k, v, du processus simulé. C'est uùe circonstance extrêmement 
importante, car le schéma permet à l'aide d'une commutation de 
déterminer l'influence éxercée sur le:processus par: les divers para- 
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mètres soit séparément. soit ensemble par variation des coefficients 
de transfert des opérateurs et des conditions initiales des intégrateurs. 


-{00cos & 0 [6] 


Arrét programmé à 
la condition 


Arrét programmé à 
a condition 


y=0 ou ÿ=0 


Fig. 3.23 


Le schéma de commutation générique s'obtient par le calcul des 
échelles à l’aide de la méthode des problèmes d'échelle. 

Si la plage de variation des paramètres du processus simulé 
est fixée, il faut concrétiser le schéma de commutation, c'est-à-dire 
déterminer les numéros des entrées utilisées des opérateurs. La 
7—0112 
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figure 3.23, b représente un tel schéma de commutation. On suppose 
que O0 <L'k < 1571, 0 < 2g/vo K 1 s71 

Attardons-nous sur le choix de l'échelle de la variable indepen- 
dante qui représente le temps dans les calculateurs analogiques. Ce 
choix est lié à. la durée de la résolution du problème par la machine. 
L'échelle temporelle est automatiquement adaptée à celle des varia- 
bles mathématiques dépendantes, de sorte que l’unité de temps est 
la seconde machine. Un cas important sur le plan terminologique 
est celui où la variable indépendante du processus simulé est le 
temps. Ceci étant, si à une seconde de marche du processus réel cor- 
respond une seconde machine, on dit que le processus est simulé 
en {emps réel. On dit par ailleurs qu'il est simulé en {emps accéléré ou 
ralenti selon qu'il est réalisé plus ou moins vite qu'en temps réel 
par la machine. La notion de temps accéléré, ralenti ou réel n’a bien 
sûr pas de sens pour les phénomènes dans lesquels la variable indé- 
pendante n’est pas le temps, mais par exemple la longueur, comme 
dans les problèmes de flexion d’une poutre. Donc, dans le cas général, 
ce n'est pas l'échelle du temps qui est essentielle, mais la durée de 
résolution du problème sur la machine. La durée maximale de réso- 
lution d’un problème par la machine MH-7M est limitée par les 
erreurs d'intégration et est comprise entre 150 et 200 secondes. Il 
faut compter 10 à 20 secondes pour une observation visuelle conve- 
nable. A signaler que la variation de la durée de réalisation du 
processus simulé ne touche pas les étapes précédentes et obéit à la 
règle suivante: 

pour modifier de m fois la durée de résolution d'un problème il faut 
modifier de m”! fois les coefficients de transfert de tous les intégrateurs 
et d'eux seuls *). 

Cette règle est une conséquence directe du fait que les variables 
de sortie et d'entrée de l’intégrateursontentre elles comme une varia- 
ble et sa vitesse de variation. Une modification identique et simul- 
tanée des coefficients de tous les intégrateurs entraîne une variation 
de la vitesse de reproduction de toutes les variables machine, et 
partant de la durée de résolution du problème. 

6. Marche à suivre. 1) Dans le tableau 3.9 repérer les valeurs 
numériques des paramètres du processus simule de la variante donnée. 

2) Simuler le mouvement du point matériel pour les paramètres 


Ld r1 4 « # e e 
donnés et a — -- d'après le schéma de commutation de la figu- 


re 3.23, b. Si le schéma de commutation est correctement monté 
et les coefficients de transfert des opérateurs exactement mis au 
point, l'erreur de la relation de contrôle n'est pas supérieure à 
+0,5 V. Si l'erreur excède cette quantité, il faut la localiser et en 
éliminer la cause. 


*) On suppose que le schéma n'utilise pas de dérivateurs. 
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Déterminer le temps {, pendant lequel le point est en l’air et les 
composantes x (f,) et y (f,) de la vitesse de chute. Si t, < 10 5, 


Tableau 3.9 


10 


FI 


Î 


to 


Paramètre 


Æfs-1110,3 [0,3 [0,3 |0,2 |0,2 |0,2 [0,4 |0.4 |o0,4 | 0.1 
0 


k{s”1] | 0,05 | 0,10 | 0,15 | 0,03 | 0,06 


= 
8 
tS 
© 
© 
= 
© 


modifier la durée de résolution du problème en la portant à 10 s. 
Relever les valeurs de x ({) et y (t{) toutes les secondes et les porter 
dans le tableau 3.10. 


Tableau 3.1V 


Construire le graphe y = œ (x). 
PL‘ 
4 
k compris entre O0 et 0,4 avec un pas de 0,05. Figer pour chaque tra- 
jectoire les valeurs des variables: r,,,, — l'abscisse de la chute, 
Ymax — la hauteur maximale, x,, —- l’abscisse correspondant à la 
hauteur maximale, {,, — le temps mis par le point pour atteindre 
la hauteur maximale, {, — durée totale du vol. Porter les résultats 
des mesures dans le tableau 3.11. 
Interpréter les résultats obtenus et construire les graphes 


3) Poser &« = — , réaliser les trajectoires du point matériel pour 


Ymax = /1(Zmax), ln = fe (À). 


4) Poser k = O0, réaliser les trajectoires du point matériel pour 
divers angles de lancement compris entre n/12 et x/2 avec un pas 
de x/12. Porter ces résultats dans le tableau 3.12. Interpréter les 
résultats obtenus. Construire le graphe de Ymax — fa (Zmax)- 


7* 
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Tableau 3.11 
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7. Généralisation. Chute d’un point matériel sur une surface 
élastique. Supposons que le point matériel tombe sur une surface 
élastique schématiquement représentée par des plates-formes impon- 
dérables non reliées entre elles, reposant sur des ressorts (fig. 3.24, a). 
En tombant sur une telle plate-forme, le point est soumis à l’action 
répulsive de la force élastique F# — —mcy. Le point rebondit et 


-S0/min{0, Foÿl 


Fig. 3.24 


poursuit son mouvement. On demande de simuler ce phénomène. 
Précisons tout d’abord le système d'équations différentielles (3.2). 
La force F agit vers le haut et uniquement lorsque le point se trouve 
au-dessous de l'horizon, i.e. lorsque y << 0, de sorte que 


0 pour y >0, 
—mcy pour y< O0, 


> 
] 
qe, 
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où m est la masse du point, c la rigidité du ressort. Le système d’équa- 
tions différentielles (3.2) s'écrit alors 


T= — kr, z(0)=0, z(0)=vocos a, G7 
25 x “ d./ 
y= —ky—g+Fim, y(0)=0, y(0)=v,sine. 


Le terme F/m est facilement càblé sur la machine avec une opéra- 
tion de choix. En effet, F/m — —c min {0, y}; les deux variantes 
de schémas fonctionnels sont représentées sur la figure 3.24, b. 
Le schéma de commutation est représenté sur la figure 3.24, c. 
La rigidité de la surface élastique vaut c — 500 c”, où c’ est le coef- 
ficient de transfert à la (26)-ième entrée de l’intégrateur [7]. 
Simuler le mouvement du point pour & = x/3, k — 0,1 et 0,10. 
Interpréter les résultats obtenus et tracer les graphes y = @ (x), 


y = # (x), F = We (x). 


$ 4. Réduction d’un problème aux limites 
à un problème de Cauchy 


Les problèmes aux limites sont envisagés pour des équations 
différentielles d'ordre supérieur au premier. Comme exemple de 
problème aux limites citons le problème 


y (D =q{t), y (0) = yo y (1) = y 


dont la solution est cherchée sur l'intervalle t € [0, 1], où (4) 
est une fonction donnée. Il est impossible d'obtenir immédiatement 
la solution de ce problème. Force est donc de lui trouver un problème 
de Cauchy équivalent, soit 


y") =œp(t), y (0) = yo y” (0) = y 


et de le simuler. 

D'une façon générale, un problème de Cauchy et un problème 
aux limites sont équivalents si leurs solutions sont confondues sur 
l'intervalle d'intégration. La recherche d’un problème de Cauchy 
équivalent à un problème aux limites donné s'appelle réduction du 
problème aux limites au problème de Cauchy. Nous allons considérer 
une méthode simple de réduction, appelée méthode des essais et 
corrections, sur l'exemple de la simulation du tir à l’obus sur une 
cible immobile. Dans le chapitre 6 des « Eléments de simulation sur 
calculateurs analogiques » on trouvera des problèmes aux limites 
rattachés à la simulation d'une flexion élastique et les méthodes 
de leur réduction à un problème de Cauchy. 
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Tir sur une cible à terre. 


1. Position du problème. Un canon placé en un point © d’une 
région défendue (RD) tire à l'instant { — 0 sur les positions enne- 
mies PE (fig. 3.25) se trouvant à une distance / de O. A sa sortie 
du fût, l’obus est animé d’une vitesse v,. On admet que l'air oppose 


ÿ 


Fig. 3.25 


une résistance proportionnelle à la vitesse de l’obus *). Simuler la 
trajectoire d’un obus qui atteint la cible. | 

2. Equations différentielles du mouvement. Le mouvement de 
l’obus est régi par le système d'équations différentielles déjà connu 
(3.4), à la seule différence qu’au lieu du problème de Cauchy on 
a affaire au problème aux limites 


zx = —kz, z(0)=0, z(t) =, 
y=hy—g, y(0)=0, y(#“*)=0, 


où {* est le temps de vol de l'obus, i.e. le temps compris entre la 
mise à feu et la destruction de la cible. 

Réduire le problème à un problème de Cauchy consiste dans 
le cas de l'équation (3.8) à résoudre le problème de rencontre de 
l’obus avec la cible. En d’autres termes il faut trouver un angle 
d’élévation à«* du canon tel qu’au bout d’un temps {* après la mise 
à feu l'on ait 


(3.8) 


z(*)=1 y(t*) = 0. 


3. Schéma fonctionnel et schéma de commutation. 
Ces schémas sont à peu de chose près ceux des figures 3.22 et 3.23. 
Nous les omettrons. 


*) Nous avons adopté une dépendance linéaire de la résistance par rappor 
à la vitesse du mouvement exceptionnellement pour simplifier le problème. 
En réalité cette dépendance est plus complexe. 
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4. Marche à suivre. 
1) Dans le tableau 3.13 repérer la variante donnée, calculer les 
coefficients de transfert du schéma de commutation et simuler la 


Tableau 3.13 


(5-11 10,014 10,014 10,020 | 0,020! 0,015 10,015 | 0,021 |0,016 | 0,018 | 0,016 
ù 
k[s”1] 0,006 |0,006|o.o0 0.10 | 0,005 0,00 |. | 0,004 | 0.008 | 0.00 


rencontre de l’obus avec la cible en admettant que celle-ci se trouve 
à une distance L — 0,4 v5/g de la pièce d'artillerie. 
2) Simuler la trajectoire de l’obus pour un angle d’élévation de 


= ; évaluer le temps de vol de l'obus; porter la durée de simu- 


lation des trajectoires jusqu'à 20 secondes machine. 

8) Réduire le problème aux limites (3.8) au problème de Cau- 
chy (3.4). A signaler que dans la plupart des cas il existe deux angles 
d'élévation, æÿ et «> qui permettent à l’obus de faire mouche. 


+ Ml 
Le Our Dame 


—V1 +! 


mpl|!] OL: 


Fig. 3.26 


La non-unicité du problème de Cauchy équivalent est une particu- 
larité importante du phénomène étudié. A chaque angle d'’élevation 
correspond une trajectoire de l’obus. La trajectoire du premier type 
est représentée sur la figure 3.25 par le chiffre romain I. Cette tra- 
jectoire est dite courbe. Le chiffre romain II représente la trajec- 
toire du second type, dite trajectoire rasante. À remarquer que 
lorsqu'on tire dans le vide les angles af et a> vérifient la relation 
aÿ + as = n/2. 
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La réduction s'opère par la méthode des essais et corrections. 
On se donne plusieurs valeurs arbitraires de l'angle &, soit 


Œjs Dep + ++ js (3.9) 


et l'on simule la trajectoire correspondant à chacune d'elles avec 
le schéma de la figure 3.23. Ce schéma comporte en plus l'inver- 
seur [1] (pour la donnée de la coordonnée x = ! au point de con- 
centration des forces ennemies) et le voltmètre VA est branché 
autrement. La figure 3.26 représente la partie du schéma de commu- 
tation notifiant ces modifications. Au cours des essais de simulation 


8 ê=-h(æ) 
4 


Fig. 3.27 


des trajectoires, le voltmètre V1 mis au mode COMPENSATION 
est utilisé pour mesurer l'erreur Ô = m, (x — Î) entre / et l’abscisse 
du point de chute de l’obus. Le signe de l'erreur indique si le tir 
est court ou long. Aux valeurs a; correspondent les erreurs 


Ô;, Ô»; . 3 Ô;, . (3.10) 


On construit le graphe de Ô = (a) (fig. 3.27, a) en se servant 
de (3.9) et (3.10). Le zéro de la fonction 1 (œ) correspond à la valeur 
cherchée a&*. Pour réduire le nombre d'essais on peut utiliser la 
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méthode suivante. Les deux premières valeurs, &, et &., sont arbi- 
traires, quant à @&.,,elle est choisie en fonction des deux précédentes 
(fig. 3.27, b). Les points Z et 2 sont reliés par une droite jusqu’à son 
intersection avec l'axe des &. Déterminer de cette façon la valeur 


Table trigonométrique 


cos sin cos Decgrés 


Dcgrés | Degrés 


0 0,0000 1,0000 90 23 0,3907 0,9205 67 

1 0,0175 0,9998 89 24 0,4067 0,9135 66 

2 0,0349 0,9994 88 25 0,4226 0,9063 65 

3 0,0523 0,9986 87 

A 0,0698 0,9976 66 26 0,4384 0,8988 64 

6) 0,0872 0,9962 85 21 0,4540 0,8910 63 
28 0,4695 0,8829 62 
29 0,4848 0,8746 61 


6 0,1045 0,9945 84 
T 0,1219 0,9925 83 30 0,5000 0,8660 60 
8 0,1392 0,9903 82 
Q 0,1564 0,9877 81 31 0,5150 0,8572 99 
10 0,1736 0,9848 80 32 0,5299 0,8480 58 
33 0,5446 0,8387 91 
34 - 0,5592 0,8290 56 
0,5736 0,8192 59 


11 0,1908 0.9816 
12 0,2079 0,9781 
13 0,2250 0,9744 
14 0,2419 0,9703 
15 0,2588 0,9659 


=) =) =) =) 
1 © © 
a 


«1 
En © 
1 


36 0,5878 0,8090 54 
1 0,6018 0,1986 93 
38 0,6157 0,7880 92 


16 0,2756 0,9613 14 39 0,6293 0,77571 91 
17 0.2924 0,9563 13 40 0,6428 0,766 ) 50 
18 0,3090 0,9511 72 
19 0,3256 0,9455 11 41 0,6561 0,1547 49 
20 0,3420 0,9397 70 2 0,6691 0,7431 48 
43 0,6820 0,7314 47 
21 0,3584 0,9336 69 44 0,6947 0,7193 46 
22 0,3746 0,9272 68 45 0,7071 0,7071 45 


suivante «, de l’angle d’élévation puis simuler la trajectoire cor- 
respondante; calculer l'erreur ô, à l’aide de V1 et représenter le 
point 3” correspondant. Mener par les points 2 et 3” une droite 
jusqu'à son intersection avec l'axe des & et trouver la valeur suivan- 
te «,, et ainsi de suite. 

Pour réduire le problème on se servira de la table des fonctions 
trigonométriques sinus et cosinus donnée plus haut. Faire la réduc- 
tion pour un tir courbe et un tir tendu. 
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4) S'assurer que l’obus a atteint la cible et vérifier le fonctionne- 
ment de la machine à l’aide d’un schéma de contrôle. Pour les deux 
types de trajectoires, définir le temps {? et {> de vol de l'obus; 
la vitesse v* et v5 de l'obus au moment où il atteint la cible; l’alti- 
tude maximale y, :nnx ©t Ye max de l’obus. Interpréter les résultats 
obtenus. 


$ 5. Simulation de l’équation de la chaleur 


1. Position du problème. Une paroi homogène, épaisse, de 
grande surface sépare deux milieux de température différente 
(fig. 3.28), la température du milieu de gauche varie comme @, (t), 
celle du milieu de droite comme ®,(t). L'épais- 
seur de la paroi est /. On demande de simu- 
ler la distribution de la température pour les 
points intérieurs de la paroi en considérant 
les fonctions œ, (ft) et œ. (t) données. 

Introduisons tout d’abord des variables 
mathématiques qui décrivent la distribution 
de la température en profondeur de la pa- 
roi. Traçons un axe Oy perpendiculaire à la pa- 
roi. La température des points intérieurs est. 
une fonction de deux variables: le temps t et 
la coordonnée y, i.e. v (f, y). Les points de la 
paroi situés sur une même verticale possèdent une même tempéra- 
ture (isothermes) par suite de la grande surface de la paroi et de 
son homogénéité. 

La fonction v (t, y) est solution de l'équation 


dv EX dr 


où À est le coefficient de conductibilité thermique, c la chaleur spé- 
cifique volumique de la paroi. La solution de l'équation (3.11) 
doit vérifier les conditions aux limites vw ({, 0) = œ, (t), v (t, !) = 
= ®p, (t) et la condition initiale v (0, y) = T (y). 

2. Approximation de l'équation de la chaleur. Pour chercher 
un système d'équations différentielles ordinaires approchant l'équa- 
tion (3.11), partageons la paroi par des lignes verticales en une 
série de couches de même épaisseur. A l'intersection des verticales 
avec l’axe des y, on obtient une série de points en lesquels la tempé- 
rature sera uniquement fonction du temps. On supposera que les 
points de l’axe Oy sont situés à la même petite distance À l'un de 
l’autre (fig. 3.29). Portons la coordonnée zu; du point j dans l’équa- 
tion (3.11): 

dv(t,yy) __ À dut, yj) 
Me or (9:22) 
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Si l’on désigne v(t,y;) par v,(t), on obtient 20.49 = 20 


dt 
Le e # e 9? {, e e A e 
La dérivée partielle ED peut approximativement être expri- 


mée en fonction des valeurs de v ({, y) au point y, et aux points 
voisins Yj+1 et Yy1: 


2 (4, 1 
0 R + [Uy4 (t)—2vi(t)+v;-:(e)], 


où 
Vs (E)=U(E, Yes), va (t)=v(t, y;-s). 
Compte tenu de ces expressions, l’équation 


(3.12) se transforme en l'équation différentielle 
ordinaire du premier ordre: 


O0 = pe {4 (4) — 2v, (8) + 05 (0), 
À 
v(0)=T(y)=Ts, p=—. 


On obtiendra autant de telles équations qu'il n'existe de points 
sur l'axe. Si les points sont au nombre de n, les équations formeront 
un système de n équations différentielles ordinaires qui approchera 
l'équation (3.12): 


0 = pH [U;+s (4) — 2v;, (t) + Uj-1 (£)], 


J= 1; 2:33: nn; 
v(0)=T;, Vo(t) = pi (t), Un+s (4) = Po (t). 
(3.13) 


Fig. 3.29 


Le système (3.13) est un |problèmeŸ'de 
Cauchy. 

3. Schéma fonctionnel. Le schéma fonc- 
tionnel résolvant le système (3.13) est re- Fig. 3.30 
présenté sur la figure 3.30. Ce schéma est écono- 
mique. Il nécessite un intégrateur en tout pour la résolution de 
chaque équation. Les rectangles indiquent les schémas générant 
les fonctions données ®, (t) et œ: (t). 

&. Adaptation des échelles des variables. En adaptant les échelles 
des variables on tiendra compte de ce que la température des points 
intérieurs ne peut dépasser les quantités œ, (ft) et œ (t). 

5. Marche à suivre*). Repérer la variante donnée sur le tableau 
3.14 et composer les schémas fonctionnel et de commutation pour 


*) Au chapitre 4 on trouvera divers cas de simulation de l'équation diffé- 
rentielle de la chaleur. 


$ 5] SIMULATION DE L'ÉQUATION DE LA CHALEUR 109 


la machine MH-7M (! = 1; toutes les données sont adaptées dans 
le tableau). 
Tableau 3.14 
Variante T (y) qe ({) 
sin Ty cxp(—t) 


1—sin y 


| | 1—exp(—1) | 
EL | 
| eXp(—1) | | 


1—cexp (—t) 


cxp(—{) 


y* 0 1—exp(—t) 


cos ñ1y exp(—t) 


Organiser le contrôle de la machine par la méthode de la variable 
redondante en utilisant la relation 


z(t)+ > vy(t)=0, (3.14) 


où z(t) est la variable redondante dont l'équation différentielle 
déterminante simple 


2 (t)= —Ivo(t)—vi(t) vus (t)— 0, (t)], 2 (0) = — À v; (0), 


1 — cos 77 | 1—exp (—t) 


se déduit de la relation de contrôle (3.14) par une dérivation et une 
substitution des dérivées v° (t) (j = 1,2, ..., n) du système (3.13). 


CHAPITRE 4 


DIVERS PROBLÈMES DE SIMULATION 


$ 1. Mouvement d’un point matériel 
sous l’action d’un centre attractif 


1. Description du phénomène. Un point matériel de masse m se 
déplace dans un milieu résistant sous l’action d’un centre attractif 
fixe O. Le point finira par tomber sur le point attractif, puisque son 
énergie diminue en raison de la résistance du milieu. La figure 4.1, a 
représente le point matériel et les forces agissant sur lui: P =—w"mr, 
force d'attraction du point vers le centre O, et Q — —kmw, résistan- 
ce du milieu. Outre ce cas, le schéma que nous donnons plus bas 
permet d'étudier le mouvement lorsque le point est soumis à l'action 
des forces P et Q et d'une force F 
(fig. 4.1, b) et lorsque le point est 
équipé en plus d’un moteur dont la 
traction $S = À;mu agit dans le sens 
de la vitesse et lui est proportion- 
nelle (fig. 4.1, c). 

2. Equations différentielles du mou- 
vement. En projetant toutes les forces 
agissant sur le point matériel sur les 
axes de coordonnées on obtient les 
systèmes d'équations différentielles 
suivants: 


pour la figure 4.1, a: Z = —w°r — kx, 
s . À (4.1) 
y = —@"y — ky, 
pour lafigure 4.1, biz — —w?r — 
s —kz + PF. (4.2) 
v #F . . 
m c) y= —oy —ky+F,, 
F pour la figure 4.1, CT = — Gr + 
0 e 7 LH —hzt+Er,, 


Fig. 4.1 y= —wy + —ky+r, (G5) 
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avec les conditions initiales communes 


z(0)=0, z(0) =, y(0) = Yw y (0) = 0, 


conditions qui supposent que le point matériel se trouve en x, de 
l’axe Ox à l'instant initial, ne possède pas de composante de Ja 
vitesse le long de l’axe Oz et que sa vitesse initiale le long de l'axe 


Oy est égale à yo. 

3. Schéma de commutation. Le schéma est représenté sur la 
figure 4.2. Î[l est composé pour les valeurs numériques suivantes 
des paramètres : £ = 0,1; #, — 2k; w — 1. Les projections F, = F, 
de la force F' qui agit le long de Ia bissectrice du premier quadrant 


Ent. F  I0V 6 CCP 
100V © Sort. 


Fig. 4.2 


sont assimilées à une tension électrique de 10 volts. Si la variable 
de sortie F' du diviseur D est appliquée aux entrées (35) et (43) 
des intégrateurs [5] et [7] et si les contacts 2PO et 4PO occupent 
la position indiquée par le schéma, on obtient la solution du systè- 
me (4.2). Si l’on change le signe de la tension appliquée à l'entrée 
du circuit de commande préprogrammée, les contacts 2PO et 4PO 
occupent la position inverse et la machine simulera le système 
d'équations (4.3). Les coefficients de transfert aux entrées (20) et (28) 
peuvent varier entre O et 1,1. La force constante F peut être modi- 
fiée à l’aide du diviseur de tension D. 

4. Résultats de la simulation. Supposons qu’on ait opéré une 
déconnexion aux entrées (35), (20), (28), (43); la machine résoudra 
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alors le système (4.1) pour 4 — 0. La trajectoire du point matériel 
est une ellipse. 

Si l'on rétablit les connexions aux entrées (20) et (28) et que 
l'on fasse varier progressivement la quantité k, on obtient sur l'écran 
une série de spirales à point asymptotique en l’origine des coordon- 
nées, où se trouve le centre attractif. À partir d’une certaine valeur 
de À dépendant de w, le point, sans avoir effectué un tour autour 
du centre attractif, commence à tendre asymptotiquement vers lui. 


Fig. 4.3 


Un tel mouvement est dit apériodique. Les trajectoires sont repré- 
sentées sur la figure 4.3. Si, ayant retenu une valeur de À (par exem- 
ple 0,1), on met la machine en marche et l’on attend que la trajectoire 
se rapproche du centre attractif, et qu'on modifie ensuite le signe 
à l'entrée du circuit de commande préprogrammée, on constate que 
la trajectoire, de spirale enroulante se transforme en spirale dérou- 
lante. ce qui traduit l'intervention de la force de traction S. 
L'établissement des connexions aux entrées (35) et (43) ne modi- 
fie en principe pas la nature du mouvement du point, mais déplace 
simplement la trajectoire en direction de la bissectrice du premier 
quadrant d’une distance qui dépend de la grandeur de la force F. 


$ 2. Pendule double 


1. Description des oscillations. Imaginons deux pendules appro- 
ximativement de même période d'oscillations, l’un lourd, par exem- 
ple un lustre, l’autre léger. Le pendule léger est suspendu au lourd 
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comme l’indique la figure 4.4. Si l’on écarte légè- 
rement le pendule lourd et qu'on le relâche, il ar- 
rête brusquement d'’osciller. Le pendule léger com- 
mence alors un mouvement rapide qui cesse au bout 
d'un certain temps; à ce moment le pendule lourd 
commence à osciller. De telles oscillations alterna- 
tives de deux systèmes, accompagnées d’un échange 
d'énergie s'appellent battements. 

Simulons ces battements. 

2. Equations différentielles du mouvement. Les 
oscillations du pendule double sont étudiées en Fig. 4.4 
physique et en mécanique. Nous ne nous arrêterons 
pas sur la description mathématique du phénomène et nous donne- 
rons les équations différentielles dans le cas de faibles oscillations. 
Ce processus est régi par le système d'équations différentielles 


X=—w2(1+2p)X+uo?z, X(0)=0, X(0)=Xo, 


z= —62r7+0%X, z(0)=z(0)=0, 
où X et x sont les écarts des pendules lourd et léger par rapport à la 


position d’équilibre, u = 7 un petit nombre, m et M les masses 


des pendules léger et lourd, w* = g/l, g l'accélération de la pesan- 
teur, ! la longueur du pendule léger. On suppose que ZL = L. 

3. Schéma simulant le système (4.4). Il est représenté sur la 
figure 4.5. Il correspond au cas w = 1, u — 0,01. Le processus est 
simulé en temps accéléré. A l'écart initial X (0) = X, du pendule 
lourd correspond une tension électrique de 10 volts (cf. conditions 
initiales de l’intégrateur [6]). Pour faciliter l'observation visuelle 
des oscillations, au lieu de la variable zx (t) on applique la variable 
X (4) = —25 + 0,1r(t) à l'entrée d.v.i de l'indicateur. Ceci 
permet de décaler verticalement les oscillogrammes X (ft) et x (t). 

4.Résultats de la simulation du système (4.4). Ils sont repré- 
sentés sur la figure 4.6. En haut on a les oscillogrammes des varia- 
bles X ({) et x (t) qui montrentquela valeur maximale de l’amplitu- 
de d’oscillations d’un pendule correspond à la valeur zéro de l’am- 
plitude de l’autre pendule. En bas on a le graphe de la solution x — 
—= % (X) du système (4.4). Il est aisé de voir que c’est une courbe 
de Lissajous qui couvre un domaine contenu dans un parallélo- 
gramme. 


(4.4) 


$ 3. Pendules liés 


1. Descriptions des oscillations. La figure 4.7 représente un système 
oscillant composé de deux pendules, dits liés ou sympathiques, reliés 
entre eux par une liaison mécanique élastique figurée par le ressort. 
Sans cette liaison les pendules oscilleraient indépendamment l’un 


8—0112 
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de l’autre. La liaison entraîne généralement un « transfert » d'énergie 

d’un pendule à l’autre, donc une modification de la forme des oscil- 

lations de chacun d'eux. On se propose de simuler ces oscillations. 
2. Equations différentielles. Les équations 

différentielles des faibles oscillations *) de 

pendules liés s'écrivent 


z: — — DT, —k; (z1— 2) —2urs, (4.5) 
2 = — Te — ko (T2 — Zi), 


où z, et x, sont les écarts des pendules par 
rapport à la position d'équilibre, w, et os: Fig. 4.7 
les fréquences déterminées par les paramètres 
des pendules, k, et k. des coefficients dépendant de la rigidité de la 
liaison entre les pendules; on suppose que k, = k., « un coefficient 
qui tient compte du frottement visqueux opposé aux oscillations 
du premier pendule. 

3. Schéma de commutation. Le schéma simulant le système 
(4.5) est représenté sur la figure 4.8. Les valeurs numériques des 


Fig. 4.8 


coefficients de transfert des opérateurs correspondent au cas de 
résonance Où. @y = @2 = 10, k, = k, = 1, à — 0,1. On visualise 
la variable X, — 25 + x, au lieu de la variable zx, (t), ceci dans 
le but de décaler verticalement les oscillogrammes des oscillations 
x, (t) et ze (t). 

4. Résultats de la simulation. Dans le cas de la figure 4.8, où 


*) Pour la déduction de ces équations, on pourra consulter un ouvrage de 
mécanique théorique. 


s # 
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Zi (0) = ze (0) = x, (0) = 0, x, (0) = 10 V les oscillations à & — 0 
(il faut rompre la liaison à l'entrée (20) de l’intégrateur [5]) sont 
des battements, les mêmes que pour le pendule double (cf. fig. 4.6). 
Si l’on rétablit la connexion à l'entrée (20) les oscillations seront 
amorties. 

Le comportement des pendules liés dépend essentiellement 


de leurs positions initiales. Si zx, (0) = x, (0) = 0 et z, (0) = 
= Z9 (0) Æ O0 ou zx, (0) = —zx, (0) -Æ 0, chaque pendule oscillera 
indépendamment de l’autre, i.e. il n’y aura pas d'échange d'énergie 
entre eux. Ces oscillations sont dites normales. 

Le schéma de la figure 4.8 permet d'étudier le comportement 
des pendules dans le cas d’un « désaccord » où w, et kX, sont légère- 
ment différents de w. et k.. Ces dernières quantités peuvent être 
modifiées à l’aide des potentiomètres placés aux entrées (26) et (28) 
de l'intégrateur [7]. 


$ 4. Oscillateur spatial 


1. Position du probième. Un oscillateur spatial anisotrope 
chargé est placé dans un champ magnétique homogène H ; on deman- 
de de simuler sa trajectoire. 

2. Equations différentielles du mouvement. On supposera que 
l’oscillateur possède une énergie potentielle u (x, y, z) = F [or + 
+ wËy? + w?z°)] et que le champ magnétique homogène est parallèle 
a l’axe des z. Si le potentiel vectoriel du champ est de la forme 
Z (y, —z, 0), le lagrangien s’écrira 


LR (a+ pe 2) 5 (0t22+ obyt+ os?) + DE (yz— ay), 


OÙ GO — À ,et e est la charge de l'oscillateur. Si l’on connaît 


le lagrangien on trouve sans peine le système d'équations différen- 
tielles cherché 


. H e 

TZ — — oiz+— y; 

.e e eH e 

Y= —0,y— TL (4.6) 
ee — (02Z. 


3. Schéma de commutation. Le schéma résolvant le système 
(4.6) est donné sur la figure 4.9. Les coefficients de transfert du sché- 
ma correspondent au cas isotrope où &, — ©, = &, — 1. Les 
valeurs w, et w, peuvent varier dans des limites assez larges aux 


$ 4) 


OSCILLATEUR SPATIAL 
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entrées (20) et (28) des intégrateurs [5] et [7]. La quantité ox — 


e # e a e e » 
ee 2, cependant sa valeur numérique peut être modifiée 


à l'aide des potentiomètres placés aux entrées (18) et (26) des intégra- 
teurs [5] et [7]. 

&. Résultats de la simulation. Les plus intéressantes sont les 
trajectoires *) situées dans le plan zOy. Elles dépendent essentielle- 
ment des conditions initiales, lesquelles définissent l’énergie initiale 
de l’oscillateur qui est la somme de l'énergie cinétique et de l’éner- 
gie potentielle. Quatre trajectoires de l'oscillateur isotrope sont 
représentées sur la figure 4.10. La trajectoire d’en haut à gauche 


correspond aux valeurs initiales x (0) = y (0) = 0, x (0) = y (0) 
Æ O (l'oscillateur ne possède à l’instant initial qu'une énergie 
potentielle). La trajectoire est une hypocycloïde. La trajectoire d'en 
haut à droite qui est une rose correspond aux conditions initiales 


z (0) = y (0) Æ 0, x (0) = y (0) = 0 (à l'instant initial l’oscillateur 
ne dispose que de l'énergie cinétique). Les deux trajectoires du bas 
qui ont la forme d'hypotrochoïdes allongées correspondent à des 
conditions initiales non nulles. La figure 4.11 représente également 


pour ce dernier cas l’hodographe de la vitesse y = Ÿ (x); pour 
l'obtenir sur l’écran il faut appliquer à l’entrée d.v. 1 (dév. verticale) 
de l'indicateur la sortie de l’intégrateur [7] et à l'entrée d.h. (dév. 
horizontale) la sortie de l’intégrateur [5]. 

La figure 4.12 représente le comportement d'un oscillateur 
spatial isotrope lorsque l’énergie se dissipe suivant la loi de Rayleigh 


R = + + Ye). La perte d'énergie modifie les deux premières 
équations du système (4.6): 


z— —oiz+ y az, 
2 cH * : (4.7) 
Y= —Oyy—- TT — AY, 
où «& est le coefficient de proportionnalité entre la vitesse et les 
projections Q, = — , Q, = — _ de la résistance généralisée 


Q. Les derniers termes des équations (4.7) ne sont pas câblés sur 
la figure 4.9. Ils le sont par l'organisation des circuits de réaction 
dans les intégrateurs [5] et [7] par l'intermédiaire de potentiomètres 
de coefficients de transfert &. La trajectoire de droite correspond 
à de petites valeurs de &, celle de gauche aux valeurs doubles. 

La figure 4.13 représente la trajectoire d’un oscillateur dissipant 
son énergie en présence d’une faible anisotropie (w, — w,) entrai- 
nant une certaine précession. 


*) L'origine des trajectoires de l’oscillateur est marquée d'un point gras 
sur tous les dessins. 
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$ 5. Charge dans des champs croisés 


1. Position du problème. Le mouvement d’une charge dans 
divers champs est étudié en physique. Nous allons examiner le 
mouvement d'une charge sous l’action de deux champs constants 
et homogènes : un champ électrique et un magnétique. On se limitera 


Fig. 4.13 


au cas non relativiste où la vitesse de la charge est de beaucoup 
inférieure à celle de la lumière. Du cours de la physique on sait 
qu'il est nécessaire pour cela que le champ électrique soit petit 
devant le champ magnétique. Nous allons simuler les trajectoires 
du mouvement qui sont des courbes gauches. 

2. Equations différentielles du mouvement. Prenons pour axe 
z la direction du champ magnétique H et pour plan yOz le plan tendu 
sur les vecteurs H et E (fig. 4.14). Les équations différentielles 
du mouvement s'écrivent 
ET ELLES Z = 2e (4.8) 

m mc m 


ce eH e ee 
T=——y, y= , 
où e est la charge, m sa masse, c la vitesse de la lumière. 
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3. Schéma de commutation. 
Le schéma simulant le système 
eH 9 Ey _ 
me “E£E, 
— 10 est représenté sur la figure 


(4.8) pour ox = 


4.15. À la quantité # corres- 


pond une tension électrique de 
20 volts. Tous les paramètres 
peuvent varier dans des limites 
déterminées : la tension du champ 
électrique Æ à l’aide du diviseur 
D; les projections £E, et E, 
à l'aide du potentiomètre place 


Fig. 4.14 


à la (28)-ième entrée de l’intégrateur [15]; la quantité OH à l’aide 
des potentiomètres disposés aux (18)-ième et (26)-ième entrées des in- 


tégrateurs [5] et [7]. 


4. Résultats de la simulation. La simulation de la trajectoire 
de la charge permet de mettre en évidence l'influence des paramètres 


Fig. 4.15 


E, H,e et des conditions initiales sur cette trajectoire. La projection 
de la trajectoire sur le plan de coordonnées yOx conduit aux trois 
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types de courbes de la figure 4.16 en fonction de la condition initiale 
zo (intégrateur ([5]) : la première trajectoire à partir du haut est une 
cycloïde ordinaire (x, = 0); la deuxième, une cycloïde allongée 


ÿ 
2 
Fig. 4.16 
{zo = —5 volts); la troisième, une cycloïde tronquée (to — 5 volts). 
Si pendant la simulation de la trajectoire on élimine le champ 
ÿ 
A [IN À A. 


Fig. 4.17 


électrique E (par exemple en coupant la tension 100 volts du divi- 
seur D) la trajectoire de la charge se transforme en l’ellipse repré- 
sentée tout en bas de la figure 4.16. 
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La figure 4.17 montre la trajectoire *) de la charge (pour To = (0) 
d'origine au point À (4,, O, À,) et ses projections sur les plans de 
coordonnées zOz et yOz. Au voisinage des points où la trajectoire 


_ 
LA 7 


TRS RIISS 
ANAAPX 


Fig. 4.18 


est tangente au plan xOz, les projections sur xOz sont difficilement 
discernables, c'est pourquoi elles ont été effectuées sur un plan 


z'Oz’ parallèle à xzOy et situé légèrement plus bas. 

La figure 4.18 représente un faisceau de trajectoires correspondant 
à diverses valeurs du champ électrique E. On voit que lorsque Æ 
croît, la période (ou base) de la cycloïde décroit. 

On n'a pas représenté sur le schéma de la figure 4.15 les deux 
sommateurs nécessaires au tracé des trajectoires en fonction des 
variables zx (f), y (t), z (t). 


*) Pour la représentation des trajectoires voir Annexe. 
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$ 6. Système bielle-manivelle 


1. Position du problème. La figure 4.19 représente un schéma 
cinématique d'un système bielle-manivelle à coulisseau. La mani- 
velle OA effectue une rotation dans le sens rétrograde. Le coulis- 
seau B est animé d’un mouvement de va-et-vient le long de l’axe Oz. 
La bielle AB effectue un mouvement plus compliqué. Chacun de ses 
points décrit une trajectoire différente de celle des autres. L'ensemble 
de toutes ces trajectoires s'appelle famille des courbes de la bielle. 


2. Déduction d’expressions mathématiques favorables à la simula- 
tion. Supposons que la longueur de la manivelle OA soit égale 
à l'unité. La manivelle tourne à une vitesse angulaire constante «. 
L'angle ® qu'elle fait avec l'axe Oz sera alors égal à ot (i.e. = wt) 
et les coordonnées du point À de la manivelle seront 

zA (t) = cos wt, yA (t) = sin wt. (4.9) 


»« 


La longueur des éléments du système étant constante, à tout 
instant on aura l'égalité 


(8 — za) + ya — L° = 0, (4.10) 


où x, est l’abscisse du point B, L la longueur de Ia bielle AB. 
Considérons sur la bielle AB un point mobile M et déplaçons-le 

de À vers B. Si la vitesse de M est suffisamment petite, sa trajectoire 
nous donnera une idée de la famille de courbes de la bielle. Soit L (t) 
la distance de B à M. Les valeurs positives de / (t) sont dirigées 
de B vers À. Introduisons la fonction non dimensionnelle u ({) — 
— l(t)/L qui permet d'exprimer aisément les coordonnées carté- 
siennes zx (£{) et y (t) du point M: 

y (t) = yau (t), 

ZT) = 28 — (28 — za) u (t). 
Les expressions (4.11) donnent la solution du problème si de (4.10) 
on tire x, et on le porte avec (4.9) dans (4.11). Or l'expression de (4.10) 
est implicite en z;. Pour expliciter zx; utilisons la méthode des 
fonctions implicites en vertu de laquelle x; se détermine comme la 
solution de l’équation différentielle 

TB = — [(ZB — za)? + ya — L?] sign(zp—Z4a), Zs (0) = 0, 

où p est un nombre positif grand. Pour les systèmes bielle-manivelle 


existants, zp — za >> 0, donc la trajectoire du point est donnée 
en définitive par les expressions 


ZA (t) = cos wo, yA (t) = sin wt, 
zB = —P (xs — za) + y4a — LI, zx8 (0) = 0, 
nu (#) = 1(GYL, 
T=Zp—(zR—za)u(t), y = ya (t). 


(4.11) 


$ 6] SYSTEME BIELLE-MANIVELLE 125 


3. Schéma de commutation. Il est représenté sur la figure 4.20. 
Il correspond au cas w© = 10, ZL = 2. Il est construit pour une fonc- 
tion linéaire pu (t) définie sur l'intervalle u € [—1, 1]. Les opéra- 
teurs [15], [5], [6] élaborent les variables z\ (t), yA (t). L'intégra- 
teur [9] délivre la fonction pu (t). L'intégrateur monté sur la base 
de l’40 [10] reproduit x, (t) en résolvant l'équation implicite (4.10). 


8-NT -z,(0)-2—| 
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Fig. 4.21 
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4. Résultats. La trajectoire du point # est représentée sur la 
figure 4.21. Elle est simulée pendant 20 secondes. La durée de 
simulation de la trajectoire ne dépend que de la vitesse du point M 
sur la bielle 4B. Cette vitesse peut être modifiée à l’aide du poten- 
tiomètre placé à l'entrée (32) de l’intégrateur [8]. Le point M peut 
être fixé en tout lieu de la bielle, i.e. on peut prendre u (t) = const. 
Pour cela, après avoir mis la machine en marche, il faut rompre la 


= 


liaison à l'entrée (32) de l'intégrateur [8]. 


$ 7. Hystérésis 


1. Position du problème. Pour fixer les idées nous allons exami- 
ner l’hystérésis magnétique qui est le type le plus connu. L’hysté- 
résis magnétique se manifeste dans les modifications irréversibles 
de l’aimantation de ferromagnétiques dans un cycle d’aimantation. 
L'aimantation des ferromagnétiques est en rétard sur le champ 
magnétique extérieur qui est une fonction du temps. On rencontre 
des phénomènes analogues en étudiant des ferro-électriques, des 
systèmes mécaniques (lien entre les tensions et la déformation du 
matériel), etc. Tous ces phénomènes sont simulés d’une façon ana- 
logue. On se propose de simuler une boucle d’hystérésis de ferro- 
magnétiques et d'étudier ces derniers en fonction du champ exté- 
rieur et de leurs propriétés internes. 

2. Modèle mathématique de l’inversion de l’aimantation. Notre 
modèle mathématique est basé sur la théorie de l’anisotropie d’'Akou- 
lov en vertu de laquelle un monocristal de ferromagnétique se com- 
porte dans un champ extérieur de sorte à minimiser son énergie 
libre totale w qui est la somme de l'énergie magnétique d’anisotropie 
Lan et de l'énergie w' par rapport au champ magnétique extérieur *), 
soit 

w (j, H) = wan (f) + w°(j, H), 
où j est le vecteur d’aimantation, H l'intensité du champ magnétique 


extérieur. Pour tout vecteur H fixe, on détermine le vecteur j à par- 
tir de la condition d'équilibre thermodynamique 


grad w (j, H) = (. (4.12) 


Pour tracer la courbe statique d’aimantation j = ÿ (H), on se 
donne une collection de valeurs de Æ et l’on calcule les j correspon- 
dants à l’aide de la condition (4.12). 

Pour simuler le mouvement d’un processus thermodynamique 
déséquilibré j (t) vers son état d'équilibre il importe de tenir compte 
de la nature de ce mouvement. 


*) On suppose que l'aimantation du ferromagnétique est définie exclusive- 
ment par la rotation du vecteur d’aimantation. 
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Nous adoptons l'hypothèse linéaire qui est valable dans un 
certain voisinage de l’état d'équilibre et qui suppose que j (f) tend 


» 


vers l’état d'équilibre à une vitesse proportionnelle à grad w, i.e. 
H = —vgradw(j, H), (4.13) 


où v est une constante positive déterminée par le type de ferroma- 
gnétique. 

A titre d'exemple concret considérons l’inversion d'aimantation 
d’un ferromagnétique pour lequel Akoulov a trouvé analytiquement 
et étudié la courbe statique d’aimantation. Ce ferromagnétique est 
un monocristal à symétrie cubique ayant la forme d’un fil fin. Sui- 
vant Akoulov, l'expression de l'énergie libre aux termes multipli- 
catifs et additifs constants près est 

w (j, H)=4(L—a)—jA (t), (4.14) 
où æ est une constante déterminée par le ferromagnétique. 

En vertu de (4.13) on a pour (4.14) 
dj ï ; 
= v(H(t)— + aj]. 

En vertu de (4.12) l'équilibre est réalisé pour des valeurs de j et des 
valeurs de À égales à S, telles que 
S — ÿ + aj = 0. (4.15) 

3. Schéma de commutation simulant une boucle d’hystérésis. 

La figure 4.22 représente un schéma de commutation permettant 
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d'étudier l'hystérésis d’un ferromagnétique pour diverses valeurs 
des paramètres lorsque l’inversion d'aimantation s'effectue sous 
l’action d’un champ de la forme 
H (t) = exp (—Bf) sin wf + y, 

où f, w, y sont des constantes arbitraires. Le schéma fonctionne 
en deux régimes selon le signe de la tension qui est appliquée à l’en- 
trée du circuit de commande préprogrammée. Le premier régime 
auquel correspond la position des contacts du relais indiquée sur 
la figure simule la boucle sous forme d'une fonction explicite ÿ (H). 
Le deuxième régime simule une courbe en forme de S correspondant 
à l’état d'équilibre, d’après l'équation (4.15). Les opérateurs [14], 
[6], [5] délivrent la fonction H (t) = exp (—B$t) sin wt pour w — 10 
et un BEÏO, 1]. L'aimantation constante y est donnée à l’aide du 
potentiomètre du sommateur [12] à l'entrée (48). Les paramètres a 
et v du ferromagnétique sont donnés respectivement aux entrées (18) 
et (38) des opérateurs [9] et [10]. 

4. Résultats de la simulation. Sur la figure 4.23 sont représen- 
tées, à droite, la boucle d’aimantation et la S-courbe du minimum 
de l’énergie libre du ferromagnétique pour &« — 0, sous l’action d’un 
champ extérieur sinusoïdal d'amplitude constante (B = 0); à gauche 
est représentée l’inversion d’ aimantation du même ferromagnétique 
sous l'action d'un champ magnétique sinusoïdal dont l'amplitude 
décroit exponentiellement (B = O). 

La figure 4.24 représente à gauche et à droite l'inversion d'ai- 
mantation d'un ferromagnétique dont & = 0 sous l’action des champs 
sinusoïdaux extérieurs, d'amplitude constante et décroissante. 
Comme & = 0 la S-courbe possède une forme plus convexe qui cor- 
respond aux ferromagnetiques dont la boucle d’hystérésis est rectan- 
gulaire. La rectangularité de l'hystérésis est nettement marquée 
lorsque l’inversion d’aimantation est effectuée dans un champ d’am- 
plitude décroissante. Dans ce cas l’aimantation du ferromagnétique 
est caractérisée par l'existence de deux états stables dont l'un marque 
la fin de l’inversion d’aimantation. 

La boucle conserve sa forme pour diverses valeurs de v et w 


pourvu que le rapport _ = const. Si ce rapport diminue la courbe 
prend la forme d'une ellipse caractéristique de l'hystérésis linéaire. 


$ 8. Représentation d’une courbe gauche 


1. Position du problème. Etant donnée une courbe gauche sous 
la forme paramétrique 


Z = & exp (—a,t) cos wi, 
Ze = Aa EXP (—Wot) COS Wot, (4.16) 
Zs = 4, exp (—@;t) sin ouf, 


9—0112 
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on demande de la visualiser pour diverses valeurs des paramètres 
Œ; Œa; (OI Do. 

2. Déduction des équations différentielles. Les équations différen- 
tielles déterminantes pour (4.16) se trouvent sans peine par une 
dérivation successive de (4.16) si l’on introduit la fonction auxiliaire 


Za —= Go EXP (—at) Sin Wat. 
En dérivant successivement et en comparant les expressions des 
dérivées on obtient 
Z,= —QiZi — WyZs, Zi (0) = a4, 
Z2= — 222 — WzZ, 22 (0) — 2, 
Zs= —Gi23—++ WiZ1, Z3(0) = 0, 
24 = —Q2y + WoZr, 2, (0) = 0. 


Pour représenter la courbe gauche (4.16) à l'aide des équations 
déterminantes (4.17) il est nécessaire (cf. Annexe) de câbler les 


FL DS" 
z,le] (22) : 
(15) 
: ! 
z5l2] (46) 


(4.17) 


variables Ÿ,, Ÿ. qui projettent le système de coordonnées Oz,z:23 
sur le plan Ÿ,O0Y, confondu avec 2,02, : 


Y, = 2 — 0,523, Yo = 22 — 0,523. (4.18) 
3. Schéma de commutation. Il est représenté sur la figure 4.25. 


I] simule les équations différentielles (4.17) et les expressions (4.18). 
Il comporte un arrêt conditionnel au changement de signe de la 
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variable d'entrée du circuit de commande préprogrammée. Dans 
la position des contacts du relais représentée sur le schéma, le cir- 
cuit délivre la variable z, pour w, = 0 et &. quelconque. Lorsque 
la variable d'entrée change de signe, le circuit de commande prépro- 
grammée reproduit le cas où w: — 0. Sur le plan « électronique », 
les sorties des opérateurs [5] et [6] sont inhabituelles en ce sens 
qu’elles ont été faites plus puissantes par adjonction de résistances 
anodiques. 

4. Résultats. On rencontre les courbes gauches du type (4.16) 
lorsqu'on étudie les faibles oscillations de systèmes oscillants à plu- 
sieurs degrés de liberté. 

Trois cas facilement simulables à l’aide du schéma de commuta- 
tion de la figure 4.25 présentent un intérêt mathématique selon les 
valeurs numériques des paramètres w&, et 2. 

Premier cas. On suppose que w: — 0. L’équation (4.16) 
s'écrit 

Z, = ajexp (—at) cos wi, Ze = Go EXP (—Kat), 
Zs = Gi EXP (—œut) sinfw;f. (3.19) 


On détermine sans peine la nature de la surface (4.19), en élevant 
la première et la troisième expression au carré, en les ajoutant et en 
effectuant quelques transformations : 


On reconnaît dans (4.20) l’expression canonique d’un cône circulaire 
droit 


a #50 4.21 
TNT HET ns ee 


La surface (4.20) sera donc appelée cône circulaire généralisé. 
L'équation (4.20) se transforme en (4.21) pour a, = &.. Fait curieux, 
lorsque @&,/&+ = 0,5, l'équation décrit un paraboloïde de révolution. 
La surface (4.20) possède de nombreuses propriétés d’une conique: 

— c'est une surface de révolution autour de l’axe Oz: 

— elle possède deux nappes correspondant aux valeurs positives 
et négatives de 2. ; 

— ces deux nappes se touchent à l'origine des coordonnées qui 
est un point asymptotique de la courbe (4.16). 

En modifiant la valeur du rapport &./@&. on simule les courbes (4.16) 
appartenant à diverses surfaces de révolution, et partant on obtient 
la représentation graphique de la surface. La figure 4.26 représente 
trois trajectoires pour les rapports &/&. << 1, &./a: > 1 et,æy/ae — 1 
(de haut en bas). On change les coefficients &, et &. en éliminant 
et en rétablissant les liaisons des intégrateurs [5], [6] et [7] respecti- 
vement aux entrées (19), (23) et (27), (28). Pour tracer une courbe 


ge 
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Fig. 4.26 


de la deuxième nappe il faut changer le signe de la condition initiale 
de l’intégrateur [7]. Les courbes de la figure 4.26 donnent une image 
assez expressive de la surface et notamment de sa courbure. D'une 
façon générale, les lignes de courbure d’une surface de révolution 
sont des méridiens et des parallèles. Un centre de courbure est situé 
sur l’axe de révolution, l'autre sur la développée du méridien. 
La normale au méridien est dirigée vers l’intérieur de la surface si 
a/as << 1 et vers l'extérieur si æ&,/&. > 1. Le rayon de courbure 
du méridien devient infini dans le cas d'un cône droit. 
Deuxième cas. On suppose æ& — 0, ce qui se traduit, 
comme le montre la figure 4.25, par une élimination de l’entrée (28) 
de l’intégrateur [7] et la mise en position basse (représentée sur le 
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Fig. 4.27 Fig. 4.28 


dessin) du contact 4PO du relais opérationnel. La machine réalise 
alors la courbe 


Z, = 4, exp (—ayt) cos mil, 29 — Aa COS Wof, 
Z3 = y eXp (—@t) sin wf, 


qui admet l'interprétation physique suivante. On possède un oscil- 
lateur spatial effectuant dans le plan z,0z, des oscillations harmoni- 
ques amorties et simultanément des oscillations harmoniques non 
amorties le long de l'axe Oz.. La courbe gauche correspondante est 
tracée sur la figure 4.27. Elle est simulée par la machine à partir 
du point À (a, —a;, 0) dans le sens de la flèche. La projection 
de cette courbe sur le plan 2,07, est la spirale logarithmique repré- 
sentée dans la partie inférieure de la figure. La courbe dégénère 
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en une droite confondue avec l’axe Oz, lorsque les oscillations s'amor- 
tissent dans le plan z,0z,; pour de grandes valeurs de &. 
Troisième cas. La figure 4.28 représente la courbe qui 
correspond aux équations (4.17). On remarque que la trajectoire 
est issue du point À (—a,;, —a2, 0) et tend asymptotiquement vers 
l'origine des coordonnées. Pour le reste le processus se poursuit 
comme dans le cas précédent. Le coefficient de transfert @&, est 
donné par une réaction branchée à l'entrée (28) de l’intégrateur [7]. 
Cette réaction est représentée en pointillé sur la figure 4.25. 


$ 9. Courbe sur une quadrique 


1. Position du problème. On demande de tracer une courbe 
gauche située sur une quadrique: ellipsoïde, hyperboloïde (à une 
et deux nappes), cône, paraboloïde elliptique. 

2. Déduction d'expressions mathématiques favorables à lasimulation. 
A signaler au départ que ce problème est indéterminé, car par un 
point d’une quadrique il passe une infinité de courbes gauches. 
Nous essayerons de lever cette indétermination de manière à ce que, 
premièrement, la courbe tracée puisse donner une idée suggestive 
de la quadrique qui la contient, et, deuxièmement, pour que le 
schèma de commutation soit facilement réalisable. Les formes cano- 
niques des quadriques citées précédemment possèdent de nombreux 
points communs ce qui permet d'étudier tous les cas avec un seul 
schéma de commutation. Une courbe gauche est représentée para- 
métriquement par trois fonctions: 2, ({f), z.(t), z4(t), dont deux 
sont arbitraires et la troisième déterminée à l’aide de l'équation 
de la quadrique. 

Les formes canoniques des quadriques mentionnées s'écrivent : 


-2 2 2 
l'ellipsoïde: ++, 
1 2 


J 
se 2 ae 
l'hyperboloïde à une nappe: nt ae 
l'hyperboloïde à deux nappes: LR A — 1, (4.22) 
a; a3 a; 


le cône: =L— 2 + 5 —(), 


le paraboloïde elliptique: —2,+<#—0. 


Les quadriques (4.22) sont toutes des surfaces de révolution autour 
de l'axe Oz,. On se donnera tout naturellement la projection de la 
courbe gauche à simuler sur le plan z,0z, sous forme d’une spirale 
enroulante dont la représentation paramétrique est 


Z, = a exp (—at) cos wt, z3 — a exp (—at) sin wl. (4.23) 
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En d’autres termes, la courbe gauche à simuler est considérée comme 
l'intersection de deux surfaces: la surface (4.22) et une surface 
cylindrique dont la directrice (4.23) est une spirale logarithmique 
et la génératrice parallèle à Oz,. Dans la suite, pour simplifier les 
schémas on se limitera au cas a, = a, — a. Plus w est prépondé- 
rant devant «, plus la spirale logarithmique remplira le disque 


2, + z, = a. On obtient le tracé de la courbe à l’aide des formules 
Y, = 2 — 0,5 23, Ya —= 22 — 0,5 z, (cf. Annexe), où 2,, z et 23 
proviennent de (4.23) et (4.22). 

3. Schéma de commutation. Il est représenté sur la figure 4.29. 
Les opérateurs [5], [6], [1], [2] délivrent les fonctions z, (t) et z, (t) 
d’après leurs équations différentielles déterminantes. Les somma- 
teurs [3] et [4] élaborent les variables Y, et Y.. Les autres opérateurs 


Fig. 4.29 


sont utilisés à la résolution des équations (4.22) par rapport à 2, (t). 
A cet effet il faut cäbler les deux fonctions supplémentaires: v -= 
— exp (—2at), u = (+) a° (+) exp (—2at). 

4. Tracé de la courbe gauche. Passons à l’étude des diverses 
surfaces. 

1) Ellipsoïde. La courbe est tracée à l’aide du schéma 
de la figure 4.29. La variable z, (ft) est élaborée à l’aide de l’expres- 
sion 


Zo = + V'a?— exp (—2at). 


Le coefficient de transfert à l'entrée (40) du sommateur [14] doit 
être légèrement supérieur à l'unité (de l'ordre de 1,1). La courbe 
située sur l’ellipsoïde est tracée sur la figure 4.30 pour les valeurs 
positives et négatives de z.. Dans les deux cas elle est issue de À. 
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Fig. 4.30 


2) Hyperboloïde à une nappe. La variable z, (t) 
est réalisée à l'aide de la formule 


Z = + Æ V —a2+exp (—2at). 


La courbe est tracée jusqu’à l'instant où l'inégalité exp (—2at) > a° 
cesse d’être valable. Ceci étant, il faudra réduire jusqu’à 0,1 le 
coefficient de transfert du sommateur [14] à l’entrée (40); changer 
le signe de la tension à cette entrée ; changer le signe de la condition 
initiale de l'intégrateur [15]. 

La courbe est tracée deux fois (cf. fig. 4.31) : à partir du point À, 
pour les z, positifs et à partir de À, pour les z, négatifs. 

Au-dessus et au-dessous de l’hyperboloïde on constate la présence 
d’anneaux plans entièrement remplis par une spirale logarithmique 
enroulante. Ces anneaux se forment lorsque exp (—2at) << a*: le 
radicande de l'expression de z, devient négatif et z, ne possède plus 
de valeurs réelles correspondant à z, et z.:. Pour éviter l’apparition de 
ces anneaux il aurait fallu prévoir une interruption du tracé de la 
courbe par un circuit de commande préprogrammée à la condition 
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NI 


1} 


—a2 +exp (—2at)= 0. On adopte cette solution en pareils cas. 
Si nous avons refusé de le faire, c’est pour montrer comment se 
comporte la machine devant des racines imaginaires. La machine 
a donc « rejeté » la valeur imaginaire de la plage +100 volts et 
a continué à délivrer les variables z, et z; comme d'habitude. En 
réalité le changement de signe de la variable u (t) à l'entrée (37) 
de l’amplificateur [10] chargé de l'extraction de racine transforme 
la réaction négative en positive. 

3) Hyperboloïde à deux nappes. Le tracé dela 
courbe sur un hyperboloïde à deux nappes s'obtient par le câblage de 


2= + Var+exp(—2ut) 
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Fig. 4.33 Fig. 4.34 


qui nécessite le changement du signe de la constante (—100 volts) 
à l'entrée (40) du sommateur [14]. La figure 4.32 représente une 
courbe d'un hyperboloïde à deux nappes. 

1 est intéressant de tracer simultanément des courbes situées 
sur un hyperboloïde à une nappe et un hyperboloïde à deux nappes. 
Ce cas est représenté sur la figure 4.33. Le circuit de commande pré- 
programmée évite l’apparition de valeurs imaginaires de z, à la 
condition exp (—2œt) = a*. 

4) Cône. On obtient un cône en câblant la variable 


Za — + (a./a) exp (—at). 


A cet effet il faut annuler le coefficient de transfert à l’entrée (40) 
du sommateur [14] (en opérant une déconnexion à cette entrée par 
exemple). La courbe est représentée sur la figure 4.34. Les sommets 
des cônes auraient dû coïncider à l’origine des coordonnées, mais 
nous avons évité cela en utilisant un traceur à deux coordonnées. 
Pour tracer la courbe il faut augmenter sensiblement le temps de 
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résolution du problème, afin que le traceur ait le temps de « suivre » 
la courbe. Les fonctions 2, (t) et z; (t) varient très lentement au 
voisinage de l'origine des coordonnées. Le temps de reproduction 


Fig. 4.35 


de la courbe par la machine risquait d’être supérieur au temps ad- 
missible d'intégration. 

5) Paraboloïde elliptique. La variable z, ({) est 
délivrée par l’intégrateur [15] sous la forme 2, ({) = a-* exp (—2at). 
La courbe correspondante est représentée sur la figure 4.35. 


$ 10. Représentation graphique d’une surface de révolution 


1. Position du problème. Une surface de révolution est donnée 
par les équations 


z = (a + b cos v) sin u, 

y = (a + b cos v) cos u, (4.24) 

Z = C Sin v, 
où u et v sont les coordonnées géographiques : u la longitude, v la 
latitude. Les équations (4.24) définissent une surface de révolution 
(un tore) engendrée par la rotation de l'ellipse x = a + b cos v, 
z = c sin v autour de l'axe Oz. Lorsque a = U, le tore se transforme 


en un ellipsoïde qui à son tour se transforme en une sphère pour 
b = c. 
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On se propose de représenter graphiquement la surface (4.24). 

2. Déduction d’expressions mathématiques favorables à la simula- 
tion. Les lignes u—const et v— const du tore sont respectivement 
appelées méridiens et parallèles. La représentation graphique des 
méridiens et des parallèles donne une vue d'ensemble de la surface 
qui les porte. Mais le câblage de (4.24) conduit à des schémas volumi- 
neux et difficilement réalisables, c’est pourquoi il faut adjoindre 
aux équations (4.24) une relation reliant u et v. On obtient cette 
relation par les raisonnements suivants. À tout couple de nombres 
(u, v) mod 2x correspond un point et un seul sur la surface du tore. 
Soit une carte plane du tore rapportée à un réseau rectangulaire 
u = const, v = const. À tout point du tore est associé un point 
et un seul de la carte. A la droite 


U — Will, U —= ol, 


où w,, &w, Sont des constantes et { un paramètre, est associée sur la 
surface une courbe coupant les méridiens (respectivement les parallè- 
les) sous un certain angle constant. Les courbes qui coupent sous 
un angle constant les méridiens (les parallèles) de la surface qui les 
porte sont appelées lorodromies. Soient deux loxodromies orthogona- 
les, la deuxième ayant pour équation 


U —= — Oo, UV = Wjé. 


Si w/ws > 1 la première loxodromie donne une idée suggestive 
de la disposition des méridiens, la seconde, de la disposition des 
parallèles. Le câblage de la relation linéaire entre u et v conjointe- 
ment avec (4.24) simplifie beaucoup le tracé des surfaces de révo- 
lution. 

3. Schéma de commutation. Il est représenté sur la figure 4.36. 
Les opérateurs [7], [8], [2] élaborent les fonctions sin w;t, cos w;f; 
les opérateurs [16], [15], [13], les fonctions sin w.t, cos w.t d’après 
leurs équations différentielles déterminantes pour les valeurs numé- 
riques des paramètres @, = 10, w, — 0,3. Les opérateurs [1], [3], 
[4], [11], [12], [14] réalisent l'expression (4.24). Les paramètres 
a, b, c du tore peuvent varier dans des limites déterminées: a et b 
respectivement aux entrées (4) et (8) du sommateur [1]. c à l'entrée 
(44) du sommateur [14]. Les sommateurs [9] et [10] projettent le 
système de coordonnées zx, y, z sur le plan Y,0Y, de l'écran de l’in- 
dicateur qui coïncide avec le plan z0r. La projection s'obtient par 
génération des expressions YŸ, = z — 0,5 y, Y, = x — 0,5y (cf. 
Annexe). 

Le schéma comporte une commutation préprogrammée qui 
permet de tracer la première ou la deuxième loxodromie en fonction 
du signe de la variable d'entrée du circuit de commande préprogram- 
mée. 
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4. Tracé de la surface. La figure 4.37 représente le tracé du tore 
pour a = 1,0, b = c — 0,25 à l’aide de deux loxodromies perpendi- 
culaires. La figure 4.38 montre une sphère, la figure 4.39, les pro- 
jections des loxodromies de la sphère sur le plan zOy, la figure 4.40, 
le tore pour b = c = 1, a = 1 + e, où e est un petit nombre. 
Lorsque e = 0 le tore est engendré par la révolution d’un cercle 
autour d'une tangente figurée par l’axe Oz. La figure 4.41 représente 
séparément un tore tracé à l’aide d’une loxodromie méridien et d’une 
loxodromie parallèle. 


$ 11. Graphe de la solution 
de l’équation de la chaleur 


1. Position du problème. Nous allons étudier la distribution non 
stationnaire de la température dans une paroi épaisse. Ce problème 
a été décrit au chapitre 3, $ 5. On se propose de tracer par la méthode 
des sections la surface qui traduit la distribution de la température 
à l'intérieur de la paroi dans le temps et dans l’espace. Nous admet- 
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trons que la fonction v ({, y) est solution de l'équation de la cha- 
leur (3.11) qui est approchée par un système d'équations différen- 
tielles ordinaires (3.13) 

2. Schéma de commutation. Le schéma de la figure 4.42 donne 
le tracé de la solution v (t, y). La surface supérieure de la paroi 
(y = 0) a une température v (f, 0) = œ, (t), la surface inférieure 
(y = 1) une température vw (t, !) = œ,(t). La paroi est partagée en 
six couches par cinq lignes horizontales continues de coordonnées 
Yys Yos Ysr Vas Ys- Les variations de la température en fonction du 
temps, aux points d’intersection de ces droites avec l'axe des y, 
sont reproduites par les intégrateurs [15], [5], [6], [7], [8] et les 
inverseurs [1] et [2]. L’intégrateur [16] est utilisé pour élaborer la 
variable indépendante t et pour interrompre la résolution à l'aide 
d’un circuit de commande préprogrammée. Les sommateurs [9] 
et [10] élaborent les variables Y,; et Y., qui donnent la solution 
v (t, y;) en projection axonométrique par les sections y; (j = 1, 2, 
ns 9). En modifiant le coefficient de transfert du potentiomètre 
à l’ entrée (42) du sommateur [11] on donne la valeur de la coordonnée 
spatiale de la section, soit y;. Les valeurs numériques de ces coor- 
données représentées en variables machine sont inscrites à gauche 
de la paroi, de sorte qu'à y, correspond une tension de 15 volts, 
à y, une tension de 30 volts, etc. Les six intégrateurs de la machine 
MH-7M ne permettent pas de résoudre l'équation de la chaleur dans 
plus de cinq sections de la paroi. Donc le schéma est complété d’un 
interpolateur linéaire qui réalise l’interpolation au milieu de l’inter- 
valle des sections. L’interpolateur est constitué des sommateurs 
[3] et [4]. Les sections médianes de la paroi sont indiquées par des 
lignes horizontales en pointillé. La température est désignée 
Par V1, dans la section qui se trouve à mi-distance des sections 


j +iet j. 


3. Résultats de] la simulation. 


1) Refroidissement «libre». Une paroi chauffée 
uniformément jusqu'à une température constante (Co = Veg = 
= Vs — Uso — 100 volts), commence à se refroidir à partir de 
l'instant t — 0, puisque œ, (f) = œ2 (t) = 0 (les liaisons sont rom- 
pues aux entrées (32) de 181 8] et (48) “ [15]). La figure 4.43 repré- 
sente la surface v (ft, y) construite par le traceur de courbes section 
par section (y, — const). Si l’on dispose d’une vue d'ensemble de la 
surface v (it, y), il est aisé de compléter à la main le tracé des sections 
de la surface v (t, y) pour t — const. Une de ces sections est repré- 
sentée sur la figure. Les familles de sections de la surface y — const 
et t — const donnent une image assez expressive du refroidissement 
de la paroi. 


Fig. 4.43 


Fig. 4.44 
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Fig. 4.45 


Fig. 4.46 
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2) Refroidissement «forcé». La paroi chauffée 
préalablement jusqu'à une température constante (u,, = eo = 
= Ugo = Vao — Vo — 100 volts) commence à refroidir à partir de 
l'instant t — 0. La température de la surface supérieure est mainte- 
nue nulle: ®, (t) = 0 (rupture de la liaison à l'entrée (32) de l’inté- 
grateur {8]), la surface inférieure est maintenue à une température 
négative constante: @, (t) — const << O0 (application d'une tension 
électrique de — 100 volts à l'entrée (48) de l'intégrateur [15]. 
La figure 4.44 représente la surface correspondante. Sur le plan yOt 
on a joint par une ligne les points v (t, y;) — 0 qui correspondent aux 
instants où la température de la j-ième couche s’annule et devient 
négative. 

3) Réchauffement unilatéral. La paroi dont la 
température est nulle à l'instant initial (v,, = Ueo = V39 = Lao = 
— Uso — 0) commence à se réchauffer à partir de l'instant { = 0. 
En outre la température de la surface y = O (la supérieure) est 
constante, i.e. @ (t) = const => 0 (une tension de 100 volts est 
appliquée à l'entrée (32) de l'intégrateur [8]). La surface inférieure 
y = Î de la paroi est maintenue à une température constante q, (t) = 
— 0 (par déconnexion à l'entrée (48) de l’intégrateur [15]). Ce cas 
est illustré sur la figure 4.45. On a construit dix sections de la surfa- 
ce v (ft, y) à l'aide du traceur de courbes: cinq d'entre elles ont été 
obtenues par câblage de l'équation de la chaleur, les cinq autres 
(Yj+1, 7) par une interpolation linéaire afin d’'accentuer l'effet visuel. 
L'interpolation n’accroît pas de beaucoup le nombre d'opérateurs. 
Ici elle a demandé deux sommateurs en plus ([3] et [4)). 

4) Réchauffement accompagné d'un re- 
froidissement. La figure 4.46 représente la simulation de la 
distribution de la température dans l'épaisseur de la paroi lorsque 
la température p, (t) de la surface y = 0 est une constante positive, 
i.e. se (t) = const > 0 (une tension de 100 volts est appliquée 
à l'entrée (32) et la température de la surface y = L est une constante 
négative : pt = const << 0 (une tension de —100 volts est appliquée 
à l'entrée (48). La surface vw (f, y) présente une symétrie axiale. 
L'’axe de symétrie est la droite v (t, ys). Cette dernière est contenue 
dans le plan yOt et représente la température de la section médiane 
de la paroi (sortie de l'intégrateur [6]). 


$ 12. Graphe de la solution de l’équation de Poisson 


1. Position du problème. On donne une équation différentielle 
de Poisson 


Zu Tor = — 1 (4.25) 


dont la solution uw (y, x) doit être simulée sur le carré —1<< x € 1, 
—1 << y << 1 avec les conditions aux limites suivantes u (—1), x) — 
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— u({,z) = 0, u (y, —1) = u (y, 1) = 0. Les équations de Pois- 
son décrivent des processus stationnaires tels que la flexion d’une 
lame, d’une membrane, la torsion d’une tige, la distribution de la 
chaleur dans une lame, etc. On peut donner à notre problème, par 
exemple, l'interprétation physique suivante: soit un tube de sec- 
tion carrée à l'extrémité libre duquel est tendue une membrane 
élastique. À l’intérieur du tube on maintient une pression constante 
légèrement supérieure à la pression atmosphérique qui bombe la 
membrane. Si l’on fait un léger vide dans le tube la membrane 
s’incurvera. 

2. Déduction d’expressions mathématiques favorables à la simula- 
tion. Pour déduire des expressions mathématiques commodes à la 
représentation graphique de la solution de l'équation de Pois- 
son (4.25), il nous faut : effectuer une approximation aux différences 
finies de (4.25) par la méthode des droites, ce qui ramène la résolu- 
tion de l’équation différentielle aux dérivées partielles à celle d’un 


6) c) 


Fig. 4.47 


système d'équations différentielles ordinaires; trouver le polynôme 
d'interpolation qui permet de simuler la solution pour des valeurs 
arbitraires x = const, y = const; tracer la surface uw (y, x) par la 
méthode des sections. 

1) Approximation de l'équation (4.25) par 
la méthode des droites. La méthode des droites appli- 
quée à l'équation générale de Poisson 


ou 


+ os = f(æ, ÿ) (4.26) 


a été améliorée par M. Slobodianski moyennant une prise 
en considération des dérivées d’un ordre supérieur que de coutume 
pendant l’approximation aux différences finies. Si la solution 
u (y, x) de l'équation (4.26) est cherchée sur les droites y, = const 
= 1,2, ..., n — 1), l'approximation suivant Slobodianski con- 
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duit aux systèmes d'équations différentielles 
” ” ” 1 
Uj+1 (x) + Oui (x) +ui-; (x) + [uj41 (x) —2u; (x) +u;(z)]= 


= fa (x) + 10f,(z)+ fu (x) (=1, 2, ..., n—1), (4.27) 
dont l'erreur d'approximation est de l’ordre de h*, où h est la distan- 
ce entre les droites y;+, et y,. S'agissant de l'équation (4.25), si dans 
le domaine donné des valeurs de x et y (cf. fig. 4.47, a) on prend 
deux droites y, — 1/3 et y, — —1/3, on obtient en vertu de (4.27) 
le système d'équations différentielles ordinaires approximant (4.25) 
sous la forme du problème aux limites 


Us + 1Our + À [us — 2us] + 12=0, 
u,(—1)=u(1)—0, 

1Ouz + 5 +1 — Que +us] + 12= 0, 
Up (—1)=u (1) = 0. 


2), Interpolation. La symétrie du problème (4.25) 
(fig. 4.47, b) nous permet de conclure que u, (x) = u, (x) et au lieu 
du système (4.28) on a une seule équation 


(4.28) 


vu, u(—1)=u(1)=0, (4.29) 


où u = u, (x) = u. (x). Désormais il est aisé de composer la fonc- 
tion d’interpolation P (y, x). A cet effet exigeons que la fonction 
cherchée P (y, x) prenne les valeurs 


0, u (x), u (x), 0 
sur les droites 


y = —1, —1/3, 1/3, 1. 
La fonction d’interpolation se détermine sans peine. Elle vaut 
P (y, x) = —9/8 (y* — 1) u (x). (4.20) 


3) Représentation graphique de u(y,x) par 
la méthode des sections. Au lieu de la fonction u (y, x), 
on trace le graphe de son approximation (4.30) en tenant compte 
du problème aux limites (4.29). La figure 4.47, c représente le 
plan Y,OY, rapporté aux axes de coordonnées Oy, Oz, Op sur lequel 
sera projetée la surface P (y, x) à l'aide des expressions Y, = P — 
— 0,5 y, Ÿ, = x — 0,5 y (cf. Annexe) pour les sections y — const 
et à l’aide des expressions Ÿ, = P — 0,5 x, Ÿ, = y — 0,5 x pour les 
sections zx — const. 
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3. Schéma de commutation. Le schéma de commutation est 
représenté sur la figure 4.48. Le problème aux limites (4.29) est 
résolu à l’aide des intégrateurs [5], [6]. Le problème a été réduit 
àa un problème équivalent de Cauchy par la méthode des essais. 


z[7].yl9] 


F, 11 v. 
zÙ1.4 (9) (86) £ k 
A, 46 ob PT Q) 
(#8) 


Fig. 4.48 


La chaîne d'opérateurs [D], [9], [1], [2], [10], [3], [4] câble (4.30). 
L’intégrateur [7] et l'inverseur [13] élaborent la variable x, les 
sommateurs [11] et [12] réalisent le tracé de certaines sections de la 
surface P (y, x). 

4. Résultats de la simulation. La figure 4.49 représente la surface 
P (y, x) construite par ses sections y — const d’abord et x — const 
ensuite, à l’aide du schéma de la figure 4.48. La section y = const 
est donnée à l’aide du diviseur D. Aux valeurs positives de y est 
associée une tension négative à l'entrée du diviseur et inversement. 
Le problème étant symétrique, pour passer du tracé des sections 
y —= const à celui des sections z — const, il suffit de modifier la 
commutation des liaisons aux entrées des sommateurs [11] et [12]. 
Sur la figure 4.49, on a découpé le quart de la surface afin de donner 
une idée de l’incurvation de la membrane élastique. La coupure 
a été réalisée à l’aide du circuit de commande préprogrammée qui 
a interrompu la résolution à la condition x = 0. La figure 4.50 
représente la surface qui correspond à l’incurvation d’une membrane 
élastique tendue sur un tube de section carrée dans lequel on a fait 
un léger vide. 
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$ 13. Minimisation d’une fonctionnelle 


1. Position du problème. Etant donnée la fonctionnelle qua- 
dratique 


J'en 22) = + (@124 + 0274) + Go (4.31) 


OÙ dy 1, de Sont des constantes positives, on demande d'illustrer 
la minimisation de cette fonctionnelle. 

2. Déduction d’expressions mathématiques favorables à la simula- 
tion. Par illustration de la minimisation de la fonctionnelle, on 
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entend la simulation de la fonctionnelle J(x,,x.), de ses 
lignes de niveau J (x;, x.) — const et de ses trajectoires de plus 
rapide pente, qui sont décrites par les équations différentielles 


dx dr 
= pars, = —par,, (4.32) 


où p est une constante positive paramétrisant le processus. 
Une ligne de niveau est l'intersection de la surface J (x;, z.) avec 
le plan J — const, i.e. son équation est 
J (xz1, 2) = const. (4.33) 
En dérivant (4.33) par rapport à { on obtient 


0J dz; ôJ dr, 
07, dt 45 Ôts dt = 0: 


D'où l’on déduit le système d'équations différentielles déterminantes 
pour la réalisation des lignes de niveau 

dz dz 

mu = — DT), = OAiT4, (4.34) 


où «w est une constante non nulle définissant la vitesse de simula- 
tion des lignes de niveau. Il est aisé de‘remarquer que les systèmes 
d'équations différentielles (4.32) et (4.34) peuvent être représen- 
tés par un même système 


dr, _f —Paiz; (minimisation), 
dt = — 4,7, (lignes de niveau), 
dr, _f —Paz,; (minimisation), 
dt = &az; '(lignes de niveau). 


(4.35) 


La réalisation simultanée de (4.35) et (4.31) en tenant compte des 
formules de projection 


Y, = J —0,52z,;, Ya = Zzo — 0,5 


(cf. Annexe) résout le problème posé. 

3. Schéma de commutation. Il est représenté sur la figure 4.51. 
Il est composé des opérateurs [5], [6], [1] qui câblent le système (4.35). 
Les lignes de niveau sont réalisées dans la position des contacts du 
relais indiquée sur le schéma. La position inverse donne les trajec- 
toires de plus rapide pente. On peut modifier la position des contacts 
du relais opérationnel soit manuellement, en inversant le signe 
de la variable d'entrée du circuit de commande préprogrammée, 
soit automatiquement, en appliquant la sortie de l’intégrateur [8] 
à l'entrée du circuit de commande préprogrammée. La chaîne d'opé- 
rateurs [14], [7], {8] délivre la fonction 30 sin 104. La fonctionnelle 
J (x,, re) est minimisée par simulation de la trajectoire de descente 
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au cours d’une demi-période de la fonction 30 sin 10; les lignes de 
niveau sont réalisées pendant l’autre demi-période. Les lignes de 
niveau et les trajectoires situées sur la surface décrite par (4.31) 
sont représentées pour donner une idée suggestive de la minimisation 
de la fonctionnelle J (x,, x.). Les opérateurs [3], [4], [11], [12], [2] 
délivrent J (x,, ze), les sommateurs [9] et [10] les variables Ÿ, 
et Ÿ:. Si l'on rompt la liaison du sommateur [9] à l'entrée (35), 
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le processus de minimisation sera représenté en projection sur le 
plan x,Ozr.. Les paramètres numériques &,s, &, Ge, p, w du schéma 
de la figure 4.51 sont choisis tels que pa, = pa, = 0,3, wa, = was = 
—= 50, a, = 50. La quantité a, peut être modifiée à l’aide du poten- 
tiomètre placé à l'entrée (8) du sommateur [2]. 

4. Résultats. La solution exacte du problème (4.31) est évi- 
dente: le minimum de J'(x,,x2) est réalisé pour x, = x, = 0 et 
vaut as. Mais c’est bien moins le résultat numérique qui nous inté- 
resse que l'illustration du processus de minimisation. 

La figure 4.52 représente la fonctionnelle (4.31) à l’aide de la 
courbe 


z, = 90 exp (—at) cos 50, zx, — 50 exp (—at) sin 501, 


qui est située sur la surface J (x,, x.). Les variables x, et x, sont 
élaborées par les opérateurs [5], [1], [6]. Le coefficient & est défini 
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Fig. 4.53 


$ 14] REDRESSEMENT AVEC FILTRAGE 157 


à l’aide des potentiomètres placés aux entrées (18) et (22) des intégra- 
teurs [5] et [6] indiqués en pointillé sur la figure 4.51. 

La figure 4.53 représente les lignes de niveau et les trajectoires 
de plus rapide pente et leurs projections sur le plan x,Oz.. A signaler 
que si le schéma est correctement monté, le voltmètre branché 
à la sortie du sommateur [2] doit indiquer une tension électrique 
constante lors de la simulation des lignes de niveau. 


$ 14. Redressement avec filtrage 


1. Position du problème. Le redressement sert à obtenir du 
courant continu pour alimenter par exemple des appareils radio- 
électriques. On redresse la tension du courant alternatif. Le filtrage 


£() ec] ut) 
D— 


Fig. 4.54 


a pour but d’affaiblir la composante alternative de la tension re- 
dressée. 

On demande de construire le schéma simulant le fonctionnement 
du filtre d’un redresseur (fig. 4.54). A l’entrée du filtre est appliquée 


Fig. 4.55 


la tension redressée Æ (t) à une ou deux alternances obtenue à partir 
d'une tension sinusoïdale. Il existe dans le schéma du filtre un élé- 
ment réactif pour affaiblir la composante alternative : c'est un con- 
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Fig. 4.56 


densateur de capacité C. Le filtre fonctionne avec une résistance de 
charge Ron. 

2. Equation différentielle. L’équation différentielle décrivant 
le fonctionnement du filtre se détermine à l'aide des lois de Kirch- 
hoff. Elle s’écrit 


d'()=T[E-(1+)ue |. u(0)=0. (4.36) 


3. Schéma de commutation. Le schéma de commutation est 
représenté sur la figure 4.55. La chaîne des opérateurs [1], [5] et [6] 
élabore les variables sin t, —sin t{ nécessaires à la formation de la 
variable E (t). Si la tension redressée est à deux alternances, 
la variable E (t)"\ est formée à l'aide de l'opération de choix 
max {—sint, sint} câblée par l’inverseur [13] avec deux diodes 
D, et D., branchées à son entrée. Si la tension est à une alternance, 
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la variable Æ (t) est formée à l'aide de l'opération de choix 
max {— sin t, 0} ou max {0, sin t}; dans ce cas il faut débrancher 
l'une des deux diodes. La chaîne d'opérateurs [7] et [2] intègre 
l'équation (4.36). Les coefficients de transfert des potentiomètres 
placés aux entrées (28) et (32) correspondent à la plus grande valeur 
de la résistance de charge R.n. Le coefficient de transfert du poten- 
tiomètre placé à l’entrée (32) de l’intégrateur [7] vaut (4 + r/Ren)/rC. 

4. Résultats. La figure 4.56 représente les tensions d’entrée 
E (t) et de sortie u (t) du filtre pour divers E (t). Le choix d'une 


Source z(t) | Emetteur | ytt)l/+mztilsinot Récepteur 
d'information \fnfor-” (modulation) Canal de liaison AM (détection) 
ation 


Fig. 4.57 
tension Æ (t) déterminée s'effectue par une commutation adéquate 
des diodes D, et D, à l'entrée de l’inverseur [13]. Le schéma de la 
figure 4.55 permet d'estimer la qualité du filtrage et l'influence 
exercée sur elle par les paramètres du filtre r, C et de la charge Rex. 


$ 15. Modulation d'amplitude et détection 


1. Position du problème. La figure 4.57 représente schémati- 
quement un système de communication. Îl est composé de quatre 
parties. 

1) Source d'information. La source choisit une in- 
formation à transmettre dans un ensemble d'informations possibles. 
En radiodiffusion les informations sont des fonctions continues de 
basse fréquence. Désignons par x (f) l'information qui sera diffusée. 

2) Emetteur. Il transforme l'information zx (f) en un signal 
radio y(t) qui est transmis par un canal. Pour transmettre une 
information par émission d'ondes électromagnétiques, il faut, au 
cours de la transformation de l'information en signal, déplacer 
le spectre de fréquence de l'information dans le domaine des hautes 
fréquences radio. Cette transformation d’un signal basse fréquence 
en un signal haute fréquence constitue une modulation. Dans la 
modulation d'amplitude le signal y (t) est formé à partir de l’infor- 
mation x (t) conformément à l'expression 


y (t) = [1 + mx (t)] sin ot, (4.37) 
où w est une fréquence fixe, appelée fréquence porteuse, m un nombre 


positif appelé faux de modulation. Le nombre m est choisi inférieur 
à l'unité et tel que l'expression se trouvant entre crochets soit 
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positive. La modulation d'amplitude traduit le fait que l’'informa- 
tion xæ(t) agit uniquement sur l’amplitude du signal radio y (t) 
dans l'expression (4.37). 

3) Canal de liaison. Le canal de liaison est le milieu 
utilisé pour, la transmission du signal de l'émetteur au récepteur. 
En radiodiffusion, le canal de liaison est une bande déterminée de 
fréquences radio. Au cours de la transmission le signal peut se 
distordre ou être perturbé. Par souci de simplicité on supposera 
que le canal utilisé est parfait, c'est-à-dire que le signal n'est ni 
distordu ni perturbé. 

4) Récepteur. Il transforme le signal reçu y (t) en l’in- 
formation de départ. Le rétablissement de l'information n'est pas 
idéal, aussi désignera-t-on par Z (t) l'information délivrée par le 
récepteur pour la distinguer de l'information initiale x (t). Sur le 
plan mathématique, le récepteur réalise la fonction inverse de l’émet- 
teur. En radiotechnique l'opération inverse de la modulation s’ap- 
pelle démodulation ou détection. Donc la démodulation a pour objet 
un signal modulé en amplitude y (ft) et pour but la transformation 


de y (t) en information x (t). La démodulation présente de nombreux 
points communs avec le redressement (cf. $ 14). Mais on relève des 


distinctions fondamentales. En effet, si le schéma du redressement 
avec filtrage est destiné à éliminer toute composante alternative 
et à obtenir un courant continu, le schéma de démodulation doit 
dégager de l'amplitude variable du signal la basse fréquence qui 
correspond à l'information. 

La figure 4.58 représente le schéma d’un récepteur élémentaire 
comprenant un détecteur (redresseur) et un filtre. Le récepteur 


délivre la variable x (t). On se propose de construire un schéma qui 
permette de dégager à l’aide de la machine l'influence des paramètres 
du récepteur et de la fréquence porteuse w sur la perturbation de 
l'information x (t). 

2. Relations mathématiques favorables à la simulation. Suppo- 
sons que x ({) — 1 — cos t. La source d’informations est alors simu- 
lée par le schéma qui câble la fonction x (t) d’après son équation 
différentielle déterminante x” — 14 — x, x’ (0) — x (0) = 0. Suppo- 
sons par ailleurs que la fréquence porteuse w du signal radio est 


D, 


D g(t)[4) Ç > Le 
| POULE 

Da 2t | r ue 1 
@ : De D 


YU) 


nono 


Ê 
Fig. 4.60 


162 DIVERS PROBLÈMES DE SIMULATION [CH. 4 


égale à 20 fois la fréquence de l'information, i.e. z ({) — sin 20 ft, 
où z (t) désigne les oscillations de la fréquence porteuse, et que le 
générateur de la fréquence porteuse est simulé par le schéma câblant 
l'équation différentielle 


z” — —400 z, 2z(0) — 0, z° (0) = 20. 


En vertu de (4.37), la simulation se réduit à une multiplication. 
L'équation différentielle décrivant le travail du récepteur est 
connue du $ 14. Dans nos notations, elle s'écrit 


r'H=[lr01-(1+5)70]. (4.38) 


3. Schéma de commutation. Il est représenté sur la figure 4.59. 
La partie gauche correspond à la modulation d'amplitude, celle 
de droite à la démodulation. Les opérateurs [2], [7], [8] élaborent 
l'information x (t). Les opérateurs 
[1], [5], [6] simulent le générateur 
de la fréquence porteuse. Les 
opérateurs [9], [3], [4] réalisent 
la modulation d'amplitude. Le 
taux de modulation m est donné 
à l’aide du potentiomètre disposé 
à l'entrée (40) du sommateur 
[9]. Les inverseurs [12], [13] et 
les diodes D, et D, délivrent la 

Fig. 4.61 variable | y (t) | par câblage de 

la relation évidente | y (t) | = 

— max {y (t), —y (t)} à l'aide d’une opération de choix. Dans les 

termes du phénomène simulé, cette opération traduit le redresse- 

ment du signal radio et la multiplication par deux de sa fréquence. 

Les opérateurs [15] et [14] simulent l'équation différentielle (4.38). 

Les potentiomètres placés à l'entrée de l’intégrateur [15] permet- 

tent d'attribuer aux paramètres r, C et À, qui caractérisent le fonc- 

tionnement du récepteur, une valeur quelconque prise dans leurs 
intervalles respectifs de variation. 

4. Résultats de la simulation. La figure 4.60 représente quatre 
oscillogrammes tracés avec un traceur de courbes d’après le schéma 
de la figure 4.59. De haut en bas on a l’oscillogramme de l'informa- 
tion x (t), l’oscillogramme du signal radio y (t) module en amplitude, 
l'oscillogramme du signal radio redressé et de fréquence double, 
l'oscillogramme de l'information 7 (t) élaborée par le récepteur. 
L'information représente la composante variable exigée et la com- 
posante continue. L’oscillogramme montre que la composante 
variable de l'information x (t) répète avec un certain retard l’infor- 
mation initiale zx (t). 
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$ 16. Génération d’oscillations 


1. Position du problème. La figure 4.61 représente le schéma 
électrique d’un générateur à tubes. On demande de simuler l'établis- 
sement des oscillations dans le générateur. Le générateur de la 
figure 4.61 a ceci de particulier que le circuit oscillant est intercalé 
dans le circuit de grille. Dans la figure 4.61 i désigne le courant 
du circuit, à, la composante variable du courant anodique, uw la 
tension du circuit (aux bornes du condensateur), u, — u — E, la 
tension de la grille, £, la polarisation de grille. 

2. Déduction des équations différentielles décrivant le fonction- 
nement du générateur. Si l’on admet l'’inexistence des courants 
de grille en raison de la polarisation suffisante Æ,, la loi de Kirchhoff 


donne 


LE+rit++ | dt = M Le. (4.39) 


Au second membre, on a une f.é.m. induite dans le circuit oscillant 
par l’inductance de couplage. Si l’on tient compte du fait que la 
caractéristique tension-courant grille-plaque peut être représentée 
au voisinage du point de fonctionnement par une fonction de la ten- 


1 
sion de grille, soit i, — ® (u), et que u + {: dt, alors de (4.39) 
on déduit que 


d°u MS (u) 


nes tt 2 — 
et | + ou = 0, (4.40) 
OÙ Oo — = désigne la fréquence de résonance du circuit. Pour 


passer de (4.39) à (4.40), on s’est servi de la relation 
dix __ d® du du 


où S (u) — Te est la pente de la caractéristique grille-plaque. 


Lorsque l'oscillateur fonctionne en régime doux, la Se tes D 
grille-plaque est bien approchée par la fonction D (u) = Su — aus. 
Compte tenu de cette fonction, (4.40) s'écrit 


y  [3aM , (MS r\Vdu , » 
mt (5o- Tr) ] à + ou —0. 
Pour simplifier l'écriture, posons 


__ 3aaf MEL ie 
HTC » LC — L) 3aM 


L'équation s’écrit alors 


= (5-7 )/30. 


d? d a 
Te HU (u2— y) À + ou = 0. 


11% 


164 DIVERS PROBLÈMES DE SIMULATION [CH. 4 
A 


La substitution 
du us 
OUT TH (vu) 


nous conduit au système d'équations différentielles du premier ordre 
équivalent 


U" = OU —p (+ vu) , 


V'= — Opu. 


(4.41) 


Si pu = 0, le système (4.41) se transforme en un système d'équations 
différentielles linéaires décrivant les oscillations d’un circuit sans 
pertes. 

L'étude analytique du système (4.41) sur le plan de phase u, v 
ne permet d'établir que l'existence d’un cycle limite et d’un point 
singulier en l'origine des coordonnées. Si y > 0 ce point est toujours 
instable. En effet, le système linéaire de première approximation 
peut être mis sous la forme 


u" = GQ0 + Yu, Lv' = —ou 


pour u et v petits. Les racines caractéristiques sont toutes deux 
réelles et positives pour uy 2, et complexes à parties réelles 
positives pour puy << 2w,. Donc l’origine est un nœud instable pour 


/ 8-W “fl rs 


uy > 20, et un foyer instable pour puy << 2w,4. Il est impossible 
d'étudier analytiquement l'établissement de l'amplitude de l'oscil- 
lateur pour diverses valeurs des paramètres u, Ÿ, w, en raison des 
difficultés entraînées par la résolution du système (4.41). 

3. Schéma de commutation. Le schéma de commutation qui est 
représenté sur la figure 4.62 simule le système (4.41). Le sommateur 
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[9] délivre la somme —(£+ vu ) . Le coefficient de transfert du 


potentiomètre placé à l'entrée (34) ne vaut pas 1/3 mais 1,5 puisque 
on a tenu compte de l'échelle des multiplieurs élaborant u*. On fait 
varier le paramètre y de l’oscillateur à l’aide du potentiomètre 
placé à l'entrée (36) du sommateur {9]. La quantité u appliquée 
à l’entrée (18) de l’intégrateur [5] peut prendre des valeurs quelcon- 
ques comprises entre 0 et 11. Le schéma de commutation de la figure 
4.62 correspond à wo = 1. 

&. Résultats de la simulation. La forme des oscillations engen- 
drées varie en fonction de u. Ces variations sont traduites aussi bien 
par les oscillogrammes u (t) que par la trajectoire de phase. La figu- 
re 4.63 représente l'allure des oscillations pour trois valeurs distinc- 
tes de u. En bas on a la courbe de phase et l’oscillogramme des 
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Fig. 4.63 


oscillations pour u relativement petit. Le cycle limite est une ellipse. 
Les oscillations générées sont presque sinusoïdales. Le point initial 
u (0) = & (0) = 0 est un foyer instable. Au milieu on a la courbe 
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de phase et l’oscillogramme des oscillations générées u ({) pour de 
plus grandes valeurs de u. Les oscillations se distinguent nettement 
d'oscillations sinusoïdales même si la forme du cycle limite ressem- 
ble encore à une ellipse. La fréquence des oscillations a diminué, 
l'amplitude prend plus rapidement ses valeurs permanentes. En haut 
on a la courbe de phase et l'oscillogramme des oscillations pour pu 
relativement grand. Le cycle limite rappelle un parallélogramme ; 
les oscillations sont essentiellement non linéaires et rappellent celles 
d'un multivibrateur. Le point initial u (0) = v (0) = O0 est un 
nœud instable; l'amplitude prend plus rapidement ses valeurs 
permanentes. La fréquence des oscillations a pratiquement diminué 
de moitié par rapport au premier cas. 


$ 17. Excitation paramétrique (paramétron) 


1. Position du problème. On demande de construire un schéma 
simulant le fonctionnement du paramétron. Dans le paramétron 
les oscillations sont excitées paramétriquement. Sur le plan électri- 
que, le paramétron est un circuit résonnant à élément réactif non 
linéaire. Le circuit oscille à la demi-fré- 
quence des oscillations génératrices. Le 
paramétron est caractérisé par l'existence 
d'états stables. Cette particularité a sug- 
y géré à son inventeur, le savant japonais 
Hoto, de l'utiliser comme élément de cal- 
culatrices numériques pour représenter 
l'unité binaire de l'information en choi- 
sissant l’un des deux états stables dis- 

Fig. 4.64 tants l’un de l'autre de x radians. La fi- 

gure 4.64 représente le circuit résonnant 

à excitation paramétrique comprenant une inductance variable L ({), 

un condensateur C et une résistance RÀ ; u désigne la tension aux bor- 
nes du circuit, à le courant traversant l’inductance. 

2. Equationsdifférentielles décrivant le fonctionnement du paramé- 
tron. La loi de Kirchhoff donne 


R 


++ 0. (4.42) 


Si l’on tient compte de ce que la tension du circuit u — Ë (Li) 


et l'inductance varie suivant la loi L(t) = L, (1 + 2y sin 2wt), 
où 2w est la fréquence d'excitation du paramétron, y le coefficient 
de modulation, L, la composante constante de l’inductance, et si 
l'on introduit la notation 


L . ; À 1 1 2 
J=ri-=({+27sin2ot)i, Ô=—R; TL, = %(1+a), 
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alors on déduit de (4.42) l'équation différentielle 


+ Bo À + w7 (1 + a — 2y sin 2wt) J = 0. (4.43) 


On admet que a & 1 et y< 1 donc (1 + a)/(1 + 2y sin 2wt) 
a été remplacé par 1 + & — 2y sin 2wt. 

Dans l'équation (4.43) Ô représente l'amortissement, & le désac- 
cord du circuit. 

Si & et Ô sont constants, l'équation (4.43) est une équation dif- 
férentielle linéaire appelée équation de Mathieu. Elle décrit 
les processus paramétriques d’excitation d'’oscillations dans les 
systèmes linéaires. Les processus se déroulant dans le paramétron 


Fig. 4.65 


sont essentiellement non linéaires. Dans le circuit de la figure 4.64, 
la non-linéarité provient de l’utilisation d’un noyau ferromagnétique 
Aans la bobine d'’inductance. Avec l'accroissement de l'amplitude 
du courant la dépendance non linéaire de l'induction B par rapport 
à l'intensité Æ du champ magnétique entraîne un désaccord supplé- 
mentaire du circuit résonnant et un accroissement simultané des 
pertes par hystérésis. Donc, l’amortissement Ô et le désaccord «& 
seront en général fonction de l'amplitude J*. Si la non-linéarité 
est prise en considération sous forme d’un désaccord supplémentaire 
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BJ* du circuit, dre s’écrit en définitive 
A Le + 80 © + @02(1+a+fJ2—2ysin2ot)J=0. (4.44) 


La décomposition 
J= J, sin ot + J, cos wt, (4.45) 


où J, est l'amplitude de la composante sinusoïdale et J. l'amplitude 
de la composante cosinusoïdale, et la substitution de (4.45) dans 
(4.44) donnent le système 


Ji + ye—[a+$ 80478 ]Je, 


2 pr. _ à  _ (4.46) 
Det e—v/s+[a+ 28342 7. 
On admet que «a, y, ô € 1 et que w > 1. 

3. Schéma de commutation. Le schéma de commutation qui 
est représenté sur la figure 4.65 correspond à B < 0, y > 6, t — 
— mut/2, t étant le temps sans dimension et m un nombre positif 

jouant le rôle d'échelle. On peut attribuer 

j: à &, P, y, Ô une valeur quelconque prise 

dans leurs intervalles de variation à l’aide 
des potentiomètres. 

4. Résultats. Les processus se déroulant 
dans le paramétron se représentent aisément 
sur le plan de phase J.OJ. (fig. 4.66). A 
tout point de ce plan est associé dans un 
système de coordonnées polaires un rayon 
vecteur À dont le module est égal à l'am- 


plitude des oscillations (R = VJE+ JE) et 
Fig. 4.66 l'angle polaire o à la phase des oscillations. 
Le cas le plus compliqué du comportement 
du paramétron est représenté sur la figure 4.67. La courbe de phase pré- 
sente cinq points singuliers : deux points selle S, et S.et trois foyers sta- 
bles alignés. Le foyer du centre correspond à l'origine des coordonnées. 
Aux deux autres foyers sont associées deux oscillations fixes du 
paramétron déphasées de x radians. Le paramétron arrive à l’un 
des trois états stables quelle que soit la position initiale. De tels 
paramétrons sont dits à trois états stables. Les valeurs numériques 
des coefficients de transfert des opérateurs de la figure 4.65 corres- 
pondent à un paramétron à trois états stables. Si l’on annule & 
(par rupture de la liaison à l'entrée (44) du sommateur [14]) le para- 
métron devient bistable. Le point singulier J, = J, = 0 qui repré- 
sente l’origine du plan de phase devient un nœud instable. 
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Fig. 4.67 


$ 18. Propagation d’une épidémie 


1. Position du problème. Construire un schéma simulant la 
propagation d'une maladie contagieuse qui s’est déclarée dans une 
localité retirée de H habitants. Au moment où l'épidémie a été 
découverte, y, étaient malades, donc agents infectieux, x, étaient 
encore sains mais potentiellement malades, les autres, soit z, = 
— H — x, — y,, relevaient de l'infection et étaient immunisés. 
Le caractère de l'épidémie et les conditions de vie de cette localité 
sont tels que chaque malade contamine quotidiennement en moyen- 
ne B personnes sur mille, chaque malade est immunisé au bout de y 
jours de convalescence. 

2. Equations différentielles. Soient x les personnes saines mais 
malades potentiels, z les personnes saines et immunisées, y les 
malades (agents infectieux). La propagation de l'épidémie est 
décrite par le système d'équations différentielles 


z= —10"Pzy,  z(0)= 20, 
y= 10 Pzy—v "y, y(0)= Yo, (4.47) 
z=% y, z (0) = 20. 


12—0112 
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On admet que l'infection se propage par contact, c'est-à-dire que 
la vitesse de contamination des personnes saines x est proportion- 
nelle à la fréquence de leurs rencontres avec les malades y. 

3. Schéma de commutation. La figure 4.68 représente le schéma 
de commutation simulant le système (4.47). Les valeurs numériques 


Fig. 4.69 


des coefficients de transfert des opérateurs correspondent aux cas H = 
— 1000 personnes, B — 1 personne sur mille, y — 10 jours. Une 
seconde machine représente une journée. On peut faire varier H 
dans un intervalle assez large à l'aide des potentiomètres placés 


aux entrées (28), (8), (24) et (20) des intégrateurs [7], [6], [5]. 
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4. Résultats de la simulation. La figure 4.69 représente les 
oscillogrammes zx (t), y (t), z (t) obtenus à l'aide du montage de la 
figure 4.68 sous l'hypothèse que x, — 900 personnes, y, = 90 per- 
sonnes, Z9 = 10 personnes. La variable zx ({) est monotone décrois- 
sante et tend asymptotiquement vers zéro; la variable y (f) présen- 
te un extremum à partir duquel elle tend asymptotiquement vers 
zéro; la variable z (t) est monotone croissante et tend asymptotique- 
ment vers H. 


$ 19. Interaction de populations 


1. Position du problème. Construire un schéma simulant l’évo- 
lution quantitative de deux populations animales coexistantes. 
A l'instant { = 0 pris pour initial, le nombre de sujets de la pre- 
mière population est de z,,, celui de la deuxième de zzo. 

2. Equations différentielles. Pour fixer les idées étudions le 
modèle correspondant au cas où les rapports entre les populations 
animales sont du type prédateur — proie, i.e. les animaux d'une 
espèce, par exemple des lièvres, constituent la pitance des animaux 
de l’autre espèce, par exemple des lynx. Pour simplifier, tous les 
facteurs supplémentaires influant sur la modification du nombre 
des animaux seront supposés constants. Ceci posé, l’évolution du 
nombre de sujets est décrite par le système d'équations différentielles 

Ti = GT —GyoTiTo y, Ayo > 0, 
: Ti diotiTo ip, Au > (4.48) 


Lo = Anal To — solo, doi; UP) > O, 


qui suppose qu'en l'absence de lynx (re — 0) l'accroissement des 
lièvres est proportionnel à leur nombre x,. La présence de lynx réduit 
la prolifération des lièvres. Cette réduction est proportionnelle 
au nombre de rencontres des lièvres et des lynx. Au contraire, le 
nombre de lynx croît avec celui des rencontres; la disparition des 
lièvres entraîne celle des lynx et ce d'autant plus vite que ces derniers 
sont nombreux. 

Sur le plan de phase r,Ozr, le système (4.48) présente un point 
singulier de coordonnées r;6q = G92/@o1r Toéq — 11/09. Si le nombre 
initial d'animaux est tel que 2z10 = Ziégr Teo — Toéqg» il restera 
constant. Analysons l’évolution quantitative des animaux en fonc- 
tion du temps tout d'abord dans un petit voisinage du point d’équi- 
libre. À cet effet, linéarisons le système (4.48). L'introduction des 
variables u, = x; — Goo/aoy et Us — Ze — aj1/a3, nous donne le 
système linéarisé 
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ou encore 


Uy—= —Gjijlqolis 
d'où il suit que 
u, = À cos of, u, = B sin ot, 


où À et B sont des constantes et w — V a@.. Maintenant il est 
aisé d'établir qu’au voisinage du point singulier la trajectoire 
z (t), y (t) est une ellipse, puisque 


Les possibilités de l'étude analytique du comportement de x (ft) 
et y (t) s'arrêtent là. Nous allons utiliser la machine pour repré- 
senter les courbes de phase du système (4.48) sur le plan de phase. 

3. Schéma de commutation. Le schéma de commutation est 
représenté sur la figure 4.70. Dans la position des contacts du relais 
indiquée sur la figure, le schéma simule le système (4.48) pour les 
valeurs numériques suivantes des paramètres : 4, = Goo = À, Aya = 
= Ge, = 2. On peut attribuer aux paramètres &@,;, Gjo, Goys Goo Une 
valeur quelconque prise dans leurs intervalles de variation à l’aide 
des potentiomètres disposés aux entrées (20), (18), (22) et (24) des 
intégrateurs [5] et [6]. Les variables zx, (ft) et x, (t) sont appliquées 
aux entrées de l'indicateur électronique pour une visualisation de 
la courbe de phase et des oscillogrammes. 

&. Résultats. La figure 4.71 représente les oscillogrammes 
x (t), Ze (e) et la courbe de phase x, = 1 (x,) dessinée par un traceur 
de courbes à l'aide du schéma de la figure 4.70. La courbe de phase 
est fermée. Les oscillogrammes zx, (t) et x, (t) sont des fonctions 
périodiques du temps. Leur examen permet de tirer des conclusions 
importantes: les lynx ne peuvent pas dévorer tous les lièvres; 
l'équilibre œcologique entre lynx et lièvres revêt un caractère 
dynamique; le nombre moyen de lynx et de lièvres est invariable 
POUT js Aer Œeys Ge Constants; les nombres initiaux x,, de lynx et 
z,0 de lièvres influent sur l’amplitude des oscillations zx, ({) et x, (t) 
autour des valeurs moyennes. 

5. Généralisation. Le modèle de coexistence de populations du 
type prédateur — proie est l’un des modèles œcologiques élémen- 
taires appelés modèles de Voltaire-Lotka. Le schéma de la figu- 
re 4.70 permet d'étudier un autre modèle de Voltaire-Lotka : 


Ti —= dyiTs — doll, Gi, A2 > O0, 


(4.49) 


T3 = — GyTiZo + Ale, ess O2 > 0. 


I1 suffit de changer dans le schéma de la figure 4.70 le signe de la 
variable appliquée à l'entrée du circuit de commande préprogram- 
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ZT 
Z (t) 
Ty (£) 
0 A 
Z 
0 ZL 


Fig. 4.71 EE — 
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mée. Le système (4.49) décrit la coexistence de deux populations 
concurrentes du type prédateur — prédateur. La disparition de l’une 
des populations se traduit par un accroissement proportionnel du 
nombre de sujets de l’autre population. Les sujets meurent à une 
vitesse proportionnelle à la fréquence des rencontres des animaux 
des diverses populations. Le système (4.49) possède le même point 
singulier que le système (4.48) à la seule différence que ce point 
correspond désormais à un équilibre instable. En effet, en linéari- 
sant (4.49) au voisinage du point singulier, on obtient l'équation 
Uj — dyjlooli, 
d'où il vient 
u, = Achot, u, = B sh wt, 


et la trajectoire est une hyperbole 


u? uà 
AT — Br = 1 


Cela veut dire que le point singulier est une selle. Le câblage de (4.48) 
montre que quels que soient &;;, js, &oi et a, et quelles que soient 
les conditions initiales, une population de prédateurs extermine 
l’autre. 

Le système (4.49) peut à son tour être généralisé par l'adjonction 
du terme —a;;ri à la première équation et du terme —a,;r à la 
deuxième, termes qui traduisent l'effet de facteurs conditionnés 
par le surpeuplement de l'espèce. Le système (4.49) devient alors 

Ti = QiiTi — GyoTiTe — GisT Ajts Gy2s dy > 0, (4.50) 


: 3 
Lo = —- oyTiLo À Apolo — loss pts Uoos A3 > 0. 


Le câblage de (4.50) sous la condition que &@;s@o > Gi3as 
conduit au même résultat que le système (4.49). Pourtant, si aies, 
<< dysdes, le système (4.50) est stable. Cette stabilité traduit la 
coexistence de deux populations du type prédateur — prédateur 
dans un même milieu. 


Représentation suggestive 
des résultats de la simulation mathématique 


1. Projection d’un système de coordonnées spatiales sur un plan. 
La figure A.1 donne une représentation suggestive d'un système de 
coordonnées cartésiennes spatiales et d'un point P défini par ses 
coordonnées. D'une façon générale, la représentation d’un objet 
sur le plan sera dite suggestive si l'observateur s'’imagine facilement 
cet objet dans l’espace. La représentation suggestive s'obtient par 
projection. Celle-ci consiste à mener mentalement des droites appe- 
lées projetantes par les points de cet objet jusqu'à leur intersection 


Fig. A.1 


avec un plan dit de projection. Dans notre cas le plan de projection 
est le plan rapporté au système de coordonnées cartésiennes Y,0Y. 
On projettera une figure géométrique avec le système de coordonnées 
par rapport auquel elle est définie. Les projetantes seront parallèles. 

Elles percent le plan de projection sous un angle tel que Oz;, 
Oz, et Oz, se projettent respectivement sur OY,, OY. et sur une 
droite faisant un angle & avec OY:. Cette projection ne conserve 
pas tous les angles et toutes les distances. Ainsi sur la figure A. 
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l'angle z.0z, n'est pas droit, il vaut x — &; les projetantes de 
P sur les plans de coordonnées du système z,, z., z; sont les arêtes 
d’un parallélépipède rectangle. La projection de ces dernières sur 
Y,O0Y, modifie la longueur des arêtes parallèles à Oz, et les angles 
de toutes les faces à l’exception de celles qui sont parallèles au 
plan de projection. Ces modifications sont indispensables. Elles 
donnent du relief à la projection du parallélépipède sur le plan. 
En effet, c’est l’image que perçoit un observateur regardant un petit 
parallélépipède placé assez loin devant lui, au-dessous du niveau 
des yeux. 

Les appareils d'enregistrement modernes: traceurs de courbes, 
autoscripteurs, potentiomètres enregistreurs etc., traduisent l’infor- 
mation reçue de la machine par un dessin, schéma ou graphe dans 
un plan. Les instruments d'observation visuelle: display, indica- 
teurs électroniques utilisent des tubes cathodiques. La représenta- 
tion d'objets avec une machine implique l'établissement d'une 
relation entre les coordonnées cartésiennes spatiales z,, Z+, z, d'un 
point de cet espace, ses projections sur les plans z,02:, 2,023, 2:023 
et ses coordonnées cartésiennes correspondantes sur le plan Y,0Y:. 

2. Relations fondamentales dans l’axonométrie parallèle. Dans 
le cas le plus général de projection parallèle, les relations entre 
les coordonnées spatiales et leurs projections sur le plan peuvent 
être représentées sous forme de la transformation affine 


Yi= Gi0 + Qy121 + Gi920 + G323, 
V2 = G20 + 2121 + 2222 + Gos2ss 


où les coefficients a;; (i = 1, 2; j = 0, 1, 2, 3) dépendent de la 
disposition respective des systèmes de coordonnées spatial et plan 
ainsi que des déformations linéaires. Cependant l'effet de relief 
n'est pas donné par toute transformation affine. On se propose 
d'étudier analytiquement les principaux cas de représentation sug- 
gestive utilisés en dessin industriel et en géométrie descriptive 
où ils sont exposés dans un esprit essentiellement géométrique et 
réalisés par des constructions géométriques. 

La perspective axonométrique correspond au cas représenté sur 
la figure A.1. Elle est définie par la transformation affine 


Vi=zitaiszs, Vo 22 + Gos2s. 


L'égalité a,; = a+ = 1 traduit l'absence de déformations linéaires 
sur les axes Oz, et Oz.. Les constantes a;; et a, s'expriment sans 
peine en fonction de l'angle & et du coefficient y de déformation 
linéaire sur l’axe Oz,. En effet, si l’on considère que la figure A. 
est plane, on obtient a;; = —"7y sin &, a, = —"ycos &. Dans la 
perspective axonométrique « classique » on prend & = x/4 et y — 1/2 
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Fig. A.2 


Table de 
commutation 


ce qui donne 


2 2 
Yi=zi— 2; Yan, 
Si y — V2/2, on obtient des expressions extrêmement simples 
Yi=u—+, Yr=n—+ (A.1) 


facilement simulables sur la machine. 

L'isométrie orthogonale est caractérisée par une même déformation 
sur les trois axes. Par ailleurs l'axe Oz, est confondu avec OY,, les 
axes Oz, et Oz; forment avec OY, des angles de +2x/3. L'isométrie 
orthogonale donne les meilleurs résultats dans la représentation 
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de surfaces de révolution. On montre qu'elle est définie par la trans- 
formation affine 


Yi=n— (ze + 2)2, Ye = (2e — 23) V3/2. 


La dymétrie orthogonale se caractérise par l'absence de déforma- 
tions sur les axes Oz,, Oz, et une déformation linéaire sur l'axe Oz, 
égale à 0,5. D'autre part &« Æ 7°10° et Lz.0z, = Lz,0z, = 131°25. 
La dymétrie orthogonale est définie par la transformation affine 


Vi = —(2 WT +28, Ye = 3 (22 — 23 V7)/8. 


Les représentations combinées sont la superposition de plusieurs 
représentations d'un même objet. Sur la figure A.2 on a quatre repré- 
sentations d’un point : la représentation suggestive IV et trois projec- 
tions ÎTT, IT, I. La représentation suggestive IV doit être prise pour 
fondamentale dans la mesure où les trois autres en découlent. En ef- 
fet, si IV est donnée par (A.1), la projection Ï découle de (A.1) par 
la substitution z, = 0, la projection II par la substitution z3 = — Î 
et la projection III par la substitution z, — —/. Les quatre repré- 
sentations sont obtenues à l’aide de deux sommateurs montés comme 
indiqué sur la figure A.3. Le montage de schémas analogues pour 
l’isométrie et la dymétrie ne soulève aucune difficulté. 
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ÉLÉMENTS DE SIMULATION 
SUR CALCULATEURS ANALOGIQUES 


par À. Ourmaev 


Alexandre Ourmacv, candidat au doctorat ès sciences techniques, 
examine la programmation des calculateurs analogiques pour la 
résolution de problèmes techniques et scientifiques. Cet ouvrage ex- 
pose les principes de fonctionnement des calculateurs analogiques et 
des principaux blocs opérationnels ; on y trouve les méthodes géné- 
rales de programmation des calculateurs analogiques, ainsi que les 
méthodes de résolution de classes concrètes de problèmes embrassant 
plusieurs domaines: mécanique générale, physique, électronique, 
chimie physique, résistance des matériaux, programmation linéaire, 
etc. Ce manuel est illustré d’un grand nombre d’exemples ainsi que 
de problèmes dont la résolution permettra au lecteur d'asseoir ses 
connaissances théoriques. 


Cet ouvrage s'adresse aux étudiants des instituts 'techniques. 
Il peut être utilisé avec le plus grand profit par les professeurs dans la 
préparation de leurs cours sur les applications des calculateurs analogi- 
ques. Il rendra également le plus grand service aux ingénieurs et te- 
chaiciens qui utilisent les calculateurs analogiques. 


MÉTHODES NUMÉRIQUES DANS LES PROBLÈMES 
EXTRÉMAUX 


par B. Pchénitchny, You. Daniline 


Cet ouvrage expose les méthodes et les algorithmes utilisés en 
programmation mathématique, en économie, en théorie de commande 
optimale, etc., partout où l’on rencontre des problèmes exigeant le 
calcul numérique des extrémums des fonctionnelles. L'attention par- 
ticulière est portée aux algorithmes à convergence rapide qu’on peut 
utiliser avec succès sur les ordinateurs. Dans le livre sont examinées 
les méthodes de minimisation des fonctions, celles qui imposent les 
restrictions aux variables indépendantes aussi bien que celles qui ne 
tiennent pas compte de ces restrictions. 


Cette monographie s'adresse aux spécialistes en programmation 
mathématique, analyse numérique et théorie de commande optimale 
et à tous ceux qui se heurtent dans leur pratique aux problèmes de mini- 
misation des fonctions. 


ÉQUATIONS DIFFÉRENTIELLES 
AUX DÉRIVÉES PARTIELLES 


par V. Mikhaïlov 


Le manuel de V. Mikhaïlov, docteur ès sciences physiques et 
mathématiques, développe un cours fait à l'Institut physico-techni- 
que de Moscou. L'objet en est, d’abord, les problèmes aux limites 
usuels pour les équations elliptiques et, ensuite, le problème de 
Cauchy et les problèmes mixtes relatifs aux équations paraboliques et 
hyperboliques du second ordre. L'auteur s'inspire de la notion de 
solution généralisée, ce qui lui permet de traiter les équations à coef- 
ficients variables d'une manière aussi claire que les relations simples 
(équation de Poisson, équation d'onde, équation de la chaleur). On 
établit l'existence et l’unicité des solutions et on s'arrête sur les 
méthodes approchées pour les problèmes aux limites correspondants. 


Le livre s'adresse aux étudiants du 3° cycle en mathématiques et 
en physique. 


ÉQUATIONS INTÉGRALES 


par M. Krasnov, A. Kissélev, G. Makarenko 


Le présent Cours est destiné aux étudiants des grandes écoles ayant 
une formation mathématique poussée et à tous ceux qui s'intéressent 
aux méthodes de résolution des principaux types d'équations inté- 
grales. On donne au début de chaque chapitre un résumé des prin- 
cipaux résultats, des formules et l’on y étudie en détail des exemples 


types. 


